----------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  H:\RESEARCH\PriceIsRight\Publication\masterlog.log
  log type:  text
 opened on:  18 Aug 2014, 14:38:21

. do master.do

. cd H:\RESEARCH\PriceIsRight\Publication
H:\RESEARCH\PriceIsRight\Publication

. 
. ** This file generates all table and figure output for "The Price is Right: Updating of Inflation Expectations in a Randomized Price Information Experiment"
. ** Note, this file can be run in two versions by changing the local called "version," below
. ** Midpoint version: SPF perception gaps is coded as distance to median SPF forecast
. ** Interval version: SPF perception gaps is coded as minimum distance to 90% interval SPF forecast
. ** The former results are reported in the paper; interval version results are a robustness check (see footnote 15)
. 
. local version "midpoint"

. *local version "interval"
. 
. *************************
. ** PART ONE: DATA PREP **
. *************************
. 
. ** Data prep for density forecasts
. use densityforecasts, clear

. reshape wide  mean median var p25 p75, i(id question) j(wave) string
(note: j = a b)

Data                               long   ->   wide
-----------------------------------------------------------------------------
Number of obs.                     1402   ->     713
Number of variables                   8   ->      12
j variable (2 values)              wave   ->   (dropped)
xij variables:
                                   mean   ->   meana meanb
                                 median   ->   mediana medianb
                                    var   ->   vara varb
                                    p25   ->   p25a p25b
                                    p75   ->   p75a p75b
-----------------------------------------------------------------------------

. rename id prim_key

. save Output\densityforecasts_reshaped, replace
file Output\densityforecasts_reshaped.dta saved

. 
. ** Data prep for financial literacy questions
. use finlit_raw, clear 

. gen fin_q1 = 0 

. replace fin_q1 = 1 if F1 == 2 
(524 real changes made)

. gen fin_q2 = 0 

. replace fin_q2 = 1 if F2 == 2 
(521 real changes made)

. gen num_q1 = 0 

. replace num_q1 = 1 if N2_100 == 10 
(696 real changes made)

. gen num_q2 = 0 

. replace num_q2 = 1 if N2_1000 == 100 
(677 real changes made)

. gen num_q3 = 0 

. replace num_q3 = 1 if N3 == 500 
(532 real changes made)

. gen num_q4 = 0 

. replace num_q4 = 1 if N4_1year == 110 
(671 real changes made)

. gen num_q5 = 0 

. replace num_q5 = 1 if N4_2years == 121 
(488 real changes made)

. gen highfinlit = 0 

. replace highfinlit = 1 if fin_q1 + fin_q2 + num_q1 + num_q2 + num_q3 + num_q4 + num_q5 >= 7
(225 real changes made)

. keep prim_key highfinlit

. save Output\finlit_scores, replace
file Output\finlit_scores.dta saved

. 
. ***********************************
. ** PART TWO: MERGES AND CLEANING **
. ***********************************
. 
. ** Merge density forecast and finlit data with point forecast raw data:
. use pointforecast_raw, clear

. set more off

. 
. ** Note 10 respondents finished the survey but do not have density forecasts:
. merge 1:1 prim_key using Output\densityforecasts_reshaped

    Result                           # of obs.
    -----------------------------------------
    not matched                            26
        from master                        22  (_merge==1)
        from using                          4  (_merge==2)

    matched                               709  (_merge==3)
    -----------------------------------------

. drop _merge

. 
. ** Some other respondents didn't finish the survey at all:
. drop if tsend==""
(30 observations deleted)

. 
. ** Merge in finlit data:
. merge 1:1 prim_key using Output\finlit_scores

    Result                           # of obs.
    -----------------------------------------
    not matched                            88
        from master                        11  (_merge==1)
        from using                         77  (_merge==2)

    matched                               694  (_merge==3)
    -----------------------------------------

. drop if _merge==2
(77 observations deleted)

. drop _merge

. 
. ** Flag the treatment groups:
. gen all=1

. gen treated = .
(705 missing values generated)

. replace treated = 1 if randomscreen2!=4  & randomscreen2!=.
(541 real changes made)

. replace treated = 0 if randomscreen2==4  & randomscreen2!=.
(164 real changes made)

. gen control = !treated

. gen treatment = ""
(705 missing values generated)

. replace treatment = "food" if randomab==1
treatment was str1 now str4
(343 real changes made)

. replace treatment = "spf" if randomab==2
(362 real changes made)

. gen food=.
(705 missing values generated)

. replace food=0 if treatment!="food"
(362 real changes made)

. replace food=1 if treatment=="food"
(343 real changes made)

. gen spf=.
(705 missing values generated)

. replace spf=0 if treatment!="spf"
(343 real changes made)

. replace spf=1 if treatment=="spf"
(362 real changes made)

. gen question_type = ""
(705 missing values generated)

. replace question_type = "PP" if randomstage1==1
question_type was str1 now str2
(338 real changes made)

. replace question_type = "RI" if randomstage1==2
(367 real changes made)

. tab question question_type, mi

           |     question_type
  question |        PP         RI |     Total
-----------+----------------------+----------
           |         5          5 |        10 
        PG |       333          0 |       333 
        RI |         0        362 |       362 
-----------+----------------------+----------
     Total |       338        367 |       705 


. drop question

. rename question_type question

. gen pp=.
(705 missing values generated)

. replace pp=0 if question=="RI"
(367 real changes made)

. replace pp=1 if question=="PP"
(338 real changes made)

. gen ri=.
(705 missing values generated)

. replace ri=0 if question=="PP"
(338 real changes made)

. replace ri=1 if question=="RI"
(367 real changes made)

. gen ppXfood = pp * food

. gen riXfood = ri * food

. gen ppXtreated = pp * treated

. gen riXtreated = ri * treated

. gen controlXfood = control * food

. gen ppXfoodXtreated = pp * food * treated

. gen riXfoodXtreated = ri * food * treated

. gen treatment_string=""
(705 missing values generated)

. replace treatment_string = "2 PP SPF Control" if randomscreen2==4 & question=="PP" & treatment=="spf"
treatment_string was str1 now str16
(38 real changes made)

. replace treatment_string = "4 PP FoodBev Control" if randomscreen2==4 & question=="PP" & treatment=="food"
treatment_string was str16 now str20
(36 real changes made)

. replace treatment_string = "1 RI SPF Control" if randomscreen2==4 & question=="RI" & treatment=="spf"
(42 real changes made)

. replace treatment_string = "3 RI FoodBev Control" if randomscreen2==4 & question=="RI" & treatment=="food"
(48 real changes made)

. replace treatment_string = "6 PP SPF Info Treatment" if randomscreen2!=4 & question=="PP" & treatment=="spf"
treatment_string was str20 now str23
(139 real changes made)

. replace treatment_string = "8 PP FoodBev Info Treatment" if randomscreen2!=4 & question=="PP" & treatment=="food"
treatment_string was str23 now str27
(125 real changes made)

. replace treatment_string = "5 RI SPF Info Treatment" if randomscreen2!=4 & question=="RI" & treatment=="spf"
(143 real changes made)

. replace treatment_string = "7 RI FoodBev Info Treatment" if randomscreen2!=4 & question=="RI" & treatment=="food"
(134 real changes made)

. 
. ** Recode education:
. recode highesteducation (4=8) (6=10) (7=11) (8=11) (9=12) (10=13) (11=14) (12=14) (13=16) (14=18) (15=18) (16=20), generate(educ_years)
(705 differences between highesteducation and educ_years)

. sum educ_years, d

       RECODE of highesteducation (HIGHEST EDUCATION)
-------------------------------------------------------------
      Percentiles      Smallest
 1%           11              8
 5%           12              8
10%           12             10       Obs                 705
25%           13             10       Sum of Wgt.         705

50%           16                      Mean           14.99007
                        Largest       Std. Dev.      2.289716
75%           16             20
90%           18             20       Variance       5.242799
95%           18             20       Skewness       .1422045
99%           20             20       Kurtosis       2.138992

. 
. ** Find initial beliefs:
. gen belief_food = .
(705 missing values generated)

. replace belief_food = e1a_higher if e1a_lower==. | (e1a_lower==0 & e1a_higher!=.)
(309 real changes made)

. replace belief_food = -e1a_lower if e1a_higher==. | (e1a_higher==0 & e1a_lower!=.)
(0 real changes made)

. gen belief_spf = .
(705 missing values generated)

. replace belief_spf = e1b_inflation if e1b_deflation==. | (e1b_deflation==0 & e1b_inflation!=.)
(346 real changes made)

. replace belief_spf = -e1b_deflation if e1b_inflation==. | (e1b_inflation==0 & e1b_deflation!=.)
(13 real changes made)

. 
. ** Calculate perception gaps:
. ** (Note randomscreen2==4 --> control group. No treatment respondent has a zero perception gap, so we can code control respondents as zero perception gaps w/o affec
> ting the regressions below)
. gen gap_food = 1.39 - belief_food if randomscreen2!=4
(472 missing values generated)

. replace gap_food = 0 if randomscreen2==4 & food==1
(84 real changes made)

. gen gap_spf = 1.96 - belief_spf if randomscreen2!=4
(425 missing values generated)

. replace gap_spf = 0 if randomscreen2==4 & food==0
(80 real changes made)

. gen gap = .
(705 missing values generated)

. replace gap = gap_food if gap_food!=.
(317 real changes made)

. replace gap = gap_spf if gap_spf!=.
(360 real changes made)

. drop gap_spf gap_food

. gen abs_gap = abs(gap)
(28 missing values generated)

. gen log_info_belief = .
(705 missing values generated)

. replace log_info_belief = ln(1.39/belief_food) if log_info_belief==.
(308 real changes made)

. replace log_info_belief = ln(1.96/belief_spf) if log_info_belief==.
(345 real changes made)

. gen recoded_neg_belief = .1 if belief_food<=0 | belief_spf<=0
(690 missing values generated)

. replace log_info_belief = ln(1.39/recoded_neg_belief) if belief_food<=0
(1 real change made)

. replace log_info_belief = ln(1.96/recoded_neg_belief) if belief_spf<=0
(14 real changes made)

. if "`version'"=="interval" {
.         replace gap = 0 if spf & (belief_spf >= 1.19 & belief_spf<=3.03) & control==0
.         replace gap = 3.03 - belief_spf if spf & belief_spf>3.03 & belief_spf!=. & control==0
.         replace gap = 1.19 - belief_spf if spf & belief_spf<1.19 & control==0
.         }

.         
. ** Error xi term:
. gen gapXpp = pp * gap
(28 missing values generated)

. gen gapXri = ri * gap
(28 missing values generated)

. 
. ** Demographics:
. gen fem = .
(705 missing values generated)

. replace fem = 0 if gender==1
(299 real changes made)

. replace fem = 1 if gender==2
(406 real changes made)

. gen incover75 = familyincome>=14 if familyincome!=.
(2 missing values generated)

. gen college = highesteducation>=13 if highesteducation!=.

. gen age55andup = rage>=55 if rage!=. 

. gen male = fem==0 if fem!=.

. gen incunder75=incover75==0 if incover75!=.
(2 missing values generated)

. gen nocollege = college==0 if college!=.

. gen lofinlit = highfinlit==0 if highfinlit!=.
(11 missing values generated)

. gen ageunder55 = !age55andup

. 
. ** SPF/FoodBev Affected:
. gen spf_affected = e5b>=5 & e5b!=. if e5b!=.
(424 missing values generated)

. gen food_affected = e5a>=5 & e5a!=. if e5a!=.
(446 missing values generated)

. gen affected = . 
(705 missing values generated)

. replace affected = food_affected if affected==.
(259 real changes made)

. replace affected = spf_affected if affected==.
(281 real changes made)

. replace affected = 0 if affected==.
(165 real changes made)

. drop spf_affected food_affected

. gen noaffected = affected==0 if affected!=.

. 
. ** Affected xi term:
. gen affectedXpp = affected * pp

. gen affectedXri = affected * ri

. gen affectedXgapXpp = affected * gap * pp
(28 missing values generated)

. gen affectedXgapXri = affected * gap * ri
(28 missing values generated)

. gen affectedXtreated = affected * treated

. gen noaffectedXpp = noaffected * pp

. gen noaffectedXri = noaffected * ri

. gen noaffectedXgapXpp = noaffected * gap * pp
(28 missing values generated)

. gen noaffectedXgapXri = noaffected * gap * ri
(28 missing values generated)

. 
. ** Baseline point forecasts:
. gen base_pf_PP1 = .
(705 missing values generated)

. replace base_pf_PP1 =  a1_12months_up if (a1_12months_down==. | a1_12months_down==0)
(329 real changes made)

. replace base_pf_PP1 = a1_12months_down if (a1_12months_up==. | a1_12months_up==0)
(6 real changes made, 2 to missing)

. gen base_pf_RI1 = .
(705 missing values generated)

. replace base_pf_RI1 =  a3_12months_inflation if (a3_12months_deflation==. | a3_12months_deflation==0)
(342 real changes made)

. replace base_pf_RI1 = a3_12months_deflation if (a3_12months_inflation==. | a3_12months_inflation==0)
(20 real changes made, 7 to missing)

. gen base_pf_PP3 = .
(705 missing values generated)

. replace base_pf_PP3 =  a1_3years_up if (a1_3years_down==. | a1_3years_down==0)
(330 real changes made)

. replace base_pf_PP3 = a1_3years_down if (a1_3years_up==. | a1_3years_up==0)
(4 real changes made, 1 to missing)

. gen base_pf_RI3 = .
(705 missing values generated)

. replace base_pf_RI3 =  a3_3years_inflation if (a3_3years_deflation==. | a3_3years_deflation==0)
(339 real changes made)

. replace base_pf_RI3 = a3_3years_deflation if (a3_3years_inflation==. | a3_3years_inflation==0)
(18 real changes made, 4 to missing)

. 
. ** Recode a  probable typo:
. replace base_pf_RI1 = 20 if base_pf_RI1==2020
(1 real change made)

. 
. ** Revised point forecasts:
. gen end_pf_PP1 = .
(705 missing values generated)

. replace end_pf_PP1 =  e1_12months_up if (e1_12months_down==. | e1_12months_down==0)
(333 real changes made)

. replace end_pf_PP1 = e1_12months_down if (e1_12months_up==. | e1_12months_up==0)
(2 real changes made)

. gen end_pf_RI1 = .
(705 missing values generated)

. replace end_pf_RI1 =  e3_12months_inflation if (e3_12months_deflation==. | e3_12months_deflation==0)
(344 real changes made)

. replace end_pf_RI1 = e3_12months_deflation if (e3_12months_inflation==. | e3_12months_inflation==0)
(18 real changes made, 3 to missing)

. gen end_pf_PP3 = .
(705 missing values generated)

. replace end_pf_PP3 =  e1_3years_up if (e1_3years_down==. | e1_3years_down==0)
(331 real changes made)

. replace end_pf_PP3 = e1_3years_down if (e1_3years_up==. | e1_3years_up==0)
(4 real changes made)

. gen end_pf_RI3 = .
(705 missing values generated)

. replace end_pf_RI3 =  e3_3years_inflation if (e3_3years_deflation==. | e3_3years_deflation==0)
(337 real changes made)

. replace end_pf_RI3 = e3_3years_deflation if (e3_3years_inflation==. | e3_3years_inflation==0)
(23 real changes made, 1 to missing)

. 
. ** Point forecast revisions:
. gen pf_revision1 = end_pf_PP1 - base_pf_PP1
(375 missing values generated)

. replace pf_revision1 = end_pf_RI1 - base_pf_RI1 if pf_revision1==.
(343 real changes made)

. gen pf_revision3 = end_pf_PP3 - base_pf_PP3
(374 missing values generated)

. replace pf_revision3 = end_pf_RI3 - base_pf_RI3 if pf_revision3==.
(345 real changes made)

. 
. ** Density forecast revisions:
. gen dens_revision = meanb - meana
(27 missing values generated)

. 
. ** Absolute revision vars:
. foreach rev in pf_revision1 pf_revision3 dens_revision {
  2.         gen abs_`rev' = abs(`rev')
  3.         }
(32 missing values generated)
(29 missing values generated)
(27 missing values generated)

.         
. ** Demographic xi terms:
. foreach var of varlist fem incover75 college highfinlit age55andup  {
  2.         gen `var'Xtreated = `var' * treated
  3.         gen `var'Xfood = `var' * food
  4.         gen `var'XgapXpp = `var' * gap * pp
  5.         gen `var'XgapXri = `var' * gap * ri
  6.         gen `var'Xpp = `var' * pp
  7.         gen `var'Xri = `var' * ri
  8.         
.         }
(28 missing values generated)
(28 missing values generated)
(2 missing values generated)
(2 missing values generated)
(30 missing values generated)
(30 missing values generated)
(2 missing values generated)
(2 missing values generated)
(28 missing values generated)
(28 missing values generated)
(11 missing values generated)
(11 missing values generated)
(39 missing values generated)
(39 missing values generated)
(11 missing values generated)
(11 missing values generated)
(28 missing values generated)
(28 missing values generated)

. foreach var of varlist male incunder75 nocollege lofinlit ageunder55 {
  2.         gen `var'XgapXpp = `var' * gap * pp
  3.         gen `var'XgapXri = `var' * gap * ri
  4.         }
(28 missing values generated)
(28 missing values generated)
(30 missing values generated)
(30 missing values generated)
(28 missing values generated)
(28 missing values generated)
(39 missing values generated)
(39 missing values generated)
(28 missing values generated)
(28 missing values generated)

. 
. ** Drop anyone with missing PF revision or missing perception gap:
. drop if pf_revision1==. | riXtreated==. | gapXri==.
(52 observations deleted)

. drop if pf_revision1==. | ppXtreated==. | gapXpp==.     
(0 observations deleted)

. 
. ** Drop outliers by initial point forecast:
. foreach pf in PP1 RI1 PP3 RI3 {
  2.         drop if abs(end_pf_`pf')>50 & end_pf_`pf'!=.
  3.         drop if abs(base_pf_`pf')>50 & base_pf_`pf'!=.
  4.         }
(0 observations deleted)
(0 observations deleted)
(1 observation deleted)
(0 observations deleted)
(1 observation deleted)
(2 observations deleted)
(0 observations deleted)
(3 observations deleted)

. 
. ** Drop outliers by info beliefs:
. keep if (abs(belief_spf)<=50 | abs(belief_food)<=50) 
(11 observations deleted)

. 
. ** Drop PP respondents (see working paper for results with PP respondents):
. keep if ri
(310 observations deleted)

. 
. ** Unertainty xi terms:
. sum vara, d

                            a var
-------------------------------------------------------------
      Percentiles      Smallest
 1%     .2351484       .2351484
 5%     .3333333       .2351484
10%     .3564723       .2351484       Obs                 319
25%     1.157058       .2351484       Sum of Wgt.         319

50%            3                      Mean           18.51402
                        Largest       Std. Dev.      41.82337
75%     13.98883       230.6946
90%     59.18045       249.0606       Variance       1749.194
95%     83.76061       288.9399       Skewness       4.065468
99%     230.6946       295.2195       Kurtosis        22.0195

. gen unc = vara>=r(p50) if vara!=.
(6 missing values generated)

. gen uncXtreated = unc * treated
(6 missing values generated)

. gen uncXgapXri = gapXri * unc
(6 missing values generated)

. gen uncXgapXpp = gapXpp * unc
(6 missing values generated)

. gen uncXfemXgapXri = unc * femXgapXri
(6 missing values generated)

. gen uncXincover75XgapXri = unc * incover75XgapXri
(7 missing values generated)

. gen uncXcollegeXgapXri = unc * collegeXgapXri
(6 missing values generated)

. gen uncXhighfinlitXgapXri = unc * highfinlitXgapXri
(7 missing values generated)

. gen uncXage55andupXgapXri = unc * age55andupXgapXri
(6 missing values generated)

. gen uncXfem = unc * fem
(6 missing values generated)

. gen uncXincover75 = unc * incover75
(7 missing values generated)

. gen uncXcollege = unc * college
(6 missing values generated)

. gen uncXhighfinlit = unc * highfinlit
(7 missing values generated)

. gen uncXage55andup = unc * age55andup
(6 missing values generated)

. gen uncXmaleXgapXri = unc * maleXgapXri
(6 missing values generated)

. gen uncXincunder75XgapXri = unc * incunder75XgapXri
(7 missing values generated)

. gen uncXnocollegeXgapXri = unc * nocollegeXgapXri
(6 missing values generated)

. gen uncXlofinlitXgapXri = unc * lofinlitXgapXri
(7 missing values generated)

. gen uncXageunder55XgapXri = unc * ageunder55XgapXri
(6 missing values generated)

. 
. **SAVE:
. save Output\mastercleaned, replace
file Output\mastercleaned.dta saved

. 
. 
. *********************************************************
. ** PART THREE: REGRESSION TABLES / SUMMARY STAT TABLES **
. *********************************************************
. 
.         * * * * * *
.         * Table 1 *
.         * * * * * *
. 
. preserve

. 
. ** Recode control-group perception gaps from zero to real values:
. replace gap = 1.96 - belief_spf if randomscreen2==4 & spf==1
(40 real changes made)

. replace gap = 1.39 - belief_food if randomscreen2==4 & food==1
(41 real changes made)

. 
. gen pf1_non_rev = pf_revision1==0 if pf_revision1!=.

. 
. if "`version'"=="interval" {
.         replace gap = 0 if spf & (belief_spf >= 1.19 & belief_spf<=3.03) /* & control==0 */
.         replace gap = 3.03 - belief_spf if spf & belief_spf>3.03 & belief_spf!=. /* & control==0 */
.         replace gap = 1.19 - belief_spf if spf & belief_spf<1.19 /* & control==0 */
.         }

. 
. table treatment_string, contents(mean base_pf_RI1 mean pf_revision1 mean abs_pf_revision1)

----------------------------------------------------------------------------
           treatment_string | mean(base_p~1)  mean(pf_rev~1)  mean(abs_pf~1)
----------------------------+-----------------------------------------------
           1 RI SPF Control |         5.4185           -.931           2.019
       3 RI FoodBev Control |       4.768293       -.7317073        1.146341
    5 RI SPF Info Treatment |       5.417754        -1.96058        2.443188
7 RI FoodBev Info Treatment |       6.138207       -1.500472        2.584434
----------------------------------------------------------------------------

. table treatment_string, contents(median base_pf_RI1 median pf_revision1 median abs_pf_revision1)

-------------------------------------------------------------------------
           treatment_string | med(base_p~1)  med(pf_rev~1)  med(abs_pf~1)
----------------------------+--------------------------------------------
           1 RI SPF Control |             3              0            .38
       3 RI FoodBev Control |             3              0              0
    5 RI SPF Info Treatment |             3              0              1
7 RI FoodBev Info Treatment |             3              0              0
-------------------------------------------------------------------------

. table treatment_string, contents(mean gap median gap)

----------------------------------------------------
           treatment_string |  mean(gap)    med(gap)
----------------------------+-----------------------
           1 RI SPF Control |      -3.84       -2.04
       3 RI FoodBev Control |  -5.268537       -3.61
    5 RI SPF Info Treatment |  -3.612464       -1.04
7 RI FoodBev Info Treatment |  -7.147736       -3.61
----------------------------------------------------

. table treatment_string, contents(mean pf1_non_rev)

--------------------------------------------
           treatment_string | mean(pf1_no~v)
----------------------------+---------------
           1 RI SPF Control |           .475
       3 RI FoodBev Control |       .6341463
    5 RI SPF Info Treatment |       .4202898
7 RI FoodBev Info Treatment |       .5283019
--------------------------------------------

. table treatment_string if pf1_non_rev==0, contents(mean gap median gap)

----------------------------------------------------
           treatment_string |  mean(gap)    med(gap)
----------------------------+-----------------------
           1 RI SPF Control |  -2.754286       -2.04
       3 RI FoodBev Control |  -4.943333       -3.61
    5 RI SPF Info Treatment |     -4.815       -2.04
7 RI FoodBev Info Treatment |      -7.65       -5.11
----------------------------------------------------

. table treatment_string if pf1_non_rev==1, contents(mean gap median gap)

----------------------------------------------------
           treatment_string |  mean(gap)    med(gap)
----------------------------+-----------------------
           1 RI SPF Control |      -5.04       -1.04
       3 RI FoodBev Control |  -5.456154       -3.61
    5 RI SPF Info Treatment |  -1.953793       -1.04
7 RI FoodBev Info Treatment |  -6.699286       -3.61
----------------------------------------------------

. table treatment_string, contents(n pf_revision1)

-----------------------------------------
           treatment_string | N(pf_rev~1)
----------------------------+------------
           1 RI SPF Control |          40
       3 RI FoodBev Control |          41
    5 RI SPF Info Treatment |         138
7 RI FoodBev Info Treatment |         106
-----------------------------------------

. 
. ttest base_pf_RI1 if treatment=="spf" & question=="RI", by(treated)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      40      5.4185    1.086674    6.872729    3.220495    7.616505
       1 |     138    5.417754    .6271334    7.367149    4.177641    6.657867
---------+--------------------------------------------------------------------
combined |     178    5.417921    .5426576    7.239955     4.34701    6.488833
---------+--------------------------------------------------------------------
    diff |            .0007464    1.303787               -2.572322    2.573815
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   0.0006
Ho: diff = 0                                     degrees of freedom =      176

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.5002         Pr(|T| > |t|) = 0.9995          Pr(T > t) = 0.4998

. ttest base_pf_RI1 if treatment=="food" & question=="RI", by(treated)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      41    4.768293     .801123     5.12969    3.149163    6.387423
       1 |     106    6.138208    .6960896    7.166681    4.757991    7.518425
---------+--------------------------------------------------------------------
combined |     147    5.756122     .550369     6.67287    4.668403    6.843842
---------+--------------------------------------------------------------------
    diff |           -1.369915     1.22619               -3.793429    1.053599
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -1.1172
Ho: diff = 0                                     degrees of freedom =      145

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1329         Pr(|T| > |t|) = 0.2658          Pr(T > t) = 0.8671

. median base_pf_RI1 if treatment=="spf" & question=="RI", by(treated)

Median test

   Greater |
  than the |        treated
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        21         80 |       101 
       yes |        19         58 |        77 
-----------+----------------------+----------
     Total |        40        138 |       178 

          Pearson chi2(1) =   0.3782   Pr = 0.539

   Continuity corrected:
          Pearson chi2(1) =   0.1881   Pr = 0.664

. median base_pf_RI1 if treatment=="food" & question=="RI", by(treated)

Median test

   Greater |
  than the |        treated
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        24         56 |        80 
       yes |        17         50 |        67 
-----------+----------------------+----------
     Total |        41        106 |       147 

          Pearson chi2(1) =   0.3881   Pr = 0.533

   Continuity corrected:
          Pearson chi2(1) =   0.1922   Pr = 0.661

. 
. ttest pf_revision1 if treatment=="spf" & question=="RI", by(treated)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      40       -.931    .7228625    4.571784   -2.393127    .5311275
       1 |     138    -1.96058     .480019    5.638947   -2.909784   -1.011375
---------+--------------------------------------------------------------------
combined |     178   -1.729213    .4064292    5.422441   -2.531284   -.9271429
---------+--------------------------------------------------------------------
    diff |             1.02958    .9733963               -.8914513    2.950611
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   1.0577
Ho: diff = 0                                     degrees of freedom =      176

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.8542         Pr(|T| > |t|) = 0.2916          Pr(T > t) = 0.1458

. ttest pf_revision1 if treatment=="food" & question=="RI", by(treated)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      41   -.7317073    .5253041    3.363587   -1.793386    .3299718
       1 |     106   -1.500472    .5614329    5.780305   -2.613689    -.387254
---------+--------------------------------------------------------------------
combined |     147   -1.286054    .4305388    5.220005   -2.136948    -.435161
---------+--------------------------------------------------------------------
    diff |            .7687644    .9612168               -1.131042     2.66857
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   0.7998
Ho: diff = 0                                     degrees of freedom =      145

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.7874         Pr(|T| > |t|) = 0.4251          Pr(T > t) = 0.2126

. median pf_revision1 if treatment=="spf" & question=="RI", by(treated)

Median test

   Greater |
  than the |        treated
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        29        110 |       139 
       yes |        11         28 |        39 
-----------+----------------------+----------
     Total |        40        138 |       178 

          Pearson chi2(1) =   0.9423   Pr = 0.332

   Continuity corrected:
          Pearson chi2(1) =   0.5680   Pr = 0.451

. median pf_revision1 if treatment=="food" & question=="RI", by(treated)

Median test

   Greater |
  than the |        treated
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        35         90 |       125 
       yes |         6         16 |        22 
-----------+----------------------+----------
     Total |        41        106 |       147 

          Pearson chi2(1) =   0.0049   Pr = 0.944

   Continuity corrected:
          Pearson chi2(1) =   0.0352   Pr = 0.851

. 
. ttest abs_pf_revision1 if treatment=="spf" & question=="RI", by(treated)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      40       2.019    .6635004    4.196345    .6769437    3.361056
       1 |     138    2.443188    .4635787    5.445816    1.526493    3.359883
---------+--------------------------------------------------------------------
combined |     178    2.347865    .3885028    5.183274    1.581171    3.114559
---------+--------------------------------------------------------------------
    diff |           -.4241884    .9328675               -2.265235    1.416858
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.4547
Ho: diff = 0                                     degrees of freedom =      176

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.3249         Pr(|T| > |t|) = 0.6499          Pr(T > t) = 0.6751

. ttest abs_pf_revision1 if treatment=="food" & question=="RI", by(treated)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      41    1.146341    .5064354    3.242769    .1227974    2.169886
       1 |     106    2.584434     .522529    5.379765    1.548355    3.620512
---------+--------------------------------------------------------------------
combined |     147    2.183333    .4050205    4.910613    1.382873    2.983794
---------+--------------------------------------------------------------------
    diff |           -1.438092    .8983336               -3.213613    .3374276
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -1.6008
Ho: diff = 0                                     degrees of freedom =      145

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0558         Pr(|T| > |t|) = 0.1116          Pr(T > t) = 0.9442

. median abs_pf_revision1 if treatment=="spf" & question=="RI", by(treated)

Median test

   Greater |
  than the |        treated
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        23         66 |        89 
       yes |        17         72 |        89 
-----------+----------------------+----------
     Total |        40        138 |       178 

          Pearson chi2(1) =   1.1609   Pr = 0.281

   Continuity corrected:
          Pearson chi2(1) =   0.8062   Pr = 0.369

. median abs_pf_revision1 if treatment=="food" & question=="RI", by(treated)

Median test

   Greater |
  than the |        treated
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        26         56 |        82 
       yes |        15         50 |        65 
-----------+----------------------+----------
     Total |        41        106 |       147 

          Pearson chi2(1) =   1.3428   Pr = 0.247

   Continuity corrected:
          Pearson chi2(1) =   0.9480   Pr = 0.330

. 
. tab pf1_non_rev treated if treatment=="spf" & question=="RI", chi2

pf1_non_re |        treated
         v |         0          1 |     Total
-----------+----------------------+----------
         0 |        21         80 |       101 
         1 |        19         58 |        77 
-----------+----------------------+----------
     Total |        40        138 |       178 

          Pearson chi2(1) =   0.3782   Pr = 0.539

. tab pf1_non_rev treated if treatment=="food" & question=="RI", chi2

pf1_non_re |        treated
         v |         0          1 |     Total
-----------+----------------------+----------
         0 |        15         50 |        65 
         1 |        26         56 |        82 
-----------+----------------------+----------
     Total |        41        106 |       147 

          Pearson chi2(1) =   1.3428   Pr = 0.247

. 
. ttest gap if treatment=="spf" & question=="RI", by(treated)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      40       -3.84    1.093219    6.914125   -6.051245   -1.628755
       1 |     138   -3.612464     .601304    7.063722   -4.801501   -2.423427
---------+--------------------------------------------------------------------
combined |     178   -3.663596     .525542    7.011605    -4.70073   -2.626461
---------+--------------------------------------------------------------------
    diff |           -.2275362    1.262549               -2.719219    2.264147
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.1802
Ho: diff = 0                                     degrees of freedom =      176

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.4286         Pr(|T| > |t|) = 0.8572          Pr(T > t) = 0.5714

. ttest gap if treatment=="food" & question=="RI", by(treated)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      41   -5.268537    .7144926    4.574985    -6.71258   -3.824493
       1 |     106   -7.147736     .691519    7.119624    -8.51889   -5.776581
---------+--------------------------------------------------------------------
combined |     147   -6.623605    .5402439    6.550109   -7.691314   -5.555897
---------+--------------------------------------------------------------------
    diff |            1.879199    1.198685               -.4899525    4.248351
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   1.5677
Ho: diff = 0                                     degrees of freedom =      145

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9404         Pr(|T| > |t|) = 0.1191          Pr(T > t) = 0.0596

. median gap if treatment=="spf" & question=="RI", by(treated)

Median test

   Greater |
  than the |        treated
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        28         94 |       122 
       yes |        12         44 |        56 
-----------+----------------------+----------
     Total |        40        138 |       178 

          Pearson chi2(1) =   0.0511   Pr = 0.821

   Continuity corrected:
          Pearson chi2(1) =   0.0011   Pr = 0.974

. median gap if treatment=="food" & question=="RI", by(treated)

Median test

   Greater |
  than the |        treated
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        29         71 |       100 
       yes |        12         35 |        47 
-----------+----------------------+----------
     Total |        41        106 |       147 

          Pearson chi2(1) =   0.1912   Pr = 0.662

   Continuity corrected:
          Pearson chi2(1) =   0.0576   Pr = 0.810

. 
. ttest gap if treated==0 & treatment=="spf" & question=="RI", by(pf1_non_rev)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      21   -2.754286    .9515473    4.360537   -4.739178   -.7693929
       1 |      19       -5.04    2.046249    8.919392   -9.339009   -.7409905
---------+--------------------------------------------------------------------
combined |      40       -3.84    1.093219    6.914125   -6.051245   -1.628755
---------+--------------------------------------------------------------------
    diff |            2.285714    2.186579               -2.140783    6.712211
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   1.0453
Ho: diff = 0                                     degrees of freedom =       38

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.8488         Pr(|T| > |t|) = 0.3025          Pr(T > t) = 0.1512

. median gap if treated==0 & treatment=="spf" & question=="RI", by(pf1_non_rev)

Median test

   Greater |
  than the |      pf1_non_rev
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        13          9 |        22 
       yes |         8         10 |        18 
-----------+----------------------+----------
     Total |        21         19 |        40 

          Pearson chi2(1) =   0.8516   Pr = 0.356

   Continuity corrected:
          Pearson chi2(1) =   0.3656   Pr = 0.545

. ttest gap if treated==0 & treatment=="food" & question=="RI", by(pf1_non_rev)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      15   -4.943333     1.27864    4.952152   -7.685744   -2.200923
       1 |      26   -5.456154     .869469    4.433439   -7.246859   -3.665449
---------+--------------------------------------------------------------------
combined |      41   -5.268537    .7144926    4.574985    -6.71258   -3.824493
---------+--------------------------------------------------------------------
    diff |            .5128204    1.500021               -2.521257    3.546898
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   0.3419
Ho: diff = 0                                     degrees of freedom =       39

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.6329         Pr(|T| > |t|) = 0.7343          Pr(T > t) = 0.3671

. median gap if treated==0 & treatment=="food" & question=="RI", by(pf1_non_rev)

Median test

   Greater |
  than the |      pf1_non_rev
    median |         0          1 |     Total
-----------+----------------------+----------
        no |         9         20 |        29 
       yes |         6          6 |        12 
-----------+----------------------+----------
     Total |        15         26 |        41 

          Pearson chi2(1) =   1.3159   Pr = 0.251

   Continuity corrected:
          Pearson chi2(1) =   0.6254   Pr = 0.429

. ttest gap if treated==1 & treatment=="spf" & question=="RI", by(pf1_non_rev)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      80      -4.815    .9627517    8.611113   -6.731309   -2.898691
       1 |      58   -1.953793    .4609917    3.510808   -2.876913   -1.030673
---------+--------------------------------------------------------------------
combined |     138   -3.612464     .601304    7.063722   -4.801501   -2.423427
---------+--------------------------------------------------------------------
    diff |           -2.861207    1.197789               -5.229907    -.492507
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -2.3887
Ho: diff = 0                                     degrees of freedom =      136

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0091         Pr(|T| > |t|) = 0.0183          Pr(T > t) = 0.9909

. median gap if treated==1 & treatment=="spf" & question=="RI", by(pf1_non_rev)

Median test

   Greater |
  than the |      pf1_non_rev
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        60         34 |        94 
       yes |        20         24 |        44 
-----------+----------------------+----------
     Total |        80         58 |       138 

          Pearson chi2(1) =   4.1534   Pr = 0.042

   Continuity corrected:
          Pearson chi2(1) =   3.4335   Pr = 0.064

. ttest gap if treated==1 & treatment=="food" & question=="RI", by(pf1_non_rev)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      50       -7.65    .9462968    6.691329   -9.551655   -5.748345
       1 |      56   -6.699286    1.003978    7.513084   -8.711303   -4.687269
---------+--------------------------------------------------------------------
combined |     106   -7.147736     .691519    7.119624    -8.51889   -5.776581
---------+--------------------------------------------------------------------
    diff |           -.9507143    1.388777               -3.704712    1.803283
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.6846
Ho: diff = 0                                     degrees of freedom =      104

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.2476         Pr(|T| > |t|) = 0.4951          Pr(T > t) = 0.7524

. median gap if treated==1 & treatment=="food" & question=="RI", by(pf1_non_rev)

Median test

   Greater |
  than the |      pf1_non_rev
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        35         36 |        71 
       yes |        15         20 |        35 
-----------+----------------------+----------
     Total |        50         56 |       106 

          Pearson chi2(1) =   0.3900   Pr = 0.532

   Continuity corrected:
          Pearson chi2(1) =   0.1744   Pr = 0.676

. 
. restore

. 
.         * * * * * *
.         * Table 2 *
.         * * * * * *
. 
. eststo clear

. foreach rev in pf_revision1 dens_revision pf_revision3 {
  2.         eststo: reg `rev' riXtreated gapXri if ri==1 & treatment=="food", robust
  3.         eststo: reg `rev' riXtreated gapXri if ri==1 & treatment=="spf", robust
  4.         }

Linear regression                                      Number of obs =     147
                                                       F(  2,   144) =    0.71
                                                       Prob > F      =  0.4924
                                                       R-squared     =  0.0052
                                                       Root MSE      =  5.2425

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |  -.9431847   .8159627    -1.16   0.250    -2.555996    .6696268
      gapXri |  -.0244022   .0897932    -0.27   0.786    -.2018852    .1530809
       _cons |  -.7317073   .5242353    -1.40   0.165    -1.767898    .3044831
------------------------------------------------------------------------------
(est1 stored)

Linear regression                                      Number of obs =     178
                                                       F(  2,   175) =    4.67
                                                       Prob > F      =  0.0106
                                                       R-squared     =  0.2092
                                                       Root MSE      =  4.8496

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |   .3900892   .8268053     0.47   0.638    -1.241704    2.021883
      gapXri |   .3929919   .1295381     3.03   0.003     .1373339    .6486498
       _cons |      -.931   .7198616    -1.29   0.198    -2.351728    .4897278
------------------------------------------------------------------------------
(est2 stored)

Linear regression                                      Number of obs =     141
                                                       F(  2,   138) =    0.39
                                                       Prob > F      =  0.6777
                                                       R-squared     =  0.0046
                                                       Root MSE      =  2.6143

------------------------------------------------------------------------------
             |               Robust
dens_revis~n |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |  -.2593106   .3600098    -0.72   0.473    -.9711592    .4525381
      gapXri |  -.0289042   .0385546    -0.75   0.455    -.1051383    .0473298
       _cons |  -.2404254   .1979921    -1.21   0.227    -.6319159    .1510651
------------------------------------------------------------------------------
(est3 stored)

Linear regression                                      Number of obs =     173
                                                       F(  2,   170) =    4.74
                                                       Prob > F      =  0.0100
                                                       R-squared     =  0.0724
                                                       Root MSE      =  2.2279

------------------------------------------------------------------------------
             |               Robust
dens_revis~n |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |  -.3454365    .385639    -0.90   0.372    -1.106694    .4158214
      gapXri |   .0955135   .0338527     2.82   0.005     .0286878    .1623393
       _cons |  -.3698361   .3362661    -1.10   0.273    -1.033631    .2939588
------------------------------------------------------------------------------
(est4 stored)

Linear regression                                      Number of obs =     145
                                                       F(  2,   142) =    2.47
                                                       Prob > F      =  0.0878
                                                       R-squared     =  0.0246
                                                       Root MSE      =  6.1261

------------------------------------------------------------------------------
             |               Robust
pf_revision3 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |  -1.513075   1.008812    -1.50   0.136    -3.507307    .4811564
      gapXri |   .0647199   .0855881     0.76   0.451    -.1044716    .2339113
       _cons |   .0243902   .6052097     0.04   0.968    -1.171995    1.220775
------------------------------------------------------------------------------
(est5 stored)

Linear regression                                      Number of obs =     176
                                                       F(  2,   173) =    0.67
                                                       Prob > F      =  0.5116
                                                       R-squared     =  0.0384
                                                       Root MSE      =   5.732

------------------------------------------------------------------------------
             |               Robust
pf_revision3 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |   .1811103   1.013627     0.18   0.858    -1.819557    2.181777
      gapXri |   .1797399    .155918     1.15   0.251    -.1280067    .4874864
       _cons |  -1.030769   .8965947    -1.15   0.252    -2.800442    .7389037
------------------------------------------------------------------------------
(est6 stored)

. eststo: reg pf_revision1 riXtreated gapXri uncXgapXri unc uncXtreated if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     144
                                                       F(  5,   138) =    3.00
                                                       Prob > F      =  0.0133
                                                       R-squared     =  0.0833
                                                       Root MSE      =  5.1367

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |   .1680353   .2970327     0.57   0.573    -.4192885    .7553591
      gapXri |  -.0168326   .0143174    -1.18   0.242    -.0451424    .0114773
  uncXgapXri |  -.1285807   .1441926    -0.89   0.374    -.4136933    .1565319
         unc |  -1.466667   1.372041    -1.07   0.287    -4.179608    1.246275
 uncXtreated |  -2.747804   1.965867    -1.40   0.164    -6.634921    1.139312
       _cons |        -.2   .2600334    -0.77   0.443     -.714165     .314165
------------------------------------------------------------------------------
(est7 stored)

. eststo: reg pf_revision1 riXtreated gapXri uncXgapXri unc uncXtreated if ri==1 & treatment=="spf", robust               

Linear regression                                      Number of obs =     175
                                                       F(  5,   169) =   27.03
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.3296
                                                       Root MSE      =  4.5388

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |  -.1630022   .1435934    -1.14   0.258    -.4464699    .1204656
      gapXri |   .1005366   .0100045    10.05   0.000     .0807868    .1202864
  uncXgapXri |   .4522307   .1446106     3.13   0.002     .1667549    .7377065
         unc |  -1.890667   1.140412    -1.66   0.099    -4.141955    .3606213
 uncXtreated |   .7323151   1.393426     0.53   0.600    -2.018447    3.483077
       _cons |   .2506667   .0928958     2.70   0.008     .0672811    .4340522
------------------------------------------------------------------------------
(est8 stored)

. test gapXri = uncXgapXri

 ( 1)  gapXri - uncXgapXri = 0

       F(  1,   169) =    5.83
            Prob > F =    0.0168

. esttab using Output\Table2.csv, se replace star(* 0.10 ** 0.05 *** 0.01) r2
(output written to Output\Table2.csv)

.         
.         * * * * * *
.         * Table 3 *
.         * * * * * *
. 
. ** Note this table is limited to non-controls except in the "control" column
. 
. preserve

. 
. gen expec_in_tail = .
(325 missing values generated)

. replace expec_in_tail = 0 if (base_pf_PP1>=0 & base_pf_PP1<5 & pp) | (base_pf_RI1>=0 & base_pf_RI1<5 & ri)
(204 real changes made)

. replace expec_in_tail = 1 if (base_pf_PP1>=5 & pp) | (base_pf_RI1>=5 & ri)
(121 real changes made)

. 
. ** Recode control-group perception gaps from zero to real values:
. replace gap = 1.39 - belief_food if randomscreen2==4 & food==1
(41 real changes made)

. replace gap = 1.96 - belief_spf if randomscreen2==4 & spf==1
(40 real changes made)

. 
. if "`version'"=="interval" {
.         replace gap = 0 if spf & (belief_spf >= 1.19 & belief_spf<=3.03) /* & control==0 */
.         replace gap = 3.03 - belief_spf if spf & belief_spf>3.03 & belief_spf!=. /* & control==0 */
.         replace gap = 1.19 - belief_spf if spf & belief_spf<1.19 /* & control==0 */
.         }

. 
. gen pf1_non_rev = pf_revision1==0 if pf_revision1!=.

. rename pf1_non_rev nonrevise

. 
. table control if question=="RI", contents(n base_pf_RI1 mean base_pf_RI1 median base_pf_RI1)

----------------------------------------------------------
  control |    N(base_p~1)  mean(base_p~1)   med(base_p~1)
----------+-----------------------------------------------
        0 |            244        5.730738               3
        1 |             81        5.089383               3
----------------------------------------------------------

. table fem if question=="RI" & control==0, contents(n base_pf_RI1 mean base_pf_RI1 median base_pf_RI1)

----------------------------------------------------------
      fem |    N(base_p~1)  mean(base_p~1)   med(base_p~1)
----------+-----------------------------------------------
        0 |            107        4.358411               3
        1 |            137        6.802555               3
----------------------------------------------------------

. table age55andup if question=="RI" & control==0, contents(n base_pf_RI1 mean base_pf_RI1 median base_pf_RI1)

----------------------------------------------------------
age55andu |
p         |    N(base_p~1)  mean(base_p~1)   med(base_p~1)
----------+-----------------------------------------------
        0 |            130           5.665               3
        1 |            114        5.805702               3
----------------------------------------------------------

. table highfinlit if question=="RI" & control==0, contents(n base_pf_RI1 mean base_pf_RI1 median base_pf_RI1)

----------------------------------------------------------
highfinli |
t         |    N(base_p~1)  mean(base_p~1)   med(base_p~1)
----------+-----------------------------------------------
        0 |            159        6.999371             3.5
        1 |             84        3.367857            2.65
----------------------------------------------------------

. table college if question=="RI" & control==0, contents(n base_pf_RI1 mean base_pf_RI1 median base_pf_RI1)

----------------------------------------------------------
  college |    N(base_p~1)  mean(base_p~1)   med(base_p~1)
----------+-----------------------------------------------
        0 |            114        7.169737               4
        1 |            130        4.468846             2.8
----------------------------------------------------------

. table incover75 if question=="RI" & control==0, contents(n base_pf_RI1 mean base_pf_RI1 median base_pf_RI1)

----------------------------------------------------------
incover75 |    N(base_p~1)  mean(base_p~1)   med(base_p~1)
----------+-----------------------------------------------
        0 |            150           6.411               3
        1 |             94        4.645213             2.5
----------------------------------------------------------

. 
. table control if question=="RI", contents(mean expec_in_tail)

--------------------------
  control | mean(expec_~l)
----------+---------------
        0 |       .3729508
        1 |       .3703704
--------------------------

. table fem if question=="RI" & control==0, contents(mean expec_in_tail)

--------------------------
      fem | mean(expec_~l)
----------+---------------
        0 |       .2803738
        1 |       .4452555
--------------------------

. table age55andup if question=="RI" & control==0, contents(mean expec_in_tail)

--------------------------
age55andu |
p         | mean(expec_~l)
----------+---------------
        0 |       .3692308
        1 |        .377193
--------------------------

. table highfinlit if question=="RI" & control==0, contents(mean expec_in_tail)

--------------------------
highfinli |
t         | mean(expec_~l)
----------+---------------
        0 |       .4528302
        1 |       .2261905
--------------------------

. table college if question=="RI" & control==0, contents(mean expec_in_tail)

--------------------------
  college | mean(expec_~l)
----------+---------------
        0 |       .4736842
        1 |       .2846154
--------------------------

. table incover75 if question=="RI" & control==0, contents(mean expec_in_tail)

--------------------------
incover75 | mean(expec_~l)
----------+---------------
        0 |       .4333333
        1 |       .2765957
--------------------------

. 
. table control if food==1 & question=="RI", contents(mean gap sd gap)

----------------------------------
  control |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -7.147736    7.119624
        1 |  -5.268537    4.574985
----------------------------------

. table fem if food==1 & question=="RI" & control==0, contents(mean gap sd gap)

----------------------------------
      fem |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -7.196957    8.204541
        1 |      -7.11    6.234132
----------------------------------

. table age55andup if food==1 & question=="RI" & control==0, contents(mean gap sd gap)

----------------------------------
age55andu |
p         |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -6.571538    7.689977
        1 |  -7.702592    6.548134
----------------------------------

. table highfinlit if food==1 & question=="RI" & control==0, contents(mean gap sd gap)

----------------------------------
highfinli |
t         |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -8.212941    7.056572
        1 |  -5.285676    7.008146
----------------------------------

. table college if food==1 & question=="RI" & control==0, contents(mean gap sd gap)

----------------------------------
  college |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -8.930755    7.839186
        1 |  -5.364717    5.866809
----------------------------------

. table incover75 if food==1 & question=="RI" & control==0, contents(mean gap sd gap)

----------------------------------
incover75 |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -8.937869    8.100867
        1 |  -4.721111    4.583678
----------------------------------

. 
. table control if spf==1 & question=="RI", contents(mean gap sd gap)

----------------------------------
  control |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -3.612464    7.063722
        1 |      -3.84    6.914125
----------------------------------

. table fem if spf==1 & question=="RI" & control==0, contents(mean gap sd gap)

----------------------------------
      fem |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -1.597377    2.674102
        1 |  -5.208831     8.85765
----------------------------------

. table age55andup if spf==1 & question=="RI" & control==0, contents(mean gap sd gap)

----------------------------------
age55andu |
p         |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -3.770769     7.72359
        1 |  -3.406667    6.158821
----------------------------------

. table highfinlit if spf==1 & question=="RI" & control==0, contents(mean gap sd gap)

----------------------------------
highfinli |
t         |  mean(gap)     sd(gap)
----------+-----------------------
        0 |   -4.71033    8.336073
        1 |  -1.486809    2.375679
----------------------------------

. table college if spf==1 & question=="RI" & control==0, contents(mean gap sd gap)

----------------------------------
  college |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -5.220328    8.419241
        1 |  -2.338701    5.501351
----------------------------------

. table incover75 if spf==1 & question=="RI" & control==0, contents(mean gap sd gap)

----------------------------------
incover75 |  mean(gap)     sd(gap)
----------+-----------------------
        0 |  -4.231011    7.760442
        1 |   -2.48898    5.477458
----------------------------------

. 
. table control if question=="RI", contents(mean vara sd vara)

----------------------------------
  control | mean(vara)    sd(vara)
----------+-----------------------
        0 |   19.96501    45.02146
        1 |   14.17918    30.15313
----------------------------------

. table fem if question=="RI" & control==0, contents(mean vara sd vara)

----------------------------------
      fem | mean(vara)    sd(vara)
----------+-----------------------
        0 |   10.78067    30.77665
        1 |    27.1617    52.60889
----------------------------------

. table age55andup if question=="RI" & control==0, contents(mean vara sd vara)

----------------------------------
age55andu |
p         | mean(vara)    sd(vara)
----------+-----------------------
        0 |   18.55018    42.66286
        1 |   21.62423    47.78328
----------------------------------

. table highfinlit if question=="RI" & control==0, contents(mean vara sd vara)

----------------------------------
highfinli |
t         | mean(vara)    sd(vara)
----------+-----------------------
        0 |   26.48526    51.75003
        1 |    8.23765     25.5551
----------------------------------

. table college if question=="RI" & control==0, contents(mean vara sd vara)

----------------------------------
  college | mean(vara)    sd(vara)
----------+-----------------------
        0 |   26.25586    48.54396
        1 |   14.60072    41.22009
----------------------------------

. table incover75 if question=="RI" & control==0, contents(mean vara sd vara)

----------------------------------
incover75 | mean(vara)    sd(vara)
----------+-----------------------
        0 |   25.17642    53.61007
        1 |    11.6381    24.00288
----------------------------------

. 
. table control if question=="RI", contents(mean pf_revision1 median pf_revision1 sd pf_revision1)

----------------------------------------------------------
  control | mean(pf_rev~1)   med(pf_rev~1)    sd(pf_rev~1)
----------+-----------------------------------------------
        0 |      -1.760697               0        5.693557
        1 |      -.8301235               0        3.981991
----------------------------------------------------------

. table fem if question=="RI" & control==0, contents(mean pf_revision1 median pf_revision1 sd pf_revision1)

----------------------------------------------------------
      fem | mean(pf_rev~1)   med(pf_rev~1)    sd(pf_rev~1)
----------+-----------------------------------------------
        0 |      -.6985981               0        3.882316
        1 |      -2.590219               0        6.677766
----------------------------------------------------------

. table age55andup if question=="RI" & control==0, contents(mean pf_revision1 median pf_revision1 sd pf_revision1)

----------------------------------------------------------
age55andu |
p         | mean(pf_rev~1)   med(pf_rev~1)    sd(pf_rev~1)
----------+-----------------------------------------------
        0 |      -2.037385               0        6.024794
        1 |      -1.445175               0        5.299416
----------------------------------------------------------

. table highfinlit if question=="RI" & control==0, contents(mean pf_revision1 median pf_revision1 sd pf_revision1)

----------------------------------------------------------
highfinli |
t         | mean(pf_rev~1)   med(pf_rev~1)    sd(pf_rev~1)
----------+-----------------------------------------------
        0 |      -2.449749               0        6.732546
        1 |       -.477381               0        2.450982
----------------------------------------------------------

. table college if question=="RI" & control==0, contents(mean pf_revision1 median pf_revision1 sd pf_revision1)

----------------------------------------------------------
  college | mean(pf_rev~1)   med(pf_rev~1)    sd(pf_rev~1)
----------+-----------------------------------------------
        0 |      -2.176842               0         6.73464
        1 |      -1.395769               0        4.587665
----------------------------------------------------------

. table incover75 if question=="RI" & control==0, contents(mean pf_revision1 median pf_revision1 sd pf_revision1)

----------------------------------------------------------
incover75 | mean(pf_rev~1)   med(pf_rev~1)    sd(pf_rev~1)
----------+-----------------------------------------------
        0 |      -2.057733               0        6.463035
        1 |      -1.286702               0        4.172395
----------------------------------------------------------

. 
. table control if question=="RI", contents(mean end_pf_RI1 median end_pf_RI1)

------------------------------------------
  control | mean(end_pf~1)   med(end_pf~1)
----------+-------------------------------
        0 |       3.970041               2
        1 |       4.259259               3
------------------------------------------

. table fem if question=="RI" & control==0, contents(mean end_pf_RI1 median end_pf_RI1)

------------------------------------------
      fem | mean(end_pf~1)   med(end_pf~1)
----------+-------------------------------
        0 |       3.659813               2
        1 |       4.212336               2
------------------------------------------

. table age55andup if question=="RI" & control==0, contents(mean end_pf_RI1 median end_pf_RI1)

------------------------------------------
age55andu |
p         | mean(end_pf~1)   med(end_pf~1)
----------+-------------------------------
        0 |       3.627615               2
        1 |       4.360526             2.5
------------------------------------------

. table highfinlit if question=="RI" & control==0, contents(mean end_pf_RI1 median end_pf_RI1)

------------------------------------------
highfinli |
t         | mean(end_pf~1)   med(end_pf~1)
----------+-------------------------------
        0 |       4.549623               3
        1 |       2.890476               2
------------------------------------------

. table college if question=="RI" & control==0, contents(mean end_pf_RI1 median end_pf_RI1)

------------------------------------------
  college | mean(end_pf~1)   med(end_pf~1)
----------+-------------------------------
        0 |       4.992895               3
        1 |       3.073077               2
------------------------------------------

. table incover75 if question=="RI" & control==0, contents(mean end_pf_RI1 median end_pf_RI1)

------------------------------------------
incover75 | mean(end_pf~1)   med(end_pf~1)
----------+-------------------------------
        0 |       4.353267             2.1
        1 |       3.358511               2
------------------------------------------

. 
. table control if question=="RI", contents(mean nonrevise)

--------------------------
  control | mean(nonrev~e)
----------+---------------
        0 |       .4672131
        1 |       .5555556
--------------------------

. table fem if question=="RI" & control==0, contents(mean nonrevise)

--------------------------
      fem | mean(nonrev~e)
----------+---------------
        0 |       .5046729
        1 |       .4379562
--------------------------

. table age55andup if question=="RI" & control==0, contents(mean nonrevise)

--------------------------
age55andu |
p         | mean(nonrev~e)
----------+---------------
        0 |       .4615385
        1 |       .4736842
--------------------------

. table highfinlit if question=="RI" & control==0, contents(mean nonrevise)

--------------------------
highfinli |
t         | mean(nonrev~e)
----------+---------------
        0 |       .4213836
        1 |        .547619
--------------------------

. table college if question=="RI" & control==0, contents(mean nonrevise)

--------------------------
  college | mean(nonrev~e)
----------+---------------
        0 |       .4736842
        1 |       .4615385
--------------------------

. table incover75 if question=="RI" & control==0, contents(mean nonrevise)

--------------------------
incover75 | mean(nonrev~e)
----------+---------------
        0 |       .4466667
        1 |             .5
--------------------------

. 
. table control if question=="RI", contents(sd base_pf_RI1 sd end_pf_RI1)

--------------------------------------
  control | sd(base_p~1)  sd(end_pf~1)
----------+---------------------------
        0 |     7.274657      4.641906
        1 |      6.02417      4.131518
--------------------------------------

. table fem if question=="RI" & control==0, contents(sd base_pf_RI1 sd end_pf_RI1)

--------------------------------------
      fem | sd(base_p~1)  sd(end_pf~1)
----------+---------------------------
        0 |     5.134151      4.015297
        1 |     8.448247      5.079264
--------------------------------------

. table age55andup if question=="RI" & control==0, contents(sd base_pf_RI1 sd end_pf_RI1)

--------------------------------------
age55andu |
p         | sd(base_p~1)  sd(end_pf~1)
----------+---------------------------
        0 |     7.762309      3.947301
        1 |     6.708744      5.316024
--------------------------------------

. table highfinlit if question=="RI" & control==0, contents(sd base_pf_RI1 sd end_pf_RI1)

--------------------------------------
highfinli |
t         | sd(base_p~1)  sd(end_pf~1)
----------+---------------------------
        0 |     8.501119      5.489685
        1 |     2.916268      1.966381
--------------------------------------

. table college if question=="RI" & control==0, contents(sd base_pf_RI1 sd end_pf_RI1)

--------------------------------------
  college | sd(base_p~1)  sd(end_pf~1)
----------+---------------------------
        0 |     7.977085      5.546181
        1 |     6.364884      3.450914
--------------------------------------

. table incover75 if question=="RI" & control==0, contents(sd base_pf_RI1 sd end_pf_RI1)

--------------------------------------
incover75 | sd(base_p~1)  sd(end_pf~1)
----------+---------------------------
        0 |     7.973473      5.176807
        1 |     5.871995      3.570651
--------------------------------------

. 
. ttest base_pf_RI1 if question=="RI" & control==0, by(fem)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     107    4.358411    .4963371    5.134151    3.374375    5.342448
       1 |     137    6.802555    .7217825    8.448247    5.375186    8.229923
---------+--------------------------------------------------------------------
combined |     244    5.730738    .4657122    7.274657     4.81339    6.648086
---------+--------------------------------------------------------------------
    diff |           -2.444144    .9272659               -4.270686   -.6176011
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -2.6359
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0045         Pr(|T| > |t|) = 0.0089          Pr(T > t) = 0.9955

. ttest base_pf_RI1 if question=="RI" & control==0, by(age55andup)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     130       5.665    .6807995    7.762309    4.318021    7.011978
       1 |     114    5.805702    .6283314    6.708744    4.560864     7.05054
---------+--------------------------------------------------------------------
combined |     244    5.730738    .4657122    7.274657     4.81339    6.648086
---------+--------------------------------------------------------------------
    diff |           -.1407018    .9353163               -1.983102    1.701698
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.1504
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.4403         Pr(|T| > |t|) = 0.8805          Pr(T > t) = 0.5597

. ttest base_pf_RI1 if question=="RI" & control==0, by(highfinlit)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     159    6.999371    .6741826    8.501119    5.667798    8.330944
       1 |      84    3.367857    .3181909    2.916268    2.734988    4.000726
---------+--------------------------------------------------------------------
combined |     243    5.744033    .4674421    7.286701    4.823258    6.664807
---------+--------------------------------------------------------------------
    diff |            3.631514    .9567232                1.746907    5.516121
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   3.7958
Ho: diff = 0                                     degrees of freedom =      241

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9999         Pr(|T| > |t|) = 0.0002          Pr(T > t) = 0.0001

. ttest base_pf_RI1 if question=="RI" & control==0, by(college)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     114    7.169737    .7471224    7.977084    5.689553    8.649921
       1 |     130    4.468846    .5582373    6.364884     3.36436    5.573332
---------+--------------------------------------------------------------------
combined |     244    5.730738    .4657122    7.274657     4.81339    6.648086
---------+--------------------------------------------------------------------
    diff |            2.700891    .9191053                .8904232    4.511358
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   2.9386
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9982         Pr(|T| > |t|) = 0.0036          Pr(T > t) = 0.0018

. ttest base_pf_RI1 if question=="RI" & control==0, by(incover75)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     150       6.411    .6510313    7.973473    5.124553    7.697447
       1 |      94    4.645213      .60565    5.871995    3.442512    5.847914
---------+--------------------------------------------------------------------
combined |     244    5.730738    .4657122    7.274657     4.81339    6.648086
---------+--------------------------------------------------------------------
    diff |            1.765787    .9522025               -.1098757     3.64145
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   1.8544
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9676         Pr(|T| > |t|) = 0.0649          Pr(T > t) = 0.0324

. 
. median base_pf_RI1 if question=="RI" & control==0, by(fem)

Median test

   Greater |
  than the |          fem
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        66         70 |       136 
       yes |        41         67 |       108 
-----------+----------------------+----------
     Total |       107        137 |       244 

          Pearson chi2(1) =   2.7296   Pr = 0.099

   Continuity corrected:
          Pearson chi2(1) =   2.3174   Pr = 0.128

. median base_pf_RI1 if question=="RI" & control==0, by(age55andup)

Median test

   Greater |
  than the |      age55andup
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        73         63 |       136 
       yes |        57         51 |       108 
-----------+----------------------+----------
     Total |       130        114 |       244 

          Pearson chi2(1) =   0.0195   Pr = 0.889

   Continuity corrected:
          Pearson chi2(1) =   0.0001   Pr = 0.992

. median base_pf_RI1 if question=="RI" & control==0, by(highfinlit)

Median test

   Greater |
  than the |      highfinlit
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        78         57 |       135 
       yes |        81         27 |       108 
-----------+----------------------+----------
     Total |       159         84 |       243 

          Pearson chi2(1) =   7.8680   Pr = 0.005

   Continuity corrected:
          Pearson chi2(1) =   7.1250   Pr = 0.008

. median base_pf_RI1 if question=="RI" & control==0, by(college)

Median test

   Greater |
  than the |        college
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        51         85 |       136 
       yes |        63         45 |       108 
-----------+----------------------+----------
     Total |       114        130 |       244 

          Pearson chi2(1) =  10.4960   Pr = 0.001

   Continuity corrected:
          Pearson chi2(1) =   9.6757   Pr = 0.002

. median base_pf_RI1 if question=="RI" & control==0, by(incover75)

Median test

   Greater |
  than the |       incover75
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        76         60 |       136 
       yes |        74         34 |       108 
-----------+----------------------+----------
     Total |       150         94 |       244 

          Pearson chi2(1) =   4.0585   Pr = 0.044

   Continuity corrected:
          Pearson chi2(1) =   3.5425   Pr = 0.060

. 
. tab expec_in_tail fem if question=="RI" & control==0, chi2

expec_in_t |          fem
       ail |         0          1 |     Total
-----------+----------------------+----------
         0 |        77         76 |       153 
         1 |        30         61 |        91 
-----------+----------------------+----------
     Total |       107        137 |       244 

          Pearson chi2(1) =   6.9840   Pr = 0.008

. tab expec_in_tail age55andup if question=="RI" & control==0, chi2

expec_in_t |      age55andup
       ail |         0          1 |     Total
-----------+----------------------+----------
         0 |        82         71 |       153 
         1 |        48         43 |        91 
-----------+----------------------+----------
     Total |       130        114 |       244 

          Pearson chi2(1) =   0.0165   Pr = 0.898

. tab expec_in_tail highfinlit if question=="RI" & control==0, chi2

expec_in_t |      highfinlit
       ail |         0          1 |     Total
-----------+----------------------+----------
         0 |        87         65 |       152 
         1 |        72         19 |        91 
-----------+----------------------+----------
     Total |       159         84 |       243 

          Pearson chi2(1) =  12.0523   Pr = 0.001

. tab expec_in_tail college if question=="RI" & control==0, chi2

expec_in_t |        college
       ail |         0          1 |     Total
-----------+----------------------+----------
         0 |        60         93 |       153 
         1 |        54         37 |        91 
-----------+----------------------+----------
     Total |       114        130 |       244 

          Pearson chi2(1) =   9.2842   Pr = 0.002

. tab expec_in_tail incover75 if question=="RI" & control==0, chi2

expec_in_t |       incover75
       ail |         0          1 |     Total
-----------+----------------------+----------
         0 |        85         68 |       153 
         1 |        65         26 |        91 
-----------+----------------------+----------
     Total |       150         94 |       244 

          Pearson chi2(1) =   6.0705   Pr = 0.014

. 
. ttest gap if food==1 & question=="RI" & control==0, by(fem)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      46   -7.196956    1.209694    8.204541   -9.633404   -4.760509
       1 |      60       -7.11     .804823    6.234133   -8.720447   -5.499553
---------+--------------------------------------------------------------------
combined |     106   -7.147736     .691519    7.119624    -8.51889   -5.776581
---------+--------------------------------------------------------------------
    diff |           -.0869565    1.401927                -2.86703    2.693117
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.0620
Ho: diff = 0                                     degrees of freedom =      104

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.4753         Pr(|T| > |t|) = 0.9507          Pr(T > t) = 0.5247

. ttest gap if food==1 & question=="RI" & control==0, by(age55andup)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      52   -6.571538    1.066408    7.689977   -8.712442   -4.430635
       1 |      54   -7.702593    .8910882    6.548134   -9.489889   -5.915296
---------+--------------------------------------------------------------------
combined |     106   -7.147736     .691519    7.119624    -8.51889   -5.776581
---------+--------------------------------------------------------------------
    diff |            1.131054    1.385487               -1.616418    3.878526
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   0.8164
Ho: diff = 0                                     degrees of freedom =      104

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.7919         Pr(|T| > |t|) = 0.4162          Pr(T > t) = 0.2081

. ttest gap if food==1 & question=="RI" & control==0, by(highfinlit)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      68   -8.212941    .8557351    7.056573   -9.920996   -6.504887
       1 |      37   -5.285676    1.152132    7.008146   -7.622308   -2.949043
---------+--------------------------------------------------------------------
combined |     105   -7.181429    .6973074    7.145274   -8.564215   -5.798642
---------+--------------------------------------------------------------------
    diff |           -2.927266    1.438111                -5.77942   -.0751107
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -2.0355
Ho: diff = 0                                     degrees of freedom =      103

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0222         Pr(|T| > |t|) = 0.0444          Pr(T > t) = 0.9778

. ttest gap if food==1 & question=="RI" & control==0, by(college)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      53   -8.930755    1.076795    7.839186    -11.0915   -6.770007
       1 |      53   -5.364717    .8058682    5.866809    -6.98181   -3.747624
---------+--------------------------------------------------------------------
combined |     106   -7.147736     .691519    7.119624    -8.51889   -5.776581
---------+--------------------------------------------------------------------
    diff |           -3.566038    1.344958               -6.233139   -.8989363
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -2.6514
Ho: diff = 0                                     degrees of freedom =      104

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0046         Pr(|T| > |t|) = 0.0093          Pr(T > t) = 0.9954

. ttest gap if food==1 & question=="RI" & control==0, by(incover75)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      61   -8.937869     1.03721    8.100867    -11.0126    -6.86314
       1 |      45   -4.721111    .6832943    4.583678     -6.0982   -3.344022
---------+--------------------------------------------------------------------
combined |     106   -7.147736     .691519    7.119624    -8.51889   -5.776581
---------+--------------------------------------------------------------------
    diff |           -4.216758    1.343592               -6.881152   -1.552364
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -3.1384
Ho: diff = 0                                     degrees of freedom =      104

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0011         Pr(|T| > |t|) = 0.0022          Pr(T > t) = 0.9989

. 
. ttest gap if spf==1 & question=="RI" & control==0, by(fem)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      61   -1.597377    .3423836    2.674102   -2.282246   -.9125078
       1 |      77   -5.208831    1.009423     8.85765   -7.219271   -3.198391
---------+--------------------------------------------------------------------
combined |     138   -3.612464     .601304    7.063722   -4.801501   -2.423427
---------+--------------------------------------------------------------------
    diff |            3.611454    1.175096                1.287631    5.935277
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   3.0733
Ho: diff = 0                                     degrees of freedom =      136

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9987         Pr(|T| > |t|) = 0.0026          Pr(T > t) = 0.0013

. ttest gap if spf==1 & question=="RI" & control==0, by(age55andup)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      78   -3.770769    .8745244     7.72359    -5.51217   -2.029369
       1 |      60   -3.406667    .7951004    6.158821   -4.997659   -1.815675
---------+--------------------------------------------------------------------
combined |     138   -3.612464     .601304    7.063722   -4.801501   -2.423427
---------+--------------------------------------------------------------------
    diff |           -.3641025    1.217021               -2.770836    2.042631
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.2992
Ho: diff = 0                                     degrees of freedom =      136

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.3826         Pr(|T| > |t|) = 0.7653          Pr(T > t) = 0.6174

. ttest gap if spf==1 & question=="RI" & control==0, by(highfinlit)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      91    -4.71033    .8738579    8.336073   -6.446401   -2.974258
       1 |      47   -1.486808    .3465283    2.375678   -2.184334   -.7892832
---------+--------------------------------------------------------------------
combined |     138   -3.612464     .601304    7.063722   -4.801501   -2.423427
---------+--------------------------------------------------------------------
    diff |           -3.223521    1.243126               -5.681879   -.7651635
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -2.5931
Ho: diff = 0                                     degrees of freedom =      136

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0053         Pr(|T| > |t|) = 0.0106          Pr(T > t) = 0.9947

. ttest gap if spf==1 & question=="RI" & control==0, by(college)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      61   -5.220328    1.077973     8.41924   -7.376596    -3.06406
       1 |      77   -2.338701    .6269372    5.501351   -3.587355   -1.090048
---------+--------------------------------------------------------------------
combined |     138   -3.612464     .601304    7.063722   -4.801501   -2.423427
---------+--------------------------------------------------------------------
    diff |           -2.881627     1.18983               -5.234587   -.5286664
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -2.4219
Ho: diff = 0                                     degrees of freedom =      136

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0084         Pr(|T| > |t|) = 0.0168          Pr(T > t) = 0.9916

. ttest gap if spf==1 & question=="RI" & control==0, by(incover75)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      89   -4.231011    .8226053    7.760442   -5.865766   -2.596256
       1 |      49    -2.48898    .7824941    5.477459   -4.062289   -.9156698
---------+--------------------------------------------------------------------
combined |     138   -3.612464     .601304    7.063722   -4.801501   -2.423427
---------+--------------------------------------------------------------------
    diff |           -1.742032    1.252284                 -4.2185    .7344369
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -1.3911
Ho: diff = 0                                     degrees of freedom =      136

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0832         Pr(|T| > |t|) = 0.1665          Pr(T > t) = 0.9168

. 
. ttest vara if question=="RI" & control==0, by(fem)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     105    10.78067    3.003493    30.77665    4.824633    16.73671
       1 |     134     27.1617    4.544716     52.6089    18.17243    36.15097
---------+--------------------------------------------------------------------
combined |     239    19.96501    2.912196    45.02146    14.22804    25.70199
---------+--------------------------------------------------------------------
    diff |           -16.38103    5.783037               -27.77375   -4.988307
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -2.8326
Ho: diff = 0                                     degrees of freedom =      237

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0025         Pr(|T| > |t|) = 0.0050          Pr(T > t) = 0.9975

. ttest vara if question=="RI" & control==0, by(age55andup)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     129    18.55018    3.756256    42.66286    11.11778    25.98257
       1 |     110    21.62423    4.555956    47.78327    12.59448    30.65399
---------+--------------------------------------------------------------------
combined |     239    19.96501    2.912196    45.02146    14.22804    25.70199
---------+--------------------------------------------------------------------
    diff |           -3.074057    5.851792               -14.60223    8.454114
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.5253
Ho: diff = 0                                     degrees of freedom =      237

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.2999         Pr(|T| > |t|) = 0.5999          Pr(T > t) = 0.7001

. ttest vara if question=="RI" & control==0, by(highfinlit)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     154    26.48526    4.170136    51.75003    18.24678    34.72374
       1 |      84     8.23765     2.78829     25.5551    2.691853    13.78345
---------+--------------------------------------------------------------------
combined |     238    20.04493    2.923357    45.09935    14.28584    25.80401
---------+--------------------------------------------------------------------
    diff |            18.24761    6.014053                6.399522    30.09569
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   3.0342
Ho: diff = 0                                     degrees of freedom =      236

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9987         Pr(|T| > |t|) = 0.0027          Pr(T > t) = 0.0013

. ttest vara if question=="RI" & control==0, by(college)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     110    26.25586    4.628485    48.54396    17.08236    35.42937
       1 |     129    14.60072    3.629226    41.22008     7.41968    21.78177
---------+--------------------------------------------------------------------
combined |     239    19.96501    2.912196    45.02146    14.22804    25.70199
---------+--------------------------------------------------------------------
    diff |            11.65514    5.806046                .2170899    23.09319
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   2.0074
Ho: diff = 0                                     degrees of freedom =      237

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9771         Pr(|T| > |t|) = 0.0458          Pr(T > t) = 0.0229

. ttest vara if question=="RI" & control==0, by(incover75)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     147    25.17642    4.421684    53.61007    16.43764    33.91519
       1 |      92     11.6381    2.502473    24.00288    6.667243    16.60895
---------+--------------------------------------------------------------------
combined |     239    19.96501    2.912196    45.02146    14.22804    25.70199
---------+--------------------------------------------------------------------
    diff |            13.53832    5.932815                 1.85053    25.22611
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   2.2819
Ho: diff = 0                                     degrees of freedom =      237

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9883         Pr(|T| > |t|) = 0.0234          Pr(T > t) = 0.0117

. 
. ttest pf_revision1 if question=="RI" & control==0, by(fem)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     107   -.6985981    .3753177    3.882316   -1.442702    .0455057
       1 |     137   -2.590219      .57052    6.677766   -3.718457   -1.461981
---------+--------------------------------------------------------------------
combined |     244   -1.760697    .3644926    5.693557   -2.478665   -1.042728
---------+--------------------------------------------------------------------
    diff |            1.891621    .7259613                .4616112     3.32163
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   2.6057
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9951         Pr(|T| > |t|) = 0.0097          Pr(T > t) = 0.0049

. ttest pf_revision1 if question=="RI" & control==0, by(age55andup)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     130   -2.037385    .5284094    6.024794   -3.082856   -.9919136
       1 |     114   -1.445175    .4963357    5.299415   -2.428506   -.4618448
---------+--------------------------------------------------------------------
combined |     244   -1.760697    .3644926    5.693557   -2.478665   -1.042728
---------+--------------------------------------------------------------------
    diff |           -.5922092     .731075               -2.032292    .8478735
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.8101
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.2094         Pr(|T| > |t|) = 0.4187          Pr(T > t) = 0.7906

. ttest pf_revision1 if question=="RI" & control==0, by(highfinlit)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     159   -2.449748    .5339256    6.732546   -3.504301   -1.395196
       1 |      84    -.477381    .2674241    2.450982   -1.009277    .0545148
---------+--------------------------------------------------------------------
combined |     243   -1.767942    .3659234    5.704181   -2.488744   -1.047141
---------+--------------------------------------------------------------------
    diff |           -1.972367    .7604652               -3.470374   -.4743605
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -2.5936
Ho: diff = 0                                     degrees of freedom =      241

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0050         Pr(|T| > |t|) = 0.0101          Pr(T > t) = 0.9950

. ttest pf_revision1 if question=="RI" & control==0, by(college)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     114   -2.176842    .6307568     6.73464   -3.426485   -.9271991
       1 |     130   -1.395769    .4023648    4.587665   -2.191858   -.5996806
---------+--------------------------------------------------------------------
combined |     244   -1.760697    .3644926    5.693557   -2.478665   -1.042728
---------+--------------------------------------------------------------------
    diff |           -.7810729    .7303417               -2.219711    .6575652
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -1.0695
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1430         Pr(|T| > |t|) = 0.2859          Pr(T > t) = 0.8570

. ttest pf_revision1 if question=="RI" & control==0, by(incover75)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     150   -2.057733    .5277046    6.463035   -3.100485   -1.014982
       1 |      94   -1.286702    .4303497    4.172395   -2.141291   -.4321129
---------+--------------------------------------------------------------------
combined |     244   -1.760697    .3644926    5.693557   -2.478665   -1.042728
---------+--------------------------------------------------------------------
    diff |           -.7710312    .7488855               -2.246197    .7041348
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -1.0296
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1521         Pr(|T| > |t|) = 0.3042          Pr(T > t) = 0.8479

. 
. median pf_revision1 if question=="RI" & control==0, by(fem)

Median test

   Greater |
  than the |          fem
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        81        119 |       200 
       yes |        26         18 |        44 
-----------+----------------------+----------
     Total |       107        137 |       244 

          Pearson chi2(1) =   5.0626   Pr = 0.024

   Continuity corrected:
          Pearson chi2(1) =   4.3357   Pr = 0.037

. median pf_revision1 if question=="RI" & control==0, by(age55andup)

Median test

   Greater |
  than the |      age55andup
    median |         0          1 |     Total
-----------+----------------------+----------
        no |       107         93 |       200 
       yes |        23         21 |        44 
-----------+----------------------+----------
     Total |       130        114 |       244 

          Pearson chi2(1) =   0.0218   Pr = 0.883

   Continuity corrected:
          Pearson chi2(1) =   0.0004   Pr = 0.985

. median pf_revision1 if question=="RI" & control==0, by(highfinlit)

Median test

   Greater |
  than the |      highfinlit
    median |         0          1 |     Total
-----------+----------------------+----------
        no |       129         70 |       199 
       yes |        30         14 |        44 
-----------+----------------------+----------
     Total |       159         84 |       243 

          Pearson chi2(1) =   0.1796   Pr = 0.672

   Continuity corrected:
          Pearson chi2(1) =   0.0618   Pr = 0.804

. median pf_revision1 if question=="RI" & control==0, by(college)

Median test

   Greater |
  than the |        college
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        96        104 |       200 
       yes |        18         26 |        44 
-----------+----------------------+----------
     Total |       114        130 |       244 

          Pearson chi2(1) =   0.7285   Pr = 0.393

   Continuity corrected:
          Pearson chi2(1) =   0.4715   Pr = 0.492

. median pf_revision1 if question=="RI" & control==0, by(incover75)

Median test

   Greater |
  than the |       incover75
    median |         0          1 |     Total
-----------+----------------------+----------
        no |       124         76 |       200 
       yes |        26         18 |        44 
-----------+----------------------+----------
     Total |       150         94 |       244 

          Pearson chi2(1) =   0.1289   Pr = 0.720

   Continuity corrected:
          Pearson chi2(1) =   0.0353   Pr = 0.851

. 
. ttest end_pf_RI1 if question=="RI" & control==0, by(fem)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     107    3.659813    .3881735    4.015298    2.890221    4.429405
       1 |     137    4.212336    .4339508    5.079264    3.354172      5.0705
---------+--------------------------------------------------------------------
combined |     244    3.970041    .2971676    4.641906    3.384688    4.555394
---------+--------------------------------------------------------------------
    diff |           -.5525227    .5990631               -1.732566    .6275208
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.9223
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1786         Pr(|T| > |t|) = 0.3573          Pr(T > t) = 0.8214

. ttest end_pf_RI1 if question=="RI" & control==0, by(age55andup)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     130    3.627615    .3462012    3.947301    2.942648    4.312583
       1 |     114    4.360526    .4978913    5.316024    3.374114    5.346939
---------+--------------------------------------------------------------------
combined |     244    3.970041    .2971676    4.641906    3.384688    4.555394
---------+--------------------------------------------------------------------
    diff |           -.7329109     .594984               -1.904919    .4390976
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -1.2318
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1096         Pr(|T| > |t|) = 0.2192          Pr(T > t) = 0.8904

. ttest end_pf_RI1 if question=="RI" & control==0, by(highfinlit)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     159    4.549623    .4353603    5.489685    3.689746    5.409499
       1 |      84    2.890476    .2145498    1.966381    2.463745    3.317207
---------+--------------------------------------------------------------------
combined |     243    3.976091    .2983312    4.650523    3.388433    4.563748
---------+--------------------------------------------------------------------
    diff |            1.659146    .6194354                .4389478    2.879345
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   2.6785
Ho: diff = 0                                     degrees of freedom =      241

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9960         Pr(|T| > |t|) = 0.0079          Pr(T > t) = 0.0040

. ttest end_pf_RI1 if question=="RI" & control==0, by(college)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     114    4.992895    .5194475    5.546181    3.963776    6.022014
       1 |     130    3.073077    .3026651    3.450914    2.474247    3.671907
---------+--------------------------------------------------------------------
combined |     244    3.970041    .2971676    4.641906    3.384688    4.555394
---------+--------------------------------------------------------------------
    diff |            1.919818    .5839482                .7695479    3.070088
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   3.2877
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9994         Pr(|T| > |t|) = 0.0012          Pr(T > t) = 0.0006

. ttest end_pf_RI1 if question=="RI" & control==0, by(incover75)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     150    4.353267    .4226845    5.176807    3.518036    5.188497
       1 |      94    3.358511    .3682846    3.570651     2.62717    4.089851
---------+--------------------------------------------------------------------
combined |     244    3.970041    .2971676    4.641906    3.384688    4.555394
---------+--------------------------------------------------------------------
    diff |             .994756    .6085449               -.2039649    2.193477
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   1.6346
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9483         Pr(|T| > |t|) = 0.1034          Pr(T > t) = 0.0517

. 
. median end_pf_RI1 if question=="RI" & control==0, by(fem)

Median test

   Greater |
  than the |          fem
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        54         71 |       125 
       yes |        53         66 |       119 
-----------+----------------------+----------
     Total |       107        137 |       244 

          Pearson chi2(1) =   0.0443   Pr = 0.833

   Continuity corrected:
          Pearson chi2(1) =   0.0066   Pr = 0.935

. median end_pf_RI1 if question=="RI" & control==0, by(age55andup)

Median test

   Greater |
  than the |      age55andup
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        70         55 |       125 
       yes |        60         59 |       119 
-----------+----------------------+----------
     Total |       130        114 |       244 

          Pearson chi2(1) =   0.7625   Pr = 0.383

   Continuity corrected:
          Pearson chi2(1) =   0.5548   Pr = 0.456

. median end_pf_RI1 if question=="RI" & control==0, by(highfinlit)

Median test

   Greater |
  than the |      highfinlit
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        74         51 |       125 
       yes |        85         33 |       118 
-----------+----------------------+----------
     Total |       159         84 |       243 

          Pearson chi2(1) =   4.4202   Pr = 0.036

   Continuity corrected:
          Pearson chi2(1) =   3.8710   Pr = 0.049

. median end_pf_RI1 if question=="RI" & control==0, by(college)

Median test

   Greater |
  than the |        college
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        48         77 |       125 
       yes |        66         53 |       119 
-----------+----------------------+----------
     Total |       114        130 |       244 

          Pearson chi2(1) =   7.1296   Pr = 0.008

   Continuity corrected:
          Pearson chi2(1) =   6.4607   Pr = 0.011

. median end_pf_RI1 if question=="RI" & control==0, by(incover75)

Median test

   Greater |
  than the |       incover75
    median |         0          1 |     Total
-----------+----------------------+----------
        no |        75         50 |       125 
       yes |        75         44 |       119 
-----------+----------------------+----------
     Total |       150         94 |       244 

          Pearson chi2(1) =   0.2356   Pr = 0.627

   Continuity corrected:
          Pearson chi2(1) =   0.1252   Pr = 0.724

. 
. tab nonrevise fem if question=="RI" & control==0, chi2

           |          fem
 nonrevise |         0          1 |     Total
-----------+----------------------+----------
         0 |        53         77 |       130 
         1 |        54         60 |       114 
-----------+----------------------+----------
     Total |       107        137 |       244 

          Pearson chi2(1) =   1.0743   Pr = 0.300

. tab nonrevise age55andup if question=="RI" & control==0, chi2

           |      age55andup
 nonrevise |         0          1 |     Total
-----------+----------------------+----------
         0 |        70         60 |       130 
         1 |        60         54 |       114 
-----------+----------------------+----------
     Total |       130        114 |       244 

          Pearson chi2(1) =   0.0360   Pr = 0.850

. tab nonrevise highfinlit if question=="RI" & control==0, chi2

           |      highfinlit
 nonrevise |         0          1 |     Total
-----------+----------------------+----------
         0 |        92         38 |       130 
         1 |        67         46 |       113 
-----------+----------------------+----------
     Total |       159         84 |       243 

          Pearson chi2(1) =   3.5207   Pr = 0.061

. tab nonrevise college if question=="RI" & control==0, chi2

           |        college
 nonrevise |         0          1 |     Total
-----------+----------------------+----------
         0 |        60         70 |       130 
         1 |        54         60 |       114 
-----------+----------------------+----------
     Total |       114        130 |       244 

          Pearson chi2(1) =   0.0360   Pr = 0.850

. tab nonrevise incover75 if question=="RI" & control==0, chi2

           |       incover75
 nonrevise |         0          1 |     Total
-----------+----------------------+----------
         0 |        83         47 |       130 
         1 |        67         47 |       114 
-----------+----------------------+----------
     Total |       150         94 |       244 

          Pearson chi2(1) =   0.6603   Pr = 0.416

. 
. ksmirnov base_pf_RI1 if question=="RI" & control==0, by(fem) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1649    0.038
 1:                  0.0000    1.000
 Combined K-S:       0.1649    0.076      0.066

Note: Ties exist in combined dataset;
      there are 29 unique values out of 244 observations.

. ksmirnov base_pf_RI1 if question=="RI" & control==0, by(age55andup) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0501    0.737
 1:                 -0.0221    0.942
 Combined K-S:       0.0501    0.998      0.995

Note: Ties exist in combined dataset;
      there are 29 unique values out of 244 observations.

. ksmirnov base_pf_RI1 if question=="RI" & control==0, by(highfinlit) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0083    0.992
 1:                 -0.2266    0.004
 Combined K-S:       0.2266    0.007      0.006

Note: Ties exist in combined dataset;
      there are 29 unique values out of 243 observations.

. ksmirnov base_pf_RI1 if question=="RI" & control==0, by(college) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0197    0.954
 1:                 -0.2065    0.006
 Combined K-S:       0.2065    0.011      0.009

Note: Ties exist in combined dataset;
      there are 29 unique values out of 244 observations.

. ksmirnov base_pf_RI1 if question=="RI" & control==0, by(incover75) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0094    0.990
 1:                 -0.1679    0.038
 Combined K-S:       0.1679    0.077      0.066

Note: Ties exist in combined dataset;
      there are 29 unique values out of 244 observations.

. 
. ksmirnov end_pf_RI1 if question=="RI" & control==0, by(fem) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0595    0.654
 1:                 -0.0220    0.943
 Combined K-S:       0.0595    0.984      0.972

Note: Ties exist in combined dataset;
      there are 25 unique values out of 244 observations.

. ksmirnov end_pf_RI1 if question=="RI" & control==0, by(age55andup) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0677    0.573
 1:                 -0.0285    0.906
 Combined K-S:       0.0677    0.943      0.918

Note: Ties exist in combined dataset;
      there are 25 unique values out of 244 observations.

. ksmirnov end_pf_RI1 if question=="RI" & control==0, by(highfinlit) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0258    0.929
 1:                 -0.1584    0.064
 Combined K-S:       0.1584    0.127      0.110

Note: Ties exist in combined dataset;
      there are 25 unique values out of 243 observations.

. ksmirnov end_pf_RI1 if question=="RI" & control==0, by(college) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0109    0.986
 1:                 -0.2036    0.006
 Combined K-S:       0.2036    0.013      0.011

Note: Ties exist in combined dataset;
      there are 25 unique values out of 244 observations.

. ksmirnov end_pf_RI1 if question=="RI" & control==0, by(incover75) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0230    0.941
 1:                 -0.1055    0.276
 Combined K-S:       0.1055    0.541      0.499

Note: Ties exist in combined dataset;
      there are 25 unique values out of 244 observations.

. 
. restore

. 
. preserve

. 
. rename base_pf_RI1 pf_RI1_init

. rename end_pf_RI1 pf_RI1_end

. keep prim_key question fem age55andup highfinlit college incover75 all pf_RI1_init pf_RI1_end control

. reshape long pf_RI1, i(prim_key question fem age55andup highfinlit college incover75 all control) j(stage) string
(note: j = _end _init)

Data                               wide   ->   long
-----------------------------------------------------------------------------
Number of obs.                      325   ->     650
Number of variables                  11   ->      11
j variable (2 values)                     ->   stage
xij variables:
                 pf_RI1_end pf_RI1_init   ->   pf_RI1
-----------------------------------------------------------------------------

. gen init = stage=="_init"

. ksmirnov pf_RI1 if question=="RI" & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1434    0.007
 1:                 -0.0574    0.448
 Combined K-S:       0.1434    0.013      0.011

Note: Ties exist in combined dataset;
      there are 35 unique values out of 488 observations.

. ksmirnov pf_RI1 if question=="RI" & control==1, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0741    0.641
 1:                 -0.0494    0.821
 Combined K-S:       0.0741    0.979      0.931

Note: Ties exist in combined dataset;
      there are 17 unique values out of 162 observations.

. ksmirnov pf_RI1 if question=="RI" & fem==0 & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1121    0.260
 1:                 -0.1028    0.323
 Combined K-S:       0.1121    0.511      0.410

Note: Ties exist in combined dataset;
      there are 26 unique values out of 214 observations.

. ksmirnov pf_RI1 if question=="RI" & fem==1 & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1898    0.007
 1:                 -0.0219    0.936
 Combined K-S:       0.1898    0.014      0.014

Note: Ties exist in combined dataset;
      there are 29 unique values out of 274 observations.

. ksmirnov pf_RI1 if question=="RI" & age55andup==0 & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1615    0.034
 1:                 -0.0769    0.463
 Combined K-S:       0.1615    0.067      0.066

Note: Ties exist in combined dataset;
      there are 32 unique values out of 260 observations.

. ksmirnov pf_RI1 if question=="RI" & age55andup==1 & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1404    0.106
 1:                 -0.0351    0.869
 Combined K-S:       0.1404    0.211      0.168

Note: Ties exist in combined dataset;
      there are 22 unique values out of 228 observations.

. ksmirnov pf_RI1 if question=="RI" & highfinlit==0 & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1698    0.010
 1:                 -0.0377    0.797
 Combined K-S:       0.1698    0.020      0.020

Note: Ties exist in combined dataset;
      there are 32 unique values out of 318 observations.

. ksmirnov pf_RI1 if question=="RI" & highfinlit==1 & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1310    0.237
 1:                 -0.0952    0.467
 Combined K-S:       0.1310    0.467      0.470

Note: Ties exist in combined dataset;
      there are 17 unique values out of 168 observations.

. ksmirnov pf_RI1 if question=="RI" & college==0 & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1491    0.079
 1:                 -0.0175    0.966
 Combined K-S:       0.1491    0.158      0.156

Note: Ties exist in combined dataset;
      there are 26 unique values out of 228 observations.

. ksmirnov pf_RI1 if question=="RI" & college==1 & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1462    0.062
 1:                 -0.0923    0.330
 Combined K-S:       0.1462    0.124      0.110

Note: Ties exist in combined dataset;
      there are 26 unique values out of 260 observations.

. ksmirnov pf_RI1 if question=="RI" & incover75==0 & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.1800    0.008
 1:                 -0.0400    0.787
 Combined K-S:       0.1800    0.016      0.015

Note: Ties exist in combined dataset;
      there are 31 unique values out of 300 observations.

. ksmirnov pf_RI1 if question=="RI" & incover75==1 & control==0, by(init) exact

Two-sample Kolmogorov-Smirnov test for equality of distribution functions

 Smaller group       D       P-value      Exact
 ----------------------------------------------
 0:                  0.0957    0.422
 1:                 -0.0851    0.506
 Combined K-S:       0.0957    0.782      0.785

Note: Ties exist in combined dataset;
      there are 23 unique values out of 188 observations.

. 
. restore

. 
.         * * * * * *
.         * Table 4 *
.         * * * * * *
. 
. eststo clear

.         eststo: reg pf_revision1 riXtreated femXtreated femXri maleXgapXri femXgapXri if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     147
                                                       F(  5,   141) =    2.11
                                                       Prob > F      =  0.0673
                                                       R-squared     =  0.0209
                                                       Root MSE      =  5.2559

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |   -1.28106   .6800432    -1.88   0.062    -2.625458    .0633391
 femXtreated |    .783787   1.608141     0.49   0.627    -2.395398    3.962972
      femXri |  -1.205314   .9336122    -1.29   0.199    -3.051001    .6403734
 maleXgapXri |  -.0677651   .0769454    -0.88   0.380    -.2198809    .0843508
  femXgapXri |   .0340166   .1758172     0.19   0.847    -.3135619    .3815951
       _cons |  -.0555556   .1494842    -0.37   0.711    -.3510756    .2399645
------------------------------------------------------------------------------
(est1 stored)

.         test maleXgapXri = femXgapXri 

 ( 1)  maleXgapXri - femXgapXri = 0

       F(  1,   141) =    0.28
            Prob > F =    0.5967

.         eststo: reg pf_revision1 riXtreated femXtreated femXri maleXgapXri femXgapXri if ri==1 & treatment=="spf", robust

Linear regression                                      Number of obs =     178
                                                       F(  5,   172) =    2.99
                                                       Prob > F      =  0.0129
                                                       R-squared     =  0.2348
                                                       Root MSE      =  4.8119

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |  -.2921607   .5939455    -0.49   0.623    -1.464521    .8801999
 femXtreated |   .5682555   1.392505     0.41   0.684    -2.180343    3.316854
      femXri |  -.8279121   1.135474    -0.73   0.467    -3.069171    1.413347
 maleXgapXri |  -.0624599   .1318311    -0.47   0.636     -.322675    .1977552
  femXgapXri |   .4042088   .1434722     2.82   0.005      .121016    .6874016
       _cons |  -.3928571   .2678651    -1.47   0.144    -.9215832    .1358689
------------------------------------------------------------------------------
(est2 stored)

.         test maleXgapXri = femXgapXri 

 ( 1)  maleXgapXri - femXgapXri = 0

       F(  1,   172) =    5.74
            Prob > F =    0.0177

.         eststo: reg pf_revision1 riXtreated incover75Xtreated incover75Xri incunder75XgapXri incover75XgapXri if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     146
                                                       F(  5,   140) =    1.60
                                                       Prob > F      =  0.1632
                                                       R-squared     =  0.0183
                                                       Root MSE      =  5.2806

-----------------------------------------------------------------------------------
                  |               Robust
     pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
------------------+----------------------------------------------------------------
       riXtreated |  -1.315992   1.376198    -0.96   0.341     -4.03681    1.404825
incover75Xtreated |   .8365058   1.484664     0.56   0.574    -2.098754    3.771766
     incover75Xri |   1.041667   .9088925     1.15   0.254    -.7552626    2.838596
incunder75XgapXri |  -.0769285   .1074083    -0.72   0.475    -.2892805    .1354235
 incover75XgapXri |   .1051923   .1428123     0.74   0.463    -.1771552    .3875398
            _cons |  -1.166667   .8821678    -1.32   0.188     -2.91076    .5774265
-----------------------------------------------------------------------------------
(est3 stored)

.         test incunder75XgapXri = incover75XgapXri 

 ( 1)  incunder75XgapXri - incover75XgapXri = 0

       F(  1,   140) =    1.04
            Prob > F =    0.3099

.         eststo: reg pf_revision1 riXtreated incover75Xtreated incover75Xri incunder75XgapXri incover75XgapXri if ri==1 & treatment=="spf", robust

Linear regression                                      Number of obs =     178
                                                       F(  5,   172) =    3.91
                                                       Prob > F      =  0.0022
                                                       R-squared     =  0.2173
                                                       Root MSE      =  4.8666

-----------------------------------------------------------------------------------
                  |               Robust
     pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
------------------+----------------------------------------------------------------
       riXtreated |   .4167866   1.323478     0.31   0.753    -2.195564    3.029137
incover75Xtreated |   .0825414   1.483246     0.06   0.956    -2.845167    3.010249
     incover75Xri |   .5449616   1.307534     0.42   0.677    -2.035918    3.125841
incunder75XgapXri |    .352618   .1439141     2.45   0.015     .0685529    .6366831
 incover75XgapXri |   .5379007   .1559162     3.45   0.001     .2301452    .8456562
            _cons |  -1.162609   1.201545    -0.97   0.335    -3.534281    1.209064
-----------------------------------------------------------------------------------
(est4 stored)

.         test incunder75XgapXri = incover75XgapXri 

 ( 1)  incunder75XgapXri - incover75XgapXri = 0

       F(  1,   172) =    0.76
            Prob > F =    0.3838

.         eststo: reg pf_revision1 riXtreated collegeXtreated collegeXri nocollegeXgapXri collegeXgapXri if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     147
                                                       F(  5,   141) =    2.10
                                                       Prob > F      =  0.0690
                                                       R-squared     =  0.0399
                                                       Root MSE      =  5.2048

----------------------------------------------------------------------------------
                 |               Robust
    pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-----------------+----------------------------------------------------------------
      riXtreated |   -1.74988   1.495013    -1.17   0.244    -4.705419    1.205659
 collegeXtreated |   1.436523   1.565039     0.92   0.360    -1.657453    4.530498
      collegeXri |    1.04902   .9018004     1.16   0.247     -.733778    2.831817
nocollegeXgapXri |  -.1480509   .1297814    -1.14   0.256    -.4046198     .108518
  collegeXgapXri |   .1818548   .1026569     1.77   0.079    -.0210909    .3848006
           _cons |  -1.166667   .8861339    -1.32   0.190    -2.918493    .5851593
----------------------------------------------------------------------------------
(est5 stored)

.         test nocollegeXgapXri = collegeXgapXri

 ( 1)  nocollegeXgapXri - collegeXgapXri = 0

       F(  1,   141) =    3.97
            Prob > F =    0.0481

.         eststo: reg pf_revision1 riXtreated collegeXtreated collegeXri nocollegeXgapXri collegeXgapXri if ri==1 & treatment=="spf", robust

Linear regression                                      Number of obs =     178
                                                       F(  5,   172) =    9.35
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.2249
                                                       Root MSE      =  4.8429

----------------------------------------------------------------------------------
                 |               Robust
    pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-----------------+----------------------------------------------------------------
      riXtreated |   .4401401   1.413838     0.31   0.756    -2.350567    3.230847
 collegeXtreated |   -.392107   1.497129    -0.26   0.794    -3.347219    2.563005
      collegeXri |   1.261074   1.262039     1.00   0.319    -1.230003    3.752152
nocollegeXgapXri |   .3172472   .1598311     1.98   0.049     .0017642    .6327302
  collegeXgapXri |   .5261306   .0895067     5.88   0.000     .3494576    .7028037
           _cons |  -1.466957    1.23692    -1.19   0.237    -3.908453    .9745398
----------------------------------------------------------------------------------
(est6 stored)

.         test nocollegeXgapXri = collegeXgapXri

 ( 1)  nocollegeXgapXri - collegeXgapXri = 0

       F(  1,   172) =    1.30
            Prob > F =    0.2558

.         eststo: reg pf_revision1 riXtreated highfinlitXtreated highfinlitXri lofinlitXgapXri highfinlitXgapXri if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     146
                                                       F(  5,   140) =    1.08
                                                       Prob > F      =  0.3754
                                                       R-squared     =  0.0234
                                                       Root MSE      =  5.2669

------------------------------------------------------------------------------------
                   |               Robust
      pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------------+----------------------------------------------------------------
        riXtreated |  -2.221766   1.134146    -1.96   0.052    -4.464034    .0205028
highfinlitXtreated |   2.862504   1.708945     1.68   0.096    -.5161725    6.241181
     highfinlitXri |   -.810049   1.251889    -0.65   0.519    -3.285101    1.665003
   lofinlitXgapXri |  -.0715276   .1345593    -0.53   0.596    -.3375586    .1945034
 highfinlitXgapXri |   .0004586   .0195089     0.02   0.981    -.0381115    .0390286
             _cons |  -.3958333   .2062238    -1.92   0.057    -.8035489    .0118822
------------------------------------------------------------------------------------
(est7 stored)

.         test lofinlitXgapXri = highfinlitXgapXri

 ( 1)  lofinlitXgapXri - highfinlitXgapXri = 0

       F(  1,   140) =    0.28
            Prob > F =    0.5973

.         eststo: reg pf_revision1 riXtreated highfinlitXtreated highfinlitXri lofinlitXgapXri highfinlitXgapXri if ri==1 & treatment=="spf", robust

Linear regression                                      Number of obs =     178
                                                       F(  5,   172) =    4.05
                                                       Prob > F      =  0.0017
                                                       R-squared     =  0.2179
                                                       Root MSE      =  4.8647

------------------------------------------------------------------------------------
                   |               Robust
      pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------------+----------------------------------------------------------------
        riXtreated |     .11658   1.125936     0.10   0.918    -2.105851    2.339011
highfinlitXtreated |   .5169371   1.255566     0.41   0.681    -1.961364    2.995239
     highfinlitXri |   .7746667    1.01664     0.76   0.447    -1.232031    2.781364
   lofinlitXgapXri |     .37263   .1339131     2.78   0.006     .1083053    .6369547
 highfinlitXgapXri |   .4640141   .5149651     0.90   0.369    -.5524508    1.480479
             _cons |  -1.124667   .9588854    -1.17   0.242    -3.017365    .7680314
------------------------------------------------------------------------------------
(est8 stored)

.         test lofinlitXgapXri = highfinlitXgapXri

 ( 1)  lofinlitXgapXri - highfinlitXgapXri = 0

       F(  1,   172) =    0.03
            Prob > F =    0.8638

.         eststo: reg pf_revision1 riXtreated age55andupXtreated age55andupXri ageunder55XgapXri age55andupXgapXri if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     147
                                                       F(  5,   141) =    1.17
                                                       Prob > F      =  0.3286
                                                       R-squared     =  0.0130
                                                       Root MSE      =  5.2772

------------------------------------------------------------------------------------
                   |               Robust
      pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------------+----------------------------------------------------------------
        riXtreated |  -.7276006   1.349964    -0.54   0.591    -3.396387    1.941186
age55andupXtreated |  -.3541134   1.723762    -0.21   0.838    -3.761872    3.053645
     age55andupXri |   1.158333   1.063566     1.09   0.278    -.9442632     3.26093
 ageunder55XgapXri |  -.0195635    .060209    -0.32   0.746    -.1385925    .0994656
 age55andupXgapXri |  -.0202254   .1886773    -0.11   0.915    -.3932274    .3527766
             _cons |     -1.325   1.004504    -1.32   0.189    -3.310836    .6608359
------------------------------------------------------------------------------------
(est9 stored)

.         test ageunder55XgapXri = age55andupXgapXri 

 ( 1)  ageunder55XgapXri - age55andupXgapXri = 0

       F(  1,   141) =    0.00
            Prob > F =    0.9973

.         eststo: reg pf_revision1 riXtreated age55andupXtreated age55andupXri ageunder55XgapXri age55andupXgapXri if ri==1 & treatment=="spf", robust

Linear regression                                      Number of obs =     178
                                                       F(  5,   172) =    5.99
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.2208
                                                       Root MSE      =  4.8557

------------------------------------------------------------------------------------
                   |               Robust
      pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------------+----------------------------------------------------------------
        riXtreated |   .0249742   1.071491     0.02   0.981    -2.089991    2.139939
age55andupXtreated |   .9736081   1.647526     0.59   0.555    -2.278365    4.225581
     age55andupXri |  -.0688491   1.490114    -0.05   0.963    -3.010114    2.872416
 ageunder55XgapXri |   .3278333   .1610595     2.04   0.043     .0099257     .645741
 age55andupXgapXri |   .5255853   .1055381     4.98   0.000     .3172687     .733902
             _cons |  -.9017391   .9170315    -0.98   0.327    -2.711824    .9083455
------------------------------------------------------------------------------------
(est10 stored)

.         test ageunder55XgapXri = age55andupXgapXri 

 ( 1)  ageunder55XgapXri - age55andupXgapXri = 0

       F(  1,   172) =    1.05
            Prob > F =    0.3059

.         eststo: reg pf_revision1 riXtreated affectedXtreated affectedXri noaffectedXgapXri affectedXgapXri if ri==1 & treatment=="food", robust
note: affectedXri omitted because of collinearity

Linear regression                                      Number of obs =     147
                                                       F(  4,   142) =    0.74
                                                       Prob > F      =  0.5693
                                                       R-squared     =  0.0232
                                                       Root MSE      =  5.2313

-----------------------------------------------------------------------------------
                  |               Robust
     pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
------------------+----------------------------------------------------------------
       riXtreated |  -.4959294   .7831023    -0.63   0.528    -2.043975    1.052116
 affectedXtreated |  -.9371254   1.345056    -0.70   0.487    -3.596046    1.721796
      affectedXri |          0  (omitted)
noaffectedXgapXri |  -.0634108   .0791308    -0.80   0.424    -.2198375    .0930158
  affectedXgapXri |   .0203391   .1775786     0.11   0.909    -.3307003    .3713785
            _cons |  -.7317073   .5279142    -1.39   0.168    -1.775294    .3118793
-----------------------------------------------------------------------------------
(est11 stored)

.         test noaffectedXgapXri = affectedXgapXri 

 ( 1)  noaffectedXgapXri - affectedXgapXri = 0

       F(  1,   142) =    0.19
            Prob > F =    0.6673

.         eststo: reg pf_revision1 riXtreated affectedXtreated affectedXri noaffectedXgapXri affectedXgapXri if ri==1 & treatment=="spf", robust
note: affectedXri omitted because of collinearity

Linear regression                                      Number of obs =     178
                                                       F(  4,   173) =    5.29
                                                       Prob > F      =  0.0005
                                                       R-squared     =  0.2501
                                                       Root MSE      =  4.7497

-----------------------------------------------------------------------------------
                  |               Robust
     pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
------------------+----------------------------------------------------------------
       riXtreated |   1.257922     .79877     1.57   0.117    -.3186674    2.834511
 affectedXtreated |  -1.706906   .7506993    -2.27   0.024    -3.188615   -.2251971
      affectedXri |          0  (omitted)
noaffectedXgapXri |     .24767   .1374049     1.80   0.073    -.0235358    .5188758
  affectedXgapXri |   .4235383   .1819516     2.33   0.021     .0644074    .7826692
            _cons |      -.931   .7240107    -1.29   0.200    -2.360032    .4980315
-----------------------------------------------------------------------------------
(est12 stored)

.         test noaffectedXgapXri = affectedXgapXri 

 ( 1)  noaffectedXgapXri - affectedXgapXri = 0

       F(  1,   173) =    0.59
            Prob > F =    0.4416

. esttab using Output\Table4.csv, se replace star(* 0.10 ** 0.05 *** 0.01) r2
(output written to Output\Table4.csv)

. 
.         * * * * * *
.         * Table 5 *
.         * * * * * *
. 
. preserve

. 
. gen delta_var = varb - vara
(11 missing values generated)

. egen med_delta_var = median(delta_var)

. gen hi_delta_var = delta_var>med_delta_var if delta_var!=.
(11 missing values generated)

. drop med_delta_var

. rename dens_revision densmean_rev

. gen densmedian_rev = medianb - mediana
(11 missing values generated)

. 
. eststo clear

. eststo: reg pf_revision1 densmean_rev, robust

Linear regression                                      Number of obs =     314
                                                       F(  1,   312) =   15.38
                                                       Prob > F      =  0.0001
                                                       R-squared     =  0.1685
                                                       Root MSE      =  4.4275

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
densmean_rev |   .8110794   .2068337     3.92   0.000     .4041142    1.218045
       _cons |   -.891134   .2381932    -3.74   0.000    -1.359802    -.422466
------------------------------------------------------------------------------
(est1 stored)

. test densmean_rev==1

 ( 1)  densmean_rev = 1

       F(  1,   312) =    0.83
            Prob > F =    0.3617

. eststo: reg pf_revision1 densmean_rev delta_var, robust

Linear regression                                      Number of obs =     314
                                                       F(  2,   311) =   12.21
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.2128
                                                       Root MSE      =   4.315

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
densmean_rev |   .7254396   .2181694     3.33   0.001     .2961648    1.154714
   delta_var |   .0426716   .0143865     2.97   0.003     .0143645    .0709787
       _cons |   -.662413   .2087557    -3.17   0.002    -1.073165    -.251661
------------------------------------------------------------------------------
(est2 stored)

. test densmean_rev==1

 ( 1)  densmean_rev = 1

       F(  1,   311) =    1.58
            Prob > F =    0.2092

. test (densmean_rev==1) (delta_var==0)

 ( 1)  densmean_rev = 1
 ( 2)  delta_var = 0

       F(  2,   311) =    4.65
            Prob > F =    0.0102

. eststo: reg pf_revision1 densmedian_rev, robust

Linear regression                                      Number of obs =     314
                                                       F(  1,   312) =   16.99
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.1854
                                                       Root MSE      =  4.3824

--------------------------------------------------------------------------------
               |               Robust
  pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
densmedian_rev |   .8432001   .2045622     4.12   0.000     .4407042    1.245696
         _cons |  -.8811054   .2339538    -3.77   0.000    -1.341432   -.4207786
--------------------------------------------------------------------------------
(est3 stored)

. test densmedian_rev==1

 ( 1)  densmedian_rev = 1

       F(  1,   312) =    0.59
            Prob > F =    0.4439

. eststo: reg pf_revision1 densmedian_rev delta_var, robust

Linear regression                                      Number of obs =     314
                                                       F(  2,   311) =   12.35
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.2436
                                                       Root MSE      =  4.2299

--------------------------------------------------------------------------------
               |               Robust
  pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
densmedian_rev |   .7886877   .2105241     3.75   0.000      .374456    1.202919
     delta_var |   .0482142   .0136076     3.54   0.000     .0214395    .0749888
         _cons |  -.5961663   .2094364    -2.85   0.005    -1.008258   -.1840749
--------------------------------------------------------------------------------
(est4 stored)

. test densmedian_rev==1

 ( 1)  densmedian_rev = 1

       F(  1,   311) =    1.01
            Prob > F =    0.3163

. test (densmedian_rev==1) (delta_var==0)

 ( 1)  densmedian_rev = 1
 ( 2)  delta_var = 0

       F(  2,   311) =    7.09
            Prob > F =    0.0010

. esttab using Output\Table5.csv, se replace star(* 0.10 ** 0.05 *** 0.01) r2
(output written to Output\Table5.csv)

. 
. restore

. 
.         * * * * * * * * * * *
.         * Appendix Table A1 *
.         * * * * * * * * * * *
. 
. eststo clear

. eststo: reg log_info_belief food fem femXfood if ri==1, robust

Linear regression                                      Number of obs =     325
                                                       F(  3,   321) =   33.42
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.1957
                                                       Root MSE      =  .87906

------------------------------------------------------------------------------
             |               Robust
log_info_b~f |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        food |  -.9604122   .1475224    -6.51   0.000    -1.250645   -.6701793
         fem |  -.4392213   .1312727    -3.35   0.001    -.6974847   -.1809578
    femXfood |   .2797376   .1957839     1.43   0.154    -.1054439    .6649192
       _cons |   -.375331   .0860986    -4.36   0.000    -.5447198   -.2059421
------------------------------------------------------------------------------
(est1 stored)

. test fem femXfood

 ( 1)  fem = 0
 ( 2)  femXfood = 0

       F(  2,   321) =    6.20
            Prob > F =    0.0023

. eststo: reg log_info_belief food age55andup age55andupXfood if ri==1, robust

Linear regression                                      Number of obs =     325
                                                       F(  3,   321) =   23.72
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.1721
                                                       Root MSE      =  .89185

---------------------------------------------------------------------------------
                |               Robust
log_info_belief |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
----------------+----------------------------------------------------------------
           food |  -.6919945   .1390572    -4.98   0.000    -.9655732   -.4184158
     age55andup |  -.0325827   .1415755    -0.23   0.818    -.3111157    .2459502
age55andupXfood |   -.199491   .1985771    -1.00   0.316    -.5901678    .1911859
          _cons |  -.6153923   .0911133    -6.75   0.000    -.7946469   -.4361377
---------------------------------------------------------------------------------
(est2 stored)

. test age55andup age55andupXfood

 ( 1)  age55andup = 0
 ( 2)  age55andupXfood = 0

       F(  2,   321) =    1.42
            Prob > F =    0.2444

. eststo: reg log_info_belief food highfinlit highfinlitXfood if ri==1, robust

Linear regression                                      Number of obs =     324
                                                       F(  3,   320) =   31.42
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.1989
                                                       Root MSE      =  .87856

---------------------------------------------------------------------------------
                |               Robust
log_info_belief |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
----------------+----------------------------------------------------------------
           food |  -.8162858   .1223134    -6.67   0.000    -1.056926   -.5756459
     highfinlit |   .3762851   .1345568     2.80   0.005     .1115575    .6410127
highfinlitXfood |   .0008256   .1985006     0.00   0.997    -.3897054    .3913566
          _cons |  -.7499829   .0888719    -8.44   0.000    -.9248298   -.5751359
---------------------------------------------------------------------------------
(est3 stored)

. test highfinlit highfinlitXfood

 ( 1)  highfinlit = 0
 ( 2)  highfinlitXfood = 0

       F(  2,   320) =    7.25
            Prob > F =    0.0008

. eststo: reg log_info_belief food college collegeXfood if ri==1, robust

Linear regression                                      Number of obs =     325
                                                       F(  3,   321) =   32.90
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.2055
                                                       Root MSE      =  .87369

------------------------------------------------------------------------------
             |               Robust
log_info_b~f |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        food |  -.7963788   .1456066    -5.47   0.000    -1.082843   -.5099151
     college |   .3709522   .1391324     2.67   0.008     .0972256    .6446788
collegeXfood |   .0405853   .1949229     0.21   0.835    -.3429025    .4240731
       _cons |  -.8253832    .113365    -7.28   0.000    -1.048415    -.602351
------------------------------------------------------------------------------
(est4 stored)

. test college collegeXfood

 ( 1)  college = 0
 ( 2)  collegeXfood = 0

       F(  2,   321) =    8.10
            Prob > F =    0.0004

. eststo: reg log_info_belief food incover75 incover75Xfood if ri==1, robust

Linear regression                                      Number of obs =     324
                                                       F(  3,   320) =   28.65
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.1920
                                                       Root MSE      =  .88233

--------------------------------------------------------------------------------
               |               Robust
log_info_bel~f |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
          food |  -.8757211   .1299208    -6.74   0.000    -1.131328   -.6201143
     incover75 |   .2423901   .1347264     1.80   0.073    -.0226714    .5074515
incover75Xfood |   .1604114   .1939054     0.83   0.409     -.221079    .5419018
         _cons |   -.719362   .0943047    -7.63   0.000    -.9048976   -.5338265
--------------------------------------------------------------------------------
(est5 stored)

. test incover75 incover75Xfood

 ( 1)  incover75 = 0
 ( 2)  incover75Xfood = 0

       F(  2,   320) =    5.79
            Prob > F =    0.0034

. eststo: reg log_info_belief food fem femXfood age55andup age55andupXfood highfinlit highfinlitXfood college collegeXfood incover75 incover75Xfood if ri==1, robust

Linear regression                                      Number of obs =     323
                                                       F( 11,   311) =   15.38
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.2668
                                                       Root MSE      =  .85243

---------------------------------------------------------------------------------
                |               Robust
log_info_belief |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
----------------+----------------------------------------------------------------
           food |   -.927187   .2359869    -3.93   0.000     -1.39152   -.4628543
            fem |  -.3452742   .1323778    -2.61   0.010    -.6057436   -.0848048
       femXfood |    .258211   .1932749     1.34   0.183    -.1220807    .6385028
     age55andup |   .0549579   .1478811     0.37   0.710     -.236016    .3459319
age55andupXfood |  -.2958403   .2012199    -1.47   0.143    -.6917649    .1000843
     highfinlit |    .244781   .1554229     1.57   0.116    -.0610324    .5505944
highfinlitXfood |   .1062866   .2100319     0.51   0.613    -.3069766    .5195498
        college |   .2648503   .1385812     1.91   0.057     -.007825    .5375257
   collegeXfood |   .0699964   .2052867     0.34   0.733    -.3339302    .4739229
      incover75 |   .1298282   .1165182     1.11   0.266    -.0994354    .3590918
 incover75Xfood |   .1094919   .1903335     0.58   0.566    -.2650124    .4839961
          _cons |  -.7198558   .1666149    -4.32   0.000    -1.047691   -.3920207
---------------------------------------------------------------------------------
(est6 stored)

. test fem femXfood age55andup age55andupXfood highfinlit highfinlitXfood college collegeXfood incover75 incover75Xfood 

 ( 1)  fem = 0
 ( 2)  femXfood = 0
 ( 3)  age55andup = 0
 ( 4)  age55andupXfood = 0
 ( 5)  highfinlit = 0
 ( 6)  highfinlitXfood = 0
 ( 7)  college = 0
 ( 8)  collegeXfood = 0
 ( 9)  incover75 = 0
 (10)  incover75Xfood = 0

       F( 10,   311) =    4.97
            Prob > F =    0.0000

. esttab using Output\TableA1.csv, se replace star(* 0.10 ** 0.05 *** 0.01) r2
(output written to Output\TableA1.csv)

. 
.         * * * * * * * * * * *
.         * Appendix Table A2 *
.         * * * * * * * * * * *
. 
. eststo clear

. 
. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated uncXfem femXgapXri maleXgapXri uncXfemXgapXri uncXmaleXgapXri if ri==1 & treatment=="food", r
> obust

Linear regression                                      Number of obs =     144
                                                       F( 10,   133) =    1.73
                                                       Prob > F      =  0.0800
                                                       R-squared     =  0.1015
                                                       Root MSE      =  5.1801

---------------------------------------------------------------------------------
                |               Robust
   pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
----------------+----------------------------------------------------------------
     riXtreated |  -.3649025   .8328502    -0.44   0.662    -2.012248    1.282443
            unc |  -.9226147   1.243378    -0.74   0.459    -3.381968    1.536739
    uncXtreated |  -2.464535   2.040585    -1.21   0.229    -6.500734    1.571664
         femXri |  -.8986412   .9196808    -0.98   0.330    -2.717734    .9204519
    femXtreated |   .9785236   1.562691     0.63   0.532    -2.112418    4.069466
        uncXfem |  -1.065029   2.032262    -0.52   0.601    -5.084764    2.954705
     femXgapXri |   .0011399   .0355612     0.03   0.974    -.0691987    .0714786
    maleXgapXri |  -.0397045   .0367543    -1.08   0.282     -.112403     .032994
 uncXfemXgapXri |  -.0799317   .2758485    -0.29   0.772    -.6255493     .465686
uncXmaleXgapXri |  -.1369871   .1248978    -1.10   0.275    -.3840301    .1100558
          _cons |   .3032391   .4845015     0.63   0.532    -.6550862    1.261564
---------------------------------------------------------------------------------
(est1 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated uncXfem femXgapXri maleXgapXri uncXfemXgapXri uncXmaleXgapXri if ri==1 & treatment=="spf", ro
> bust

Linear regression                                      Number of obs =     175
                                                       F( 10,   164) =   21.58
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.3629
                                                       Root MSE      =  4.4916

---------------------------------------------------------------------------------
                |               Robust
   pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
----------------+----------------------------------------------------------------
     riXtreated |  -.6534954   .7209332    -0.91   0.366    -2.077003     .770012
            unc |  -2.013741   1.161071    -1.73   0.085    -4.306316     .278835
    uncXtreated |  -.1248529    1.35093    -0.09   0.926    -2.792312    2.542606
         femXri |  -.6812746    1.08332    -0.63   0.530    -2.820327    1.457778
    femXtreated |   .8733782   1.343272     0.65   0.516    -1.778959    3.525716
        uncXfem |   .3475626   1.511496     0.23   0.818    -2.636938    3.332063
     femXgapXri |   .1113688   .0096194    11.58   0.000      .092375    .1303625
    maleXgapXri |   .0098004   .0293914     0.33   0.739     -.048234    .0678348
 uncXfemXgapXri |   .4492167   .1503937     2.99   0.003     .1522592    .7461743
uncXmaleXgapXri |  -.3581053   .3567466    -1.00   0.317    -1.062514    .3463032
          _cons |   .6140131   .5907762     1.04   0.300     -.552495    1.780521
---------------------------------------------------------------------------------
(est2 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated uncXfem uncXgapXri femXgapXri maleXgapXri uncXfemXgapXri if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     144
                                                       F( 10,   133) =    1.73
                                                       Prob > F      =  0.0800
                                                       R-squared     =  0.1015
                                                       Root MSE      =  5.1801

--------------------------------------------------------------------------------
               |               Robust
  pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
    riXtreated |  -.3649025   .8328502    -0.44   0.662    -2.012248    1.282443
           unc |  -.9226147   1.243378    -0.74   0.459    -3.381968    1.536739
   uncXtreated |  -2.464535   2.040585    -1.21   0.229    -6.500734    1.571664
        femXri |  -.8986412   .9196808    -0.98   0.330    -2.717734    .9204519
   femXtreated |   .9785236   1.562691     0.63   0.532    -2.112418    4.069466
       uncXfem |  -1.065029   2.032262    -0.52   0.601    -5.084764    2.954705
    uncXgapXri |  -.1369871   .1248978    -1.10   0.275    -.3840301    .1100558
    femXgapXri |   .0011399   .0355612     0.03   0.974    -.0691987    .0714786
   maleXgapXri |  -.0397045   .0367543    -1.08   0.282     -.112403     .032994
uncXfemXgapXri |   .0570555     .27875     0.20   0.838    -.4943011    .6084121
         _cons |   .3032391   .4845015     0.63   0.532    -.6550862    1.261564
--------------------------------------------------------------------------------
(est3 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated uncXfem uncXgapXri femXgapXri maleXgapXri uncXfemXgapXri if ri==1 & treatment=="spf", robust

Linear regression                                      Number of obs =     175
                                                       F( 10,   164) =   21.58
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.3629
                                                       Root MSE      =  4.4916

--------------------------------------------------------------------------------
               |               Robust
  pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
    riXtreated |  -.6534954   .7209332    -0.91   0.366    -2.077003     .770012
           unc |  -2.013741   1.161071    -1.73   0.085    -4.306316     .278835
   uncXtreated |  -.1248529    1.35093    -0.09   0.926    -2.792312    2.542606
        femXri |  -.6812746    1.08332    -0.63   0.530    -2.820327    1.457778
   femXtreated |   .8733782   1.343272     0.65   0.516    -1.778959    3.525716
       uncXfem |   .3475626   1.511496     0.23   0.818    -2.636938    3.332063
    uncXgapXri |  -.3581053   .3567466    -1.00   0.317    -1.062514    .3463032
    femXgapXri |   .1113688   .0096194    11.58   0.000      .092375    .1303625
   maleXgapXri |   .0098004   .0293914     0.33   0.739     -.048234    .0678348
uncXfemXgapXri |    .807322   .3765356     2.14   0.033     .0638394    1.550805
         _cons |   .6140131   .5907762     1.04   0.300     -.552495    1.780521
--------------------------------------------------------------------------------
(est4 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated femXgapXri maleXgapXri uncXfemXgapXri uncXmaleXgapXri if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     144
                                                       F(  9,   134) =    1.80
                                                       Prob > F      =  0.0737
                                                       R-squared     =  0.0999
                                                       Root MSE      =  5.1652

---------------------------------------------------------------------------------
                |               Robust
   pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
----------------+----------------------------------------------------------------
     riXtreated |  -.3768232    .832837    -0.45   0.652     -2.02403    1.270383
            unc |  -1.501348    1.38318    -1.09   0.280    -4.237037    1.234341
    uncXtreated |  -2.428516   2.047722    -1.19   0.238    -6.478553    1.621522
         femXri |  -1.300539   .9957787    -1.31   0.194    -3.270016    .6689377
    femXtreated |   .9436617   1.587078     0.59   0.553    -2.195303    4.082626
     femXgapXri |  -.0182839   .0570293    -0.32   0.749    -.1310779    .0945101
    maleXgapXri |  -.0236757   .0506534    -0.47   0.641    -.1238594    .0765079
 uncXfemXgapXri |  -.0375827   .2397057    -0.16   0.876    -.5116788    .4365133
uncXmaleXgapXri |  -.1715834   .1323191    -1.30   0.197    -.4332875    .0901207
          _cons |   .5283019   .5546695     0.95   0.343    -.5687377    1.625342
---------------------------------------------------------------------------------
(est5 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated femXgapXri maleXgapXri uncXfemXgapXri uncXmaleXgapXri if ri==1 & treatment=="spf", robust

Linear regression                                      Number of obs =     175
                                                       F(  9,   165) =   20.17
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.3627
                                                       Root MSE      =  4.4786

---------------------------------------------------------------------------------
                |               Robust
   pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
----------------+----------------------------------------------------------------
     riXtreated |  -.6090182   .6361825    -0.96   0.340    -1.865126    .6470895
            unc |  -1.800766   .9936908    -1.81   0.072    -3.762754    .1612226
    uncXtreated |  -.1304929   1.355116    -0.10   0.923    -2.806096    2.545111
         femXri |  -.4816109   .9492257    -0.51   0.613    -2.355806    1.392584
    femXtreated |   .8011252   1.204243     0.67   0.507    -1.576587    3.178837
     femXgapXri |   .1138534   .0165612     6.87   0.000     .0811542    .1465525
    maleXgapXri |   .0015814   .0520016     0.03   0.976    -.1010929    .1042557
 uncXfemXgapXri |   .4426436   .1479088     2.99   0.003     .1506057    .7346815
uncXmaleXgapXri |  -.3207898   .2726773    -1.18   0.241    -.8591763    .2175967
          _cons |   .5075258   .5169606     0.98   0.328    -.5131847    1.528236
---------------------------------------------------------------------------------
(est6 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated uncXgapXri femXgapXri maleXgapXri uncXfemXgapXri if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     144
                                                       F(  9,   134) =    1.80
                                                       Prob > F      =  0.0737
                                                       R-squared     =  0.0999
                                                       Root MSE      =  5.1652

--------------------------------------------------------------------------------
               |               Robust
  pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
    riXtreated |  -.3768232    .832837    -0.45   0.652     -2.02403    1.270383
           unc |  -1.501348    1.38318    -1.09   0.280    -4.237037    1.234341
   uncXtreated |  -2.428516   2.047722    -1.19   0.238    -6.478553    1.621522
        femXri |  -1.300539   .9957787    -1.31   0.194    -3.270016    .6689377
   femXtreated |   .9436617   1.587078     0.59   0.553    -2.195303    4.082626
    uncXgapXri |  -.1715834   .1323191    -1.30   0.197    -.4332875    .0901207
    femXgapXri |  -.0182839   .0570293    -0.32   0.749    -.1310779    .0945101
   maleXgapXri |  -.0236757   .0506534    -0.47   0.641    -.1238594    .0765079
uncXfemXgapXri |   .1340006    .220218     0.61   0.544    -.3015523    .5695535
         _cons |   .5283019   .5546695     0.95   0.343    -.5687377    1.625342
--------------------------------------------------------------------------------
(est7 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated uncXgapXri femXgapXri maleXgapXri uncXfemXgapXri if ri==1 & treatment=="spf", robust     

Linear regression                                      Number of obs =     175
                                                       F(  9,   165) =   20.17
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.3627
                                                       Root MSE      =  4.4786

--------------------------------------------------------------------------------
               |               Robust
  pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
    riXtreated |  -.6090182   .6361825    -0.96   0.340    -1.865126    .6470895
           unc |  -1.800766   .9936908    -1.81   0.072    -3.762754    .1612226
   uncXtreated |  -.1304929   1.355116    -0.10   0.923    -2.806096    2.545111
        femXri |  -.4816109   .9492257    -0.51   0.613    -2.355806    1.392584
   femXtreated |   .8011252   1.204243     0.67   0.507    -1.576587    3.178837
    uncXgapXri |  -.3207898   .2726773    -1.18   0.241    -.8591763    .2175967
    femXgapXri |   .1138534   .0165612     6.87   0.000     .0811542    .1465525
   maleXgapXri |   .0015814   .0520016     0.03   0.976    -.1010929    .1042557
uncXfemXgapXri |   .7634333   .2753483     2.77   0.006     .2197731    1.307094
         _cons |   .5075258   .5169606     0.98   0.328    -.5131847    1.528236
--------------------------------------------------------------------------------
(est8 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated uncXfem uncXgapXri femXgapXri maleXgapXri if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     144
                                                       F(  9,   134) =    1.86
                                                       Prob > F      =  0.0637
                                                       R-squared     =  0.1012
                                                       Root MSE      =  5.1616

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |  -.3657097   .8302849    -0.44   0.660    -2.007869    1.276449
         unc |  -.7836761   1.256915    -0.62   0.534    -3.269636    1.702283
 uncXtreated |  -2.500055   2.006999    -1.25   0.215     -6.46955     1.46944
      femXri |   -.802156   .9211756    -0.87   0.385    -2.624081    1.019769
 femXtreated |   1.009599   1.599272     0.63   0.529    -2.153483    4.172681
     uncXfem |  -1.320715   1.598324    -0.83   0.410    -4.481921    1.840491
  uncXgapXri |   -.108595   .1619768    -0.67   0.504     -.428957     .211767
  femXgapXri |   .0183242   .0968725     0.19   0.850    -.1732728    .2099211
 maleXgapXri |  -.0603039   .1111462    -0.54   0.588    -.2801318     .159524
       _cons |   .2492074   .4914043     0.51   0.613    -.7227047    1.221119
------------------------------------------------------------------------------
(est9 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated uncXfem uncXgapXri femXgapXri maleXgapXri if ri==1 & treatment=="spf", robust

Linear regression                                      Number of obs =     175
                                                       F(  9,   165) =    5.54
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.3508
                                                       Root MSE      =  4.5203

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |  -.5053097   .6965177    -0.73   0.469    -1.880546    .8699267
         unc |  -1.228707   1.053169    -1.17   0.245    -3.308133     .850718
 uncXtreated |   .3057887    1.26774     0.24   0.810    -2.197295    2.808873
      femXri |    .054694   .9851824     0.06   0.956    -1.890495    1.999883
 femXtreated |   .3887092    1.26854     0.31   0.760    -2.115953    2.893372
     uncXfem |  -.9335679   1.195266    -0.78   0.436    -3.293555    1.426419
  uncXgapXri |   .4003508   .1464092     2.73   0.007     .1112737    .6894279
  femXgapXri |   .1419646    .039103     3.63   0.000     .0647579    .2191714
 maleXgapXri |  -.2904994   .2047927    -1.42   0.158    -.6948516    .1138527
       _cons |   .2214965   .5290845     0.42   0.676     -.823152    1.266145
------------------------------------------------------------------------------
(est10 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated uncXgapXri femXgapXri maleXgapXri if ri==1 & treatment=="food", robust

Linear regression                                      Number of obs =     144
                                                       F(  8,   135) =    1.77
                                                       Prob > F      =  0.0887
                                                       R-squared     =  0.0977
                                                       Root MSE      =  5.1525

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |  -.3895662   .8356408    -0.47   0.642    -2.042207    1.263074
         unc |  -1.501348   1.378048    -1.09   0.278    -4.226702    1.224006
 uncXtreated |  -2.521843   2.005373    -1.26   0.211    -6.487853    1.444168
      femXri |  -1.300539   .9920838    -1.31   0.192    -3.262576    .6614974
 femXtreated |   1.022522   1.600671     0.64   0.524    -2.143112    4.188156
  uncXgapXri |  -.1018219   .1633543    -0.62   0.534    -.4248865    .2212427
  femXgapXri |   .0252024   .0980296     0.26   0.798    -.1686699    .2190747
 maleXgapXri |  -.0818541   .1153344    -0.71   0.479      -.30995    .1462418
       _cons |   .5283019   .5526113     0.96   0.341    -.5645933    1.621197
------------------------------------------------------------------------------
(est11 stored)

. eststo: reg pf_revision1 riXtreated unc uncXtreated femXri femXtreated uncXgapXri femXgapXri maleXgapXri if ri==1 & treatment=="spf", robust    

Linear regression                                      Number of obs =     175
                                                       F(  8,   166) =    5.88
                                                       Prob > F      =  0.0000
                                                       R-squared     =  0.3492
                                                       Root MSE      =  4.5122

------------------------------------------------------------------------------
             |               Robust
pf_revision1 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  riXtreated |  -.6276553   .6469601    -0.97   0.333    -1.904986    .6496755
         unc |  -1.800766   .9906933    -1.82   0.071    -3.756749     .155217
 uncXtreated |   .4033762   1.300572     0.31   0.757    -2.164418     2.97117
      femXri |  -.4816109   .9463623    -0.51   0.611    -2.350069    1.386847
 femXtreated |   .5429115   1.197979     0.45   0.651    -1.822328     2.90815
  uncXgapXri |   .4135074   .1482498     2.79   0.006     .1208093    .7062055
  femXgapXri |   .1392058    .041359     3.37   0.001     .0575483    .2208633
 maleXgapXri |  -.3177297   .2040986    -1.56   0.121    -.7206933    .0852339
       _cons |   .5075258   .5154011     0.98   0.326    -.5100603    1.525112
------------------------------------------------------------------------------
(est12 stored)

. 
. esttab using Output\TableA2.csv, se replace star(* 0.10 ** 0.05 *** 0.01) r2
(output written to Output\TableA2.csv)

.         
. *****************************************************
. ** PART FOUR: RESULTS FOR PAPER TEXT AND FOOTNOTES **
. *****************************************************
. 
. ** Free-response answers:
. outsheet e4_prices e4_inflation using Output\openanswertexts.csv, comma replace

. 
. ** Perception gap distributions
. preserve

. replace gap = 1.39 - belief_food if randomscreen2==4 & food==1
(41 real changes made)

. replace gap = 1.96 - belief_spf if randomscreen2==4 & spf==1
(40 real changes made)

. if "`version'"=="interval" {
.         replace gap = 0 if spf & (belief_spf >= 1.19 & belief_spf<=3.03)
.         replace gap = 3.03 - belief_spf if spf & belief_spf>3.03 & belief_spf!=.
.         replace gap = 1.19 - belief_spf if spf & belief_spf<1.19
.         }

. sum gap if ri==1 & spf==1 & control==0, d

                             gap
-------------------------------------------------------------
      Percentiles      Smallest
 1%       -33.04         -48.04
 5%       -18.04         -33.04
10%        -8.04         -33.04       Obs                 138
25%        -4.04         -32.04       Sum of Wgt.         138

50%        -1.04                      Mean          -3.612464
                        Largest       Std. Dev.      7.063722
75%         -.04            .96
90%         -.04            .96       Variance       49.89617
95%          .96           1.96       Skewness      -3.678043
99%         1.96           1.96       Kurtosis       18.74991

. sum gap if ri==1 & food==1 & control==0, d

                             gap
-------------------------------------------------------------
      Percentiles      Smallest
 1%       -28.61         -38.61
 5%       -18.61         -28.61
10%       -18.61         -23.61       Obs                 106
25%        -8.61         -23.61       Sum of Wgt.         106

50%        -3.61                      Mean          -7.147736
                        Largest       Std. Dev.      7.119624
75%        -1.61            .39
90%         -.61            .39       Variance       50.68904
95%          .39            .39       Skewness        -1.4932
99%          .39            .39       Kurtosis       5.867083

. sum gap if ri==1 & spf==1, d

                             gap
-------------------------------------------------------------
      Percentiles      Smallest
 1%       -38.04         -48.04
 5%       -18.04         -38.04
10%        -8.04         -33.04       Obs                 178
25%        -4.04         -33.04       Sum of Wgt.         178

50%        -1.04                      Mean          -3.663596
                        Largest       Std. Dev.      7.011605
75%         -.04            .96
90%         -.04           1.96       Variance       49.16261
95%          .96           1.96       Skewness      -3.561048
99%         1.96           6.96       Kurtosis       18.19928

. sum gap if ri==1 & food==1, d

                             gap
-------------------------------------------------------------
      Percentiles      Smallest
 1%       -28.61         -38.61
 5%       -18.61         -28.61
10%       -16.61         -23.61       Obs                 147
25%        -8.61         -23.61       Sum of Wgt.         147

50%        -3.61                      Mean          -6.623605
                        Largest       Std. Dev.      6.550109
75%        -1.61            .39
90%         -.61            .39       Variance       42.90392
95%          .39            .39       Skewness      -1.621386
99%          .39            .39       Kurtosis       6.577086

. restore

. 
. ** Distribution of expectations:
. sum base_pf_PP1, d

                         base_pf_PP1
-------------------------------------------------------------
no observations

. sum base_pf_RI1, d

                         base_pf_RI1
-------------------------------------------------------------
      Percentiles      Smallest
 1%            0              0
 5%            1              0
10%            1              0       Obs                 325
25%            2              0       Sum of Wgt.         325

50%            3                      Mean           5.570892
                        Largest       Std. Dev.      6.980568
75%            5             30
90%           12             30       Variance       48.72833
95%           20             45       Skewness       2.922285
99%           30             50       Kurtosis       13.26124

. 
. ** Distribution of non-updaters in shaded and unshaded quadrants (Fig 2):
. preserve

. replace gap = 1.39 - belief_food if randomscreen2==4 & food==1
(41 real changes made)

. replace gap = 1.96 - belief_spf if randomscreen2==4 & spf==1
(40 real changes made)

. if "`version'"=="interval" {
.         replace gap = 0 if spf & (belief_spf >= 1.19 & belief_spf<=3.03) /* & control==0 */
.         replace gap = 3.03 - belief_spf if spf & belief_spf>3.03 & belief_spf!=. /* & control==0 */
.         replace gap = 1.19 - belief_spf if spf & belief_spf<1.19 /* & control==0 */
.         }

. gen sensible = gap * pf_revision1 > 0 if pf_revision1!=0
(159 missing values generated)

. bysort food control: sum sensible

----------------------------------------------------------------------------------------------------------------------------------------------------------------------
-> food = 0, control = 0

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    sensible |        80        .675    .4713299          0          1

----------------------------------------------------------------------------------------------------------------------------------------------------------------------
-> food = 0, control = 1

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    sensible |        21    .5238095    .5117663          0          1

----------------------------------------------------------------------------------------------------------------------------------------------------------------------
-> food = 1, control = 0

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    sensible |        50          .7      .46291          0          1

----------------------------------------------------------------------------------------------------------------------------------------------------------------------
-> food = 1, control = 1

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
    sensible |        15          .6    .5070926          0          1


. restore

. 
. ** Changes in individual uncertainty:
. preserve

. replace gap = 1.39 - belief_food if randomscreen2==4 & food==1
(41 real changes made)

. replace gap = 1.96 - belief_spf if randomscreen2==4 & spf==1
(40 real changes made)

. if "`version'"=="interval" {
.         replace gap = 0 if spf & (belief_spf >= 1.19 & belief_spf<=3.03) /* & control==0 */
.         replace gap = 3.03 - belief_spf if spf & belief_spf>3.03 & belief_spf!=. /* & control==0 */
.         replace gap = 1.19 - belief_spf if spf & belief_spf<1.19 /* & control==0 */
.         }

. gen low_gap = abs(gap)<=2 if gap

. table low_gap control, contents(mean vara mean varb)

------------------------------
          |      control      
  low_gap |        0         1
----------+-------------------
        0 | 28.37107  17.20053
          | 16.49735  12.98348
          | 
        1 | 9.768781   7.53221
          | 6.012055  6.461976
------------------------------

. restore

. 
. ** Forecast accuracy:
. preserve

. local realized_spf = 2.93

. gen base_accurate = abs(base_pf_PP1-`realized_spf')<1 if base_pf_PP1!=.
(325 missing values generated)

. replace base_accurate = abs(base_pf_RI1-`realized_spf')<1 if base_pf_RI1!=.
(325 real changes made)

. gen end_accurate = abs(end_pf_PP1-`realized_spf')<1 if end_pf_PP1!=.
(325 missing values generated)

. replace end_accurate = abs(end_pf_RI1-`realized_spf')<1 if end_pf_RI1!=.
(325 real changes made)

. gen acc_change = end_accurate - base_accurate

. sum base_accurate if !control

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
base_accur~e |       244    .3934426    .4895177          0          1

. sum end_accurate if !control

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
end_accurate |       244    .5655738    .4967003          0          1

. restore

. 
. ** Gender differences in info sources and strength of priors:
. egen source_weight = rowtotal(b4_a- b4_g)

. gen shop_pref = b4_g / source_weight
(2 missing values generated)

. gen tv_pref = b4_a / source_weight

. gen paper_pref = b4_b / source_weight

. gen famfriends = b4_f / source_weight

. gen groceries=b3_l==7

. foreach pref in shop_pref tv_pref paper_pref famfriends groceries {
  2.         ttest `pref', by(fem)
  3.         }

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     139    .1709385     .007174    .0845798    .1567534    .1851236
       1 |     184    .1798279    .0060673    .0823007    .1678571    .1917987
---------+--------------------------------------------------------------------
combined |     323    .1760024    .0046336    .0832755    .1668865    .1851183
---------+--------------------------------------------------------------------
    diff |           -.0088894    .0093598               -.0273038    .0095249
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.9497
Ho: diff = 0                                     degrees of freedom =      321

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1715         Pr(|T| > |t|) = 0.3430          Pr(T > t) = 0.8285

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     139    .1898338    .0056255    .0663242    .1787104    .2009572
       1 |     186    .1782993     .005034    .0686548    .1683678    .1882307
---------+--------------------------------------------------------------------
combined |     325    .1832325    .0037612    .0678057    .1758331    .1906319
---------+--------------------------------------------------------------------
    diff |            .0115345    .0075869               -.0033916    .0264606
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   1.5203
Ho: diff = 0                                     degrees of freedom =      323

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9353         Pr(|T| > |t|) = 0.1294          Pr(T > t) = 0.0647

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     139     .162687    .0057446    .0677274    .1513282    .1740457
       1 |     186    .1615127    .0045791    .0624509    .1524787    .1705468
---------+--------------------------------------------------------------------
combined |     325    .1620149    .0035867    .0646605    .1549587    .1690711
---------+--------------------------------------------------------------------
    diff |            .0011742    .0072606               -.0131097    .0154582
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   0.1617
Ho: diff = 0                                     degrees of freedom =      323

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.5642         Pr(|T| > |t|) = 0.8716          Pr(T > t) = 0.4358

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     139    .1118393    .0046422    .0547302    .1026603    .1210182
       1 |     186    .1343745    .0050924    .0694511    .1243279    .1444212
---------+--------------------------------------------------------------------
combined |     325    .1247364    .0035754    .0644565    .1177024    .1317703
---------+--------------------------------------------------------------------
    diff |           -.0225352    .0071285               -.0365594   -.0085111
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -3.1613
Ho: diff = 0                                     degrees of freedom =      323

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0009         Pr(|T| > |t|) = 0.0017          Pr(T > t) = 0.9991

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     139    .1798561     .032694    .3854566    .1152101    .2445021
       1 |     186    .3333333    .0346583    .4726769    .2649569    .4017097
---------+--------------------------------------------------------------------
combined |     325    .2676923    .0245976     .443439    .2193012    .3160834
---------+--------------------------------------------------------------------
    diff |           -.1534772    .0490569               -.2499887   -.0569657
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -3.1286
Ho: diff = 0                                     degrees of freedom =      323

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0010         Pr(|T| > |t|) = 0.0019          Pr(T > t) = 0.9990

. bysort question: ttest vara, by(fem)

----------------------------------------------------------------------------------------------------------------------------------------------------------------------
-> question = RI

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     137    9.715088     2.37374    27.78392    5.020873     14.4093
       1 |     182    25.13739     3.62412    48.89205    17.98643    32.28835
---------+--------------------------------------------------------------------
combined |     319    18.51402    2.341659    41.82337    13.90692    23.12112
---------+--------------------------------------------------------------------
    diff |            -15.4223    4.658228               -24.58725   -6.257354
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -3.3108
Ho: diff = 0                                     degrees of freedom =      317

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0005         Pr(|T| > |t|) = 0.0010          Pr(T > t) = 0.9995

. bysort question: ttest varb, by(fem)

----------------------------------------------------------------------------------------------------------------------------------------------------------------------
-> question = RI

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     134    7.484503    2.338383    27.06874    2.859273    12.10973
       1 |     182    14.47077    2.264136    30.54486    10.00327    18.93827
---------+--------------------------------------------------------------------
combined |     316    11.50824    1.647232    29.28183    8.267272    14.74921
---------+--------------------------------------------------------------------
    diff |           -6.986267    3.315083               -13.50885   -.4636833
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -2.1074
Ho: diff = 0                                     degrees of freedom =      314

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0179         Pr(|T| > |t|) = 0.0359          Pr(T > t) = 0.9821

. 
. ** Non-updater regression analysis:
. preserve

. gen non_revise = .
(325 missing values generated)

. gen riXspf = ri * spf

.         foreach x in 0 10 20 {
  2.         replace non_revise = pf_revision1==0 if abs_gap>`x'/10 & !control
  3.         reg non_revise riXfood fem age55andup incunder75 highfinlit college, robust
  4.         foreach pref in shop_pref tv_pref paper_pref famfriends {
  5.                 ttest `pref', by(non_revise)
  6.                 }
  7.         replace non_revise = .
  8.         }       
(244 real changes made)

Linear regression                                      Number of obs =     243
                                                       F(  6,   236) =    1.25
                                                       Prob > F      =  0.2799
                                                       R-squared     =  0.0309
                                                       Root MSE      =  .49823

------------------------------------------------------------------------------
             |               Robust
  non_revise |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     riXfood |    .095529   .0651773     1.47   0.144    -.0328746    .2239327
         fem |  -.0585598   .0651297    -0.90   0.370    -.1868697    .0697501
  age55andup |   .0227899   .0652295     0.35   0.727    -.1057167    .1512964
  incunder75 |  -.0447326   .0685239    -0.65   0.515    -.1797293    .0902641
  highfinlit |    .122224   .0690054     1.77   0.078    -.0137212    .2581692
     college |  -.0315533   .0672745    -0.47   0.639    -.1640886     .100982
       _cons |    .447884   .0947411     4.73   0.000     .2612376    .6345303
------------------------------------------------------------------------------

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     129    .1737553    .0070143     .079667    .1598763    .1876342
       1 |     113    .1809583    .0078564    .0835148    .1653918    .1965247
---------+--------------------------------------------------------------------
combined |     242    .1771186    .0052323    .0813957    .1668117    .1874256
---------+--------------------------------------------------------------------
    diff |            -.007203    .0104991               -.0278852    .0134792
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.6861
Ho: diff = 0                                     degrees of freedom =      240

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.2467         Pr(|T| > |t|) = 0.4933          Pr(T > t) = 0.7533

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     130    .1867852    .0056782    .0647412    .1755508    .1980196
       1 |     114    .1813965    .0065461     .069893    .1684275    .1943654
---------+--------------------------------------------------------------------
combined |     244    .1842675    .0042964    .0671116    .1758046    .1927304
---------+--------------------------------------------------------------------
    diff |            .0053887    .0086221               -.0115953    .0223727
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   0.6250
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.7337         Pr(|T| > |t|) = 0.5326          Pr(T > t) = 0.2663

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     130    .1539151    .0050481    .0575569    .1439274    .1639029
       1 |     114    .1732156    .0066811    .0713345    .1599792    .1864521
---------+--------------------------------------------------------------------
combined |     244    .1629326    .0041578    .0649466    .1547427    .1711225
---------+--------------------------------------------------------------------
    diff |           -.0193005     .008258               -.0355672   -.0030337
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -2.3372
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0101         Pr(|T| > |t|) = 0.0202          Pr(T > t) = 0.9899

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     130     .134801    .0055577     .063367     .123805    .1457969
       1 |     114    .1137061    .0052773    .0563456    .1032509    .1241613
---------+--------------------------------------------------------------------
combined |     244    .1249452    .0039042    .0609853    .1172548    .1326355
---------+--------------------------------------------------------------------
    diff |            .0210949    .0077232                .0058816    .0363082
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   2.7314
Ho: diff = 0                                     degrees of freedom =      242

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9966         Pr(|T| > |t|) = 0.0068          Pr(T > t) = 0.0034
(244 real changes made, 244 to missing)
(180 real changes made)

Linear regression                                      Number of obs =     179
                                                       F(  6,   172) =    1.35
                                                       Prob > F      =  0.2372
                                                       R-squared     =  0.0469
                                                       Root MSE      =   .4914

------------------------------------------------------------------------------
             |               Robust
  non_revise |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     riXfood |   .1023903   .0742266     1.38   0.170    -.0441219    .2489026
         fem |  -.1500419   .0764639    -1.96   0.051    -.3009703    .0008866
  age55andup |   .0263682   .0746696     0.35   0.724    -.1210185     .173755
  incunder75 |  -.0304001   .0796453    -0.38   0.703    -.1876082    .1268079
  highfinlit |   .0888212     .08232     1.08   0.282    -.0736663    .2513087
     college |  -.0456427   .0778629    -0.59   0.559    -.1993325    .1080472
       _cons |   .4631668   .1115007     4.15   0.000     .2430809    .6832527
------------------------------------------------------------------------------

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     103    .1719868    .0079882    .0810712    .1561422    .1878313
       1 |      75    .1823351    .0102843    .0890644    .1618433     .202827
---------+--------------------------------------------------------------------
combined |     178     .176347    .0063291    .0844406    .1638568    .1888372
---------+--------------------------------------------------------------------
    diff |           -.0103484    .0128304               -.0356697    .0149729
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.8066
Ho: diff = 0                                     degrees of freedom =      176

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.2105         Pr(|T| > |t|) = 0.4210          Pr(T > t) = 0.7895

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     104    .1850619    .0065089    .0663781     .172153    .1979708
       1 |      76    .1725135    .0084066    .0732873    .1557667    .1892604
---------+--------------------------------------------------------------------
combined |     180    .1797637    .0051771    .0694578    .1695477    .1899797
---------+--------------------------------------------------------------------
    diff |            .0125484     .010469               -.0081109    .0332077
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   1.1986
Ho: diff = 0                                     degrees of freedom =      178

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.8839         Pr(|T| > |t|) = 0.2323          Pr(T > t) = 0.1161

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     104    .1511513    .0057217      .05835    .1398037    .1624989
       1 |      76     .182451    .0083119    .0724617    .1658928    .1990092
---------+--------------------------------------------------------------------
combined |     180    .1643667    .0049438    .0663286     .154611    .1741224
---------+--------------------------------------------------------------------
    diff |           -.0312996    .0097596                -.050559   -.0120402
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -3.2071
Ho: diff = 0                                     degrees of freedom =      178

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0008         Pr(|T| > |t|) = 0.0016          Pr(T > t) = 0.9992

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |     104    .1391215    .0062469    .0637066    .1267322    .1515109
       1 |      76    .1130542    .0065042    .0567022    .1000972    .1260112
---------+--------------------------------------------------------------------
combined |     180    .1281153    .0046243    .0620417    .1189901    .1372405
---------+--------------------------------------------------------------------
    diff |            .0260673    .0091833                .0079451    .0441895
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   2.8386
Ho: diff = 0                                     degrees of freedom =      178

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9975         Pr(|T| > |t|) = 0.0051          Pr(T > t) = 0.0025
(180 real changes made, 180 to missing)
(135 real changes made)

Linear regression                                      Number of obs =     134
                                                       F(  6,   127) =    1.56
                                                       Prob > F      =  0.1636
                                                       R-squared     =  0.0679
                                                       Root MSE      =  .48325

------------------------------------------------------------------------------
             |               Robust
  non_revise |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     riXfood |    .195228   .0851302     2.29   0.023     .0267706    .3636853
         fem |  -.0842174   .0908523    -0.93   0.356    -.2639977     .095563
  age55andup |    .035034   .0860202     0.41   0.684    -.1351845    .2052525
  incunder75 |  -.0418855   .0900865    -0.46   0.643    -.2201505    .1363795
  highfinlit |   .0864725   .0946333     0.91   0.363    -.1007897    .2737347
     college |  -.0503726    .086345    -0.58   0.561    -.2212337    .1204886
       _cons |   .3425585   .1269864     2.70   0.008     .0912754    .5938416
------------------------------------------------------------------------------

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      81     .169307    .0090437    .0813931    .1513095    .1873045
       1 |      52    .1851909    .0140424    .1012613    .1569996    .2133822
---------+--------------------------------------------------------------------
combined |     133    .1755172    .0077737     .089651    .1601401    .1908944
---------+--------------------------------------------------------------------
    diff |           -.0158839    .0159311               -.0473994    .0156317
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -0.9970
Ho: diff = 0                                     degrees of freedom =      131

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1603         Pr(|T| > |t|) = 0.3206          Pr(T > t) = 0.8397

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      82    .1817199    .0075241    .0681337    .1667493    .1966905
       1 |      53    .1723527    .0108671    .0791134    .1505463     .194159
---------+--------------------------------------------------------------------
combined |     135    .1780424    .0062397    .0724985    .1657014    .1903834
---------+--------------------------------------------------------------------
    diff |            .0093672    .0127998               -.0159504    .0346848
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   0.7318
Ho: diff = 0                                     degrees of freedom =      133

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.7672         Pr(|T| > |t|) = 0.4656          Pr(T > t) = 0.2328

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      82    .1496781    .0062265     .056383    .1372894    .1620668
       1 |      53    .1852969    .0109793    .0799309    .1632652    .2073285
---------+--------------------------------------------------------------------
combined |     135    .1636617     .005904    .0685985    .1519846    .1753389
---------+--------------------------------------------------------------------
    diff |           -.0356188     .011736               -.0588323   -.0124054
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =  -3.0350
Ho: diff = 0                                     degrees of freedom =      133

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0014         Pr(|T| > |t|) = 0.0029          Pr(T > t) = 0.9986

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
       0 |      82    .1430875    .0064189    .0581256    .1303159    .1558591
       1 |      53    .1185525     .008176    .0595223    .1021461    .1349588
---------+--------------------------------------------------------------------
combined |     135    .1334552    .0051365    .0596804    .1232962    .1436143
---------+--------------------------------------------------------------------
    diff |             .024535    .0103414                .0040801    .0449899
------------------------------------------------------------------------------
    diff = mean(0) - mean(1)                                      t =   2.3725
Ho: diff = 0                                     degrees of freedom =      133

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9904         Pr(|T| > |t|) = 0.0191          Pr(T > t) = 0.0096
(135 real changes made, 135 to missing)

. replace non_revise = pf_revision1==0 if abs_gap>1 & !control
(180 real changes made)

. restore

. 
. ** Correlations of Table 4 regressors
. corr fem college incover75 highfinlit age55andup
(obs=323)

             |      fem  college incov~75 highfi~t age55~up
-------------+---------------------------------------------
         fem |   1.0000
     college |  -0.0993   1.0000
   incover75 |  -0.0992   0.2855   1.0000
  highfinlit |  -0.1262   0.1257   0.0581   1.0000
  age55andup |   0.0494  -0.0642   0.0030  -0.0684   1.0000


. 
. ** Gap between density and PF forecasts:
. preserve

. gen base_pf_dens_gap = base_pf_RI1 - meana
(6 missing values generated)

. replace base_pf_dens_gap = base_pf_PP1 - meana if base_pf_dens_gap==.
(0 real changes made)

. count if meana==.
    6

. gen end_pf_dens_gap = end_pf_RI1 - meanb
(9 missing values generated)

. replace  end_pf_dens_gap = end_pf_PP1 - meanb if end_pf_dens_gap==.
(0 real changes made)

. count if meanb==.
    9

. gen base_pf_dens_gap2 = base_pf_RI1 - mediana
(6 missing values generated)

. replace base_pf_dens_gap2 = base_pf_PP1 - mediana if base_pf_dens_gap2==.
(0 real changes made)

. count if mediana==.
    6

. gen end_pf_dens_gap2 = end_pf_RI1 - medianb
(9 missing values generated)

. replace end_pf_dens_gap2 = end_pf_PP1 - medianb if end_pf_dens_gap2==.
(0 real changes made)

. count if medianb==.
    9

. tab base_pf_dens_gap

base_pf_den |
      s_gap |      Freq.     Percent        Cum.
------------+-----------------------------------
        -10 |          1        0.31        0.31
  -9.323057 |          1        0.31        0.63
  -8.986798 |          1        0.31        0.94
  -8.237783 |          1        0.31        1.25
  -5.587743 |          1        0.31        1.57
  -5.429499 |          1        0.31        1.88
  -5.348971 |          1        0.31        2.19
  -5.182182 |          1        0.31        2.51
         -5 |          1        0.31        2.82
  -4.582358 |          1        0.31        3.13
  -3.998669 |          1        0.31        3.45
  -3.510189 |          1        0.31        3.76
  -3.414917 |          1        0.31        4.08
  -3.402034 |          1        0.31        4.39
  -3.274693 |          1        0.31        4.70
    -3.2213 |          1        0.31        5.02
  -3.116398 |          1        0.31        5.33
  -3.043982 |          1        0.31        5.64
  -2.884361 |          1        0.31        5.96
  -2.716315 |          1        0.31        6.27
  -2.650198 |          1        0.31        6.58
   -2.42331 |          1        0.31        6.90
   -2.36729 |          1        0.31        7.21
  -2.313309 |          1        0.31        7.52
  -2.262143 |          1        0.31        7.84
  -2.179023 |          1        0.31        8.15
  -2.159967 |          1        0.31        8.46
  -2.107406 |          1        0.31        8.78
  -2.051035 |          1        0.31        9.09
  -2.046669 |          1        0.31        9.40
         -2 |          1        0.31        9.72
  -1.987562 |          1        0.31       10.03
  -1.963086 |          1        0.31       10.34
  -1.888098 |          1        0.31       10.66
   -1.71948 |          1        0.31       10.97
  -1.691754 |          1        0.31       11.29
   -1.62331 |          1        0.31       11.60
  -1.618034 |          1        0.31       11.91
  -1.617358 |          1        0.31       12.23
  -1.599607 |          1        0.31       12.54
  -1.572991 |          1        0.31       12.85
  -1.567113 |          1        0.31       13.17
  -1.486081 |          1        0.31       13.48
  -1.462475 |          1        0.31       13.79
  -1.381966 |          1        0.31       14.11
  -1.358372 |          1        0.31       14.42
  -1.344923 |          1        0.31       14.73
  -1.257098 |          1        0.31       15.05
  -1.240567 |          1        0.31       15.36
  -1.235626 |          1        0.31       15.67
  -1.235308 |          1        0.31       15.99
  -1.221133 |          1        0.31       16.30
  -1.217237 |          1        0.31       16.61
   -1.19129 |          1        0.31       16.93
  -1.190983 |          1        0.31       17.24
  -1.179964 |          1        0.31       17.55
  -1.174911 |          1        0.31       17.87
  -1.166625 |          1        0.31       18.18
  -1.166564 |          1        0.31       18.50
  -1.149837 |          1        0.31       18.81
  -1.128897 |          1        0.31       19.12
  -1.121313 |          1        0.31       19.44
  -1.075049 |          1        0.31       19.75
  -1.046669 |          1        0.31       20.06
  -1.037368 |          1        0.31       20.38
  -1.015459 |          2        0.63       21.00
         -1 |          3        0.94       21.94
  -.9972816 |          1        0.31       22.26
  -.9916453 |          1        0.31       22.57
  -.9808419 |          3        0.94       23.51
  -.9619601 |          1        0.31       23.82
  -.9468609 |          1        0.31       24.14
  -.9442911 |          1        0.31       24.45
  -.9165087 |          1        0.31       24.76
  -.8673916 |          1        0.31       25.08
  -.8137271 |          1        0.31       25.39
  -.8132937 |          1        0.31       25.71
   -.812191 |          1        0.31       26.02
  -.7621815 |          1        0.31       26.33
       -.75 |          1        0.31       26.65
   -.742111 |          1        0.31       26.96
  -.7345829 |          1        0.31       27.27
   -.723074 |          1        0.31       27.59
  -.7119927 |          1        0.31       27.90
  -.6880074 |          1        0.31       28.21
  -.6730304 |          1        0.31       28.53
  -.6616378 |          1        0.31       28.84
  -.6414132 |          1        0.31       29.15
  -.6295598 |          1        0.31       29.47
  -.5930905 |          2        0.63       30.09
  -.5814366 |          1        0.31       30.41
  -.5469182 |          2        0.63       31.03
  -.5375247 |          1        0.31       31.35
  -.5307865 |          1        0.31       31.66
        -.5 |          1        0.31       31.97
  -.4624753 |          3        0.94       32.92
  -.4230027 |          1        0.31       33.23
  -.4179411 |          1        0.31       33.54
  -.4168448 |          1        0.31       33.86
  -.4143574 |          1        0.31       34.17
  -.4053369 |          1        0.31       34.48
  -.3840704 |          1        0.31       34.80
  -.3772302 |          1        0.31       35.11
  -.3765996 |          1        0.31       35.42
  -.3591697 |          1        0.31       35.74
  -.3583419 |          1        0.31       36.05
  -.3405898 |          1        0.31       36.36
  -.3201599 |          1        0.31       36.68
  -.3196054 |          1        0.31       36.99
  -.2957478 |          1        0.31       37.30
  -.2880073 |          1        0.31       37.62
  -.2810936 |          1        0.31       37.93
  -.2775943 |          1        0.31       38.24
  -.2775307 |          1        0.31       38.56
       -.26 |          1        0.31       38.87
  -.2363321 |          1        0.31       39.18
  -.2326421 |          1        0.31       39.50
  -.2119929 |          1        0.31       39.81
  -.2068129 |          1        0.31       40.13
  -.2015791 |          1        0.31       40.44
  -.1974452 |          1        0.31       40.75
  -.1974404 |          1        0.31       41.07
  -.1952472 |          1        0.31       41.38
  -.1878091 |          3        0.94       42.32
   -.126503 |          1        0.31       42.63
  -.1136959 |          1        0.31       42.95
  -.1134979 |          1        0.31       43.26
  -.1129402 |          1        0.31       43.57
   -.110189 |          1        0.31       43.89
  -.1086645 |          1        0.31       44.20
  -.0928864 |          1        0.31       44.51
  -.0673641 |          1        0.31       44.83
  -.0466568 |          1        0.31       45.14
  -.0375247 |          1        0.31       45.45
  -.0286889 |          1        0.31       45.77
  -.0226941 |          1        0.31       46.08
  -.0154591 |          1        0.31       46.39
   -.009017 |          1        0.31       46.71
  -.0020881 |          1        0.31       47.02
  -2.22e-09 |          1        0.31       47.34
          0 |          9        2.82       50.16
   .0028193 |          1        0.31       50.47
   .0082077 |          1        0.31       50.78
   .0191581 |          2        0.63       51.41
   .0202298 |          1        0.31       51.72
   .0327857 |          1        0.31       52.04
   .0433803 |          1        0.31       52.35
   .0538626 |          1        0.31       52.66
   .0738615 |          1        0.31       52.98
   .0919006 |          1        0.31       53.29
   .0983713 |          1        0.31       53.61
   .1104563 |          1        0.31       53.92
   .1127589 |          1        0.31       54.23
   .1305606 |          1        0.31       54.55
    .163168 |          1        0.31       54.86
   .1634083 |          1        0.31       55.17
   .1798072 |          1        0.31       55.49
    .210744 |          1        0.31       55.80
   .2288862 |          1        0.31       56.11
   .2306275 |          1        0.31       56.43
   .2370681 |          1        0.31       56.74
   .2530146 |          1        0.31       57.05
   .2705609 |          1        0.31       57.37
   .2815733 |          1        0.31       57.68
   .2925053 |          1        0.31       57.99
   .3348122 |          1        0.31       58.31
   .3722285 |          1        0.31       58.62
    .375618 |          1        0.31       58.93
   .3819656 |          1        0.31       59.25
   .3968625 |          1        0.31       59.56
   .3987248 |          1        0.31       59.87
   .4162655 |          1        0.31       60.19
   .4530818 |          2        0.63       60.82
   .4665489 |          1        0.31       61.13
   .4758615 |          1        0.31       61.44
   .4845419 |          1        0.31       61.76
    .529579 |          1        0.31       62.07
   .5469182 |          1        0.31       62.38
   .5474262 |          1        0.31       62.70
   .5647149 |          1        0.31       63.01
   .6065769 |          1        0.31       63.32
   .6100283 |          1        0.31       63.64
   .6173837 |          1        0.31       63.95
   .6180339 |          2        0.63       64.58
   .6407607 |          1        0.31       64.89
   .6769152 |          1        0.31       65.20
   .6799936 |          1        0.31       65.52
   .6895013 |          1        0.31       65.83
   .6998297 |          1        0.31       66.14
   .7119927 |          1        0.31       66.46
   .7161703 |          1        0.31       66.77
   .7188776 |          1        0.31       67.08
   .7538399 |          1        0.31       67.40
   .7688274 |          1        0.31       67.71
   .7739716 |          1        0.31       68.03
   .7950568 |          1        0.31       68.34
   .7974055 |          1        0.31       68.65
   .7988856 |          1        0.31       68.97
   .8047528 |          1        0.31       69.28
   .8090169 |          1        0.31       69.59
    .809347 |          1        0.31       69.91
   .9029692 |          1        0.31       70.22
   .9091156 |          1        0.31       70.53
   .9135907 |          1        0.31       70.85
   .9494693 |          1        0.31       71.16
   .9624753 |          1        0.31       71.47
   .9647006 |          1        0.31       71.79
   .9825349 |          1        0.31       72.10
          1 |          3        0.94       73.04
   1.002626 |          1        0.31       73.35
   1.037525 |          1        0.31       73.67
   1.051112 |          1        0.31       73.98
   1.061028 |          1        0.31       74.29
   1.088264 |          1        0.31       74.61
   1.089377 |          1        0.31       74.92
   1.113612 |          1        0.31       75.24
   1.239299 |          1        0.31       75.55
    1.26576 |          1        0.31       75.86
   1.288376 |          1        0.31       76.18
   1.295864 |          1        0.31       76.49
   1.379573 |          1        0.31       76.80
   1.450189 |          1        0.31       77.12
        1.5 |          1        0.31       77.43
   1.517488 |          1        0.31       77.74
   1.578971 |          1        0.31       78.06
   1.587671 |          1        0.31       78.37
   1.653785 |          1        0.31       78.68
   1.673832 |          1        0.31       79.00
   1.711993 |          1        0.31       79.31
    1.76741 |          1        0.31       79.62
   1.902969 |          1        0.31       79.94
    1.92541 |          1        0.31       80.25
   2.090209 |          1        0.31       80.56
   2.102572 |          1        0.31       80.88
   2.254886 |          1        0.31       81.19
   2.255458 |          1        0.31       81.50
   2.299099 |          1        0.31       81.82
   2.324213 |          1        0.31       82.13
   2.490275 |          1        0.31       82.45
   2.542754 |          1        0.31       82.76
   2.582493 |          1        0.31       83.07
   2.653988 |          1        0.31       83.39
   2.722993 |          1        0.31       83.70
    2.94549 |          1        0.31       84.01
   2.976931 |          1        0.31       84.33
   3.055729 |          1        0.31       84.64
   3.142857 |          1        0.31       84.95
   3.393032 |          1        0.31       85.27
   3.436637 |          1        0.31       85.58
   3.889673 |          1        0.31       85.89
   3.906663 |          1        0.31       86.21
   3.947875 |          1        0.31       86.52
          4 |          1        0.31       86.83
   4.170963 |          1        0.31       87.15
   4.332482 |          1        0.31       87.46
   4.367932 |          1        0.31       87.77
   4.401412 |          1        0.31       88.09
    4.57281 |          1        0.31       88.40
   4.599374 |          1        0.31       88.71
          5 |          2        0.63       89.34
   5.777871 |          1        0.31       89.66
   5.784859 |          1        0.31       89.97
   5.812921 |          1        0.31       90.28
   6.379319 |          1        0.31       90.60
    6.67187 |          1        0.31       90.91
    6.70469 |          1        0.31       91.22
   6.740735 |          1        0.31       91.54
   6.836022 |          1        0.31       91.85
          8 |          1        0.31       92.16
   8.737311 |          1        0.31       92.48
   9.029656 |          1        0.31       92.79
   9.793032 |          1        0.31       93.10
    9.87855 |          1        0.31       93.42
   9.902969 |          1        0.31       93.73
    10.3163 |          1        0.31       94.04
   10.80146 |          1        0.31       94.36
   11.43592 |          1        0.31       94.67
   11.91486 |          1        0.31       94.98
   12.20572 |          1        0.31       95.30
   12.97693 |          1        0.31       95.61
   13.16969 |          1        0.31       95.92
   14.48348 |          1        0.31       96.24
   15.52033 |          1        0.31       96.55
   17.85347 |          1        0.31       96.87
   18.16173 |          1        0.31       97.18
   19.90297 |          1        0.31       97.49
   21.38121 |          1        0.31       97.81
   22.10549 |          1        0.31       98.12
   27.33819 |          1        0.31       98.43
   27.99567 |          1        0.31       98.75
   28.86777 |          1        0.31       99.06
   29.14048 |          1        0.31       99.37
   30.67545 |          1        0.31       99.69
   49.90297 |          1        0.31      100.00
------------+-----------------------------------
      Total |        319      100.00

. tab end_pf_dens_gap

end_pf_dens |
       _gap |      Freq.     Percent        Cum.
------------+-----------------------------------
        -10 |          2        0.63        0.63
  -9.987562 |          1        0.32        0.95
  -7.472136 |          1        0.32        1.27
   -6.72819 |          1        0.32        1.58
   -5.84317 |          1        0.32        1.90
  -5.266542 |          1        0.32        2.22
         -5 |          2        0.63        2.85
  -4.884361 |          1        0.32        3.16
  -4.801476 |          1        0.32        3.48
  -3.316714 |          1        0.32        3.80
    -3.2213 |          1        0.32        4.11
         -3 |          1        0.32        4.43
  -2.568528 |          1        0.32        4.75
  -2.514328 |          1        0.32        5.06
  -2.414917 |          1        0.32        5.38
  -2.406704 |          1        0.32        5.70
  -2.358372 |          1        0.32        6.01
  -2.166625 |          1        0.32        6.33
  -2.119022 |          1        0.32        6.65
  -2.107406 |          1        0.32        6.96
   -2.06111 |          1        0.32        7.28
  -2.013203 |          1        0.32        7.59
         -2 |          1        0.32        7.91
  -1.930356 |          1        0.32        8.23
  -1.875621 |          1        0.32        8.54
  -1.857874 |          1        0.32        8.86
  -1.817822 |          1        0.32        9.18
  -1.691754 |          1        0.32        9.49
  -1.614166 |          1        0.32        9.81
  -1.599607 |          2        0.63       10.44
  -1.416845 |          1        0.32       10.76
   -1.41577 |          1        0.32       11.08
  -1.374642 |          1        0.32       11.39
  -1.312637 |          1        0.32       11.71
  -1.292698 |          1        0.32       12.03
  -1.264231 |          1        0.32       12.34
  -1.236591 |          1        0.32       12.66
  -1.235626 |          1        0.32       12.97
  -1.235308 |          1        0.32       13.29
  -1.223626 |          1        0.32       13.61
  -1.217237 |          1        0.32       13.92
  -1.199833 |          1        0.32       14.24
   -1.19129 |          1        0.32       14.56
  -1.166464 |          1        0.32       14.87
  -1.149167 |          1        0.32       15.19
   -1.14159 |          1        0.32       15.51
  -1.103327 |          1        0.32       15.82
  -1.076137 |          1        0.32       16.14
  -1.075049 |          1        0.32       16.46
  -1.046669 |          1        0.32       16.77
  -1.032797 |          1        0.32       17.09
  -1.024731 |          1        0.32       17.41
         -1 |          4        1.27       18.67
  -.9973738 |          1        0.32       18.99
  -.9808419 |          1        0.32       19.30
  -.9730303 |          1        0.32       19.62
  -.9724979 |          1        0.32       19.94
  -.9619601 |          1        0.32       20.25
  -.9377465 |          1        0.32       20.57
  -.8955483 |          1        0.32       20.89
  -.8765864 |          1        0.32       21.20
  -.8624932 |          1        0.32       21.52
   -.812191 |          2        0.63       22.15
  -.8090169 |          1        0.32       22.47
  -.7935739 |          2        0.63       23.10
  -.7750652 |          1        0.32       23.42
  -.7607012 |          1        0.32       23.73
  -.7499309 |          1        0.32       24.05
   -.742111 |          1        0.32       24.37
  -.7345829 |          1        0.32       24.68
   -.723074 |          1        0.32       25.00
  -.7119927 |          1        0.32       25.32
  -.6917543 |          1        0.32       25.63
  -.6878091 |          1        0.32       25.95
  -.6871858 |          1        0.32       26.27
  -.6701527 |          1        0.32       26.58
  -.6560621 |          1        0.32       26.90
  -.6414132 |          1        0.32       27.22
  -.6321157 |          1        0.32       27.53
  -.6136959 |          2        0.63       28.16
  -.6044995 |          1        0.32       28.48
   -.572417 |          1        0.32       28.80
  -.5469182 |          2        0.63       29.43
  -.5375247 |          1        0.32       29.75
  -.5142376 |          1        0.32       30.06
  -.4995265 |          1        0.32       30.38
  -.4808419 |          1        0.32       30.70
  -.4655429 |          1        0.32       31.01
  -.4624753 |          1        0.32       31.33
  -.4530818 |          1        0.32       31.65
  -.3957162 |          1        0.32       31.96
  -.3561745 |          1        0.32       32.28
  -.3440948 |          1        0.32       32.59
  -.3293586 |          1        0.32       32.91
  -.3225048 |          1        0.32       33.23
  -.3060246 |          1        0.32       33.54
  -.2880073 |          4        1.27       34.81
  -.2879641 |          1        0.32       35.13
  -.2754025 |          1        0.32       35.44
   -.261728 |          1        0.32       35.76
  -.2610965 |          1        0.32       36.08
  -.2239394 |          1        0.32       36.39
  -.2192502 |          1        0.32       36.71
  -.2101364 |          1        0.32       37.03
  -.1998371 |          1        0.32       37.34
  -.1974452 |          1        0.32       37.66
  -.1952477 |          1        0.32       37.97
  -.1878091 |          3        0.95       38.92
   -.187809 |          1        0.32       39.24
  -.1811564 |          1        0.32       39.56
  -.1761985 |          1        0.32       39.87
  -.1625218 |          1        0.32       40.19
   -.157927 |          1        0.32       40.51
  -.1192765 |          1        0.32       40.82
  -.1134979 |          1        0.32       41.14
  -.1129402 |          1        0.32       41.46
  -.0719876 |          1        0.32       41.77
  -.0469182 |          1        0.32       42.09
  -.0375247 |          1        0.32       42.41
  -.0154591 |          1        0.32       42.72
   -.009017 |          1        0.32       43.04
  -.0020881 |          1        0.32       43.35
  -2.22e-09 |          1        0.32       43.67
          0 |         14        4.43       48.10
   .0006442 |          1        0.32       48.42
   .0191581 |          1        0.32       48.73
   .0202298 |          1        0.32       49.05
   .0300534 |          1        0.32       49.37
   .0327857 |          2        0.63       50.00
   .0406924 |          1        0.32       50.32
   .0433803 |          1        0.32       50.63
   .0469631 |          1        0.32       50.95
   .0503377 |          1        0.32       51.27
   .0533142 |          1        0.32       51.58
   .0538626 |          2        0.63       52.22
   .0672071 |          1        0.32       52.53
   .0862678 |          1        0.32       52.85
    .092545 |          1        0.32       53.16
   .0976237 |          1        0.32       53.48
   .0983713 |          2        0.63       54.11
   .1081557 |          1        0.32       54.43
   .1127589 |          1        0.32       54.75
   .1305606 |          1        0.32       55.06
    .163168 |          1        0.32       55.38
   .2064261 |          2        0.63       56.01
   .2113273 |          1        0.32       56.33
   .2249348 |          1        0.32       56.65
   .2315617 |          1        0.32       56.96
   .2370681 |          1        0.32       57.28
   .2429018 |          1        0.32       57.59
   .2656802 |          1        0.32       57.91
   .2808121 |          2        0.63       58.54
   .2860346 |          1        0.32       58.86
   .2976676 |          1        0.32       59.18
   .3301845 |          1        0.32       59.49
   .3526423 |          1        0.32       59.81
   .3678843 |          1        0.32       60.13
    .370002 |          1        0.32       60.44
   .3722285 |          1        0.32       60.76
    .389811 |          1        0.32       61.08
   .3955005 |          1        0.32       61.39
   .3980577 |          1        0.32       61.71
    .402458 |          1        0.32       62.03
   .4530818 |          2        0.63       62.66
   .4565569 |          1        0.32       62.97
   .4624753 |          1        0.32       63.29
   .4780896 |          1        0.32       63.61
   .4866768 |          1        0.32       63.92
         .5 |          2        0.63       64.56
   .5000001 |          1        0.32       64.87
   .5081614 |          1        0.32       65.19
   .5375247 |          2        0.63       65.82
   .5452344 |          1        0.32       66.14
   .5469182 |          2        0.63       66.77
   .5616026 |          1        0.32       67.09
   .5784037 |          1        0.32       67.41
   .5990288 |          1        0.32       67.72
   .6065769 |          1        0.32       68.04
   .6100283 |          1        0.32       68.35
   .6104563 |          1        0.32       68.67
   .6180339 |          1        0.32       68.99
   .6360432 |          1        0.32       69.30
   .6512642 |          1        0.32       69.62
   .6537855 |          1        0.32       69.94
   .6909831 |          1        0.32       70.25
   .7119927 |          1        0.32       70.57
    .719188 |          1        0.32       70.89
   .7979983 |          1        0.32       71.20
   .8034034 |          1        0.32       71.52
   .8047528 |          1        0.32       71.84
   .8090131 |          1        0.32       72.15
   .8090169 |          1        0.32       72.47
   .8121909 |          2        0.63       73.10
   .8370707 |          1        0.32       73.42
   .8474405 |          1        0.32       73.73
   .8554537 |          1        0.32       74.05
   .9029692 |          1        0.32       74.37
   .9043168 |          1        0.32       74.68
   .9061637 |          1        0.32       75.00
   .9380488 |          1        0.32       75.32
   .9940587 |          1        0.32       75.63
          1 |          3        0.95       76.58
          1 |          1        0.32       76.90
   1.000001 |          1        0.32       77.22
   1.019158 |          1        0.32       77.53
   1.027727 |          1        0.32       77.85
   1.036589 |          1        0.32       78.16
   1.059032 |          1        0.32       78.48
   1.074134 |          1        0.32       78.80
   1.108238 |          1        0.32       79.11
   1.146271 |          1        0.32       79.43
   1.211085 |          1        0.32       79.75
   1.211993 |          1        0.32       80.06
   1.244188 |          1        0.32       80.38
   1.285028 |          1        0.32       80.70
   1.307564 |          1        0.32       81.01
   1.317973 |          1        0.32       81.33
   1.324769 |          1        0.32       81.65
   1.379043 |          1        0.32       81.96
    1.38189 |          1        0.32       82.28
   1.433361 |          1        0.32       82.59
   1.462475 |          1        0.32       82.91
   1.578971 |          1        0.32       83.23
    1.63271 |          1        0.32       83.54
    1.64083 |          1        0.32       83.86
   1.668137 |          1        0.32       84.18
   1.673832 |          1        0.32       84.49
   1.711993 |          1        0.32       84.81
   1.902969 |          2        0.63       85.44
   1.911114 |          1        0.32       85.76
   1.984541 |          1        0.32       86.08
          2 |          1        0.32       86.39
   2.000001 |          1        0.32       86.71
   2.008355 |          1        0.32       87.03
   2.071964 |          1        0.32       87.34
   2.090209 |          1        0.32       87.66
   2.116089 |          1        0.32       87.97
   2.280949 |          1        0.32       88.29
   2.295864 |          1        0.32       88.61
   2.727027 |          1        0.32       88.92
   2.816901 |          1        0.32       89.24
          3 |          1        0.32       89.56
   3.055729 |          1        0.32       89.87
     3.3086 |          1        0.32       90.19
   3.368086 |          1        0.32       90.51
   3.521715 |          1        0.32       90.82
   3.582059 |          1        0.32       91.14
   3.801783 |          1        0.32       91.46
   3.889673 |          1        0.32       91.77
   4.015435 |          1        0.32       92.09
   4.261664 |          1        0.32       92.41
   4.304763 |          1        0.32       92.72
   4.367885 |          1        0.32       93.04
   4.440812 |          1        0.32       93.35
   4.756403 |          1        0.32       93.67
   4.828108 |          1        0.32       93.99
   4.878368 |          1        0.32       94.30
   5.163408 |          1        0.32       94.62
   5.170513 |          1        0.32       94.94
   5.497416 |          1        0.32       95.25
   5.712245 |          1        0.32       95.57
    6.26576 |          1        0.32       95.89
    6.26888 |          1        0.32       96.20
   6.350212 |          1        0.32       96.52
   6.443104 |          1        0.32       96.84
          7 |          1        0.32       97.15
   9.070368 |          1        0.32       97.47
   9.802817 |          1        0.32       97.78
   9.840514 |          1        0.32       98.10
    10.3833 |          1        0.32       98.42
   11.20732 |          1        0.32       98.73
   12.90297 |          1        0.32       99.05
   15.67544 |          2        0.63       99.68
   34.90297 |          1        0.32      100.00
------------+-----------------------------------
      Total |        316      100.00

. sum base_pf_dens_gap end_pf_dens_gap

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
base_pf_de~p |       319     1.45928    5.989639        -10   49.90297
end_pf_den~p |       316    .5057486    3.451994        -10   34.90297

. tab base_pf_dens_gap2

base_pf_den |
     s_gap2 |      Freq.     Percent        Cum.
------------+-----------------------------------
        -10 |          1        0.31        0.31
  -9.757781 |          1        0.31        0.63
  -9.445333 |          1        0.31        0.94
   -7.33689 |          1        0.31        1.25
  -7.194073 |          1        0.31        1.57
   -6.22834 |          1        0.31        1.88
  -5.292004 |          1        0.31        2.19
         -5 |          1        0.31        2.51
  -4.095685 |          1        0.31        2.82
   -3.81988 |          1        0.31        3.13
  -3.762981 |          1        0.31        3.45
  -3.671712 |          1        0.31        3.76
  -3.609202 |          1        0.31        4.08
  -3.414135 |          1        0.31        4.39
  -3.187449 |          1        0.31        4.70
  -2.994953 |          1        0.31        5.02
  -2.659754 |          1        0.31        5.33
  -2.598338 |          1        0.31        5.64
  -2.546517 |          1        0.31        5.96
   -2.35058 |          1        0.31        6.27
  -2.223725 |          1        0.31        6.58
  -2.197112 |          1        0.31        6.90
  -2.179827 |          1        0.31        7.21
  -2.116518 |          1        0.31        7.52
  -2.093211 |          1        0.31        7.84
  -2.052653 |          1        0.31        8.15
   -2.05041 |          1        0.31        8.46
  -2.025753 |          1        0.31        8.78
  -2.016717 |          1        0.31        9.09
         -2 |          1        0.31        9.40
   -1.99753 |          1        0.31        9.72
  -1.955655 |          1        0.31       10.03
  -1.951488 |          1        0.31       10.34
  -1.915169 |          1        0.31       10.66
  -1.882538 |          1        0.31       10.97
  -1.649209 |          1        0.31       11.29
  -1.618898 |          1        0.31       11.60
  -1.618034 |          1        0.31       11.91
  -1.596719 |          1        0.31       12.23
  -1.564191 |          1        0.31       12.54
  -1.477407 |          1        0.31       12.85
  -1.462475 |          1        0.31       13.17
  -1.423725 |          1        0.31       13.48
  -1.417243 |          1        0.31       13.79
  -1.381966 |          1        0.31       14.11
  -1.255735 |          1        0.31       14.42
  -1.225613 |          1        0.31       14.73
  -1.220621 |          1        0.31       15.05
  -1.190983 |          1        0.31       15.36
  -1.167796 |          1        0.31       15.67
  -1.137557 |          1        0.31       15.99
  -1.120125 |          1        0.31       16.30
  -1.095543 |          1        0.31       16.61
  -1.075049 |          1        0.31       16.93
  -1.054591 |          1        0.31       17.24
  -1.030039 |          1        0.31       17.55
  -1.005141 |          1        0.31       17.87
         -1 |          4        1.25       19.12
  -.9522696 |          3        0.94       20.06
  -.9518759 |          2        0.63       20.69
   -.951488 |          1        0.31       21.00
  -.9426603 |          1        0.31       21.32
  -.9191842 |          1        0.31       21.63
  -.9175758 |          1        0.31       21.94
  -.9134599 |          1        0.31       22.26
  -.9035991 |          1        0.31       22.57
  -.8946247 |          1        0.31       22.88
  -.8851435 |          1        0.31       23.20
  -.8662915 |          1        0.31       23.51
   -.812191 |          1        0.31       23.82
  -.8098021 |          1        0.31       24.14
  -.8027401 |          1        0.31       24.45
  -.7960544 |          1        0.31       24.76
  -.7500277 |          1        0.31       25.08
       -.75 |          1        0.31       25.39
  -.7403698 |          1        0.31       25.71
  -.7234042 |          1        0.31       26.02
  -.7127016 |          1        0.31       26.33
  -.7119927 |          1        0.31       26.65
  -.6880074 |          1        0.31       26.96
  -.6815913 |          1        0.31       27.27
  -.6798892 |          1        0.31       27.59
  -.6350851 |          1        0.31       27.90
  -.6335182 |          1        0.31       28.21
  -.6313124 |          1        0.31       28.53
  -.6194565 |          1        0.31       28.84
  -.5651903 |          1        0.31       29.15
  -.5553842 |          1        0.31       29.47
  -.5543437 |          1        0.31       29.78
  -.5483224 |          1        0.31       30.09
  -.5469182 |          2        0.63       30.72
  -.5453572 |          1        0.31       31.03
  -.5375247 |          1        0.31       31.35
  -.5031943 |          1        0.31       31.66
        -.5 |          1        0.31       31.97
  -.4789619 |          1        0.31       32.29
  -.4624753 |          3        0.94       33.23
   -.445075 |          2        0.63       33.86
  -.3864622 |          1        0.31       34.17
  -.3834636 |          1        0.31       34.48
  -.3527169 |          1        0.31       34.80
  -.3351889 |          1        0.31       35.11
    -.32639 |          1        0.31       35.42
  -.3128123 |          1        0.31       35.74
  -.3036624 |          1        0.31       36.05
  -.2895222 |          1        0.31       36.36
  -.2880073 |          1        0.31       36.68
  -.2633457 |          1        0.31       36.99
       -.26 |          1        0.31       37.30
  -.2556958 |          1        0.31       37.62
  -.2483938 |          1        0.31       37.93
  -.2387288 |          1        0.31       38.24
  -.2119929 |          1        0.31       38.56
  -.2049079 |          1        0.31       38.87
  -.2015791 |          1        0.31       39.18
  -.1952472 |          1        0.31       39.50
  -.1878091 |          3        0.94       40.44
  -.1795244 |          1        0.31       40.75
  -.1747549 |          1        0.31       41.07
  -.1715965 |          1        0.31       41.38
  -.1616936 |          1        0.31       41.69
  -.1331317 |          1        0.31       42.01
  -.1331315 |          1        0.31       42.32
  -.1127111 |          1        0.31       42.63
  -.0774071 |          1        0.31       42.95
  -.0375745 |          1        0.31       43.26
  -.0375247 |          1        0.31       43.57
  -.0340834 |          1        0.31       43.89
  -.0215921 |          1        0.31       44.20
  -.0143466 |          1        0.31       44.51
   -.009017 |          1        0.31       44.83
  -.0077164 |          1        0.31       45.14
          0 |         12        3.76       48.90
   .0070217 |          1        0.31       49.22
   .0156064 |          1        0.31       49.53
   .0266417 |          1        0.31       49.84
   .0350406 |          1        0.31       50.16
   .0363393 |          1        0.31       50.47
   .0407543 |          1        0.31       50.78
   .0477017 |          1        0.31       51.10
   .0477304 |          2        0.63       51.72
   .0481241 |          1        0.31       52.04
    .067454 |          1        0.31       52.35
    .070021 |          1        0.31       52.66
   .0737091 |          1        0.31       52.98
   .0780818 |          1        0.31       53.29
   .1106623 |          1        0.31       53.61
    .118309 |          1        0.31       53.92
   .1244781 |          1        0.31       54.23
   .1440226 |          1        0.31       54.55
   .1518874 |          1        0.31       54.86
   .1565304 |          1        0.31       55.17
   .1903874 |          1        0.31       55.49
    .234098 |          1        0.31       55.80
   .2347105 |          1        0.31       56.11
   .2381479 |          1        0.31       56.43
   .2565105 |          1        0.31       56.74
    .286906 |          1        0.31       57.05
    .315811 |          1        0.31       57.37
   .3533559 |          1        0.31       57.68
   .3539746 |          1        0.31       57.99
    .375618 |          1        0.31       58.31
   .3819656 |          1        0.31       58.62
   .4137731 |          1        0.31       58.93
   .4138461 |          1        0.31       59.25
    .424611 |          1        0.31       59.56
   .4273925 |          1        0.31       59.87
   .4370852 |          1        0.31       60.19
   .4530818 |          2        0.63       60.82
   .4590381 |          1        0.31       61.13
   .4830303 |          1        0.31       61.44
   .4999187 |          1        0.31       61.76
   .5469182 |          1        0.31       62.07
   .5477409 |          1        0.31       62.38
   .5479803 |          1        0.31       62.70
   .5674182 |          1        0.31       63.01
   .5894786 |          1        0.31       63.32
   .5956879 |          1        0.31       63.64
   .5986648 |          1        0.31       63.95
   .6101494 |          1        0.31       64.26
   .6180339 |          2        0.63       64.89
   .7068813 |          1        0.31       65.20
   .7119927 |          1        0.31       65.52
    .733839 |          1        0.31       65.83
   .7363539 |          1        0.31       66.14
   .7372576 |          1        0.31       66.46
   .7548428 |          1        0.31       66.77
   .7556705 |          1        0.31       67.08
   .7586739 |          1        0.31       67.40
   .7726612 |          1        0.31       67.71
   .7889693 |          1        0.31       68.03
   .7945302 |          1        0.31       68.34
   .8047528 |          1        0.31       68.65
   .8090169 |          1        0.31       68.97
   .8234656 |          1        0.31       69.28
   .9076388 |          1        0.31       69.59
   .9146494 |          1        0.31       69.91
   .9327241 |          1        0.31       70.22
   .9404025 |          1        0.31       70.53
   .9528971 |          1        0.31       70.85
   .9624753 |          1        0.31       71.16
   .9804956 |          1        0.31       71.47
   .9833593 |          1        0.31       71.79
          1 |          3        0.94       72.73
   1.009282 |          1        0.31       73.04
   1.030602 |          1        0.31       73.35
   1.035521 |          1        0.31       73.67
   1.037525 |          1        0.31       73.98
   1.045246 |          1        0.31       74.29
   1.052857 |          1        0.31       74.61
   1.053078 |          1        0.31       74.92
   1.073868 |          1        0.31       75.24
   1.143516 |          1        0.31       75.55
   1.170224 |          1        0.31       75.86
   1.174276 |          1        0.31       76.18
   1.257029 |          1        0.31       76.49
   1.293351 |          1        0.31       76.80
   1.355998 |          1        0.31       77.12
   1.375868 |          1        0.31       77.43
   1.379986 |          1        0.31       77.74
    1.41956 |          1        0.31       78.06
   1.480008 |          1        0.31       78.37
   1.482158 |          1        0.31       78.68
        1.5 |          1        0.31       79.00
   1.602243 |          1        0.31       79.31
   1.629955 |          1        0.31       79.62
   1.663989 |          1        0.31       79.94
   1.711993 |          1        0.31       80.25
   1.751568 |          1        0.31       80.56
   1.774231 |          1        0.31       80.88
   2.045246 |          1        0.31       81.19
    2.05427 |          1        0.31       81.50
   2.059016 |          1        0.31       81.82
   2.286496 |          1        0.31       82.13
   2.616261 |          1        0.31       82.45
   2.637015 |          1        0.31       82.76
   2.700286 |          1        0.31       83.07
   2.876683 |          1        0.31       83.39
          3 |          1        0.31       83.70
   3.036462 |          1        0.31       84.01
   3.055729 |          1        0.31       84.33
   3.142857 |          1        0.31       84.64
    3.21514 |          1        0.31       84.95
   3.357055 |          1        0.31       85.27
   3.491388 |          1        0.31       85.58
   3.612233 |          1        0.31       85.89
   3.615162 |          1        0.31       86.21
          4 |          1        0.31       86.52
   4.039225 |          1        0.31       86.83
   4.343797 |          1        0.31       87.15
   4.521832 |          1        0.31       87.46
   4.542677 |          1        0.31       87.77
   4.860994 |          1        0.31       88.09
   4.925896 |          1        0.31       88.40
          5 |          2        0.63       89.03
   5.307984 |          1        0.31       89.34
   5.374727 |          1        0.31       89.66
   5.613784 |          1        0.31       89.97
   6.608578 |          1        0.31       90.28
   6.836022 |          1        0.31       90.60
    7.32396 |          1        0.31       90.91
   7.936744 |          1        0.31       91.22
          8 |          1        0.31       91.54
   8.031191 |          1        0.31       91.85
   8.079985 |          1        0.31       92.16
   8.129813 |          1        0.31       92.48
   8.508043 |          1        0.31       92.79
   8.988444 |          1        0.31       93.10
   9.973679 |          1        0.31       93.42
   10.04525 |          1        0.31       93.73
   10.24068 |          1        0.31       94.04
   10.62138 |          1        0.31       94.36
   11.23419 |          1        0.31       94.67
   11.58575 |          1        0.31       94.98
    11.6186 |          1        0.31       95.30
         13 |          1        0.31       95.61
   13.98156 |          1        0.31       95.92
   16.87288 |          1        0.31       96.24
   17.78453 |          1        0.31       96.55
   17.86004 |          1        0.31       96.87
   19.19776 |          1        0.31       97.18
   20.04525 |          1        0.31       97.49
   21.45662 |          1        0.31       97.81
   21.65118 |          1        0.31       98.12
   26.74667 |          1        0.31       98.43
   27.63198 |          1        0.31       98.75
   28.61103 |          1        0.31       99.06
   29.77257 |          1        0.31       99.37
   30.67545 |          1        0.31       99.69
   50.04525 |          1        0.31      100.00
------------+-----------------------------------
      Total |        319      100.00

. tab end_pf_dens_gap2

end_pf_dens |
      _gap2 |      Freq.     Percent        Cum.
------------+-----------------------------------
  -10.65975 |          1        0.32        0.32
        -10 |          2        0.63        0.95
  -7.472136 |          1        0.32        1.27
  -6.498382 |          1        0.32        1.58
  -6.334432 |          1        0.32        1.90
  -6.004538 |          1        0.32        2.22
  -5.769951 |          1        0.32        2.53
  -5.187449 |          1        0.32        2.85
         -5 |          2        0.63        3.48
         -3 |          2        0.63        4.11
  -2.994953 |          1        0.32        4.43
  -2.671712 |          1        0.32        4.75
  -2.618898 |          1        0.32        5.06
  -2.540733 |          1        0.32        5.38
  -2.431842 |          1        0.32        5.70
  -2.067066 |          1        0.32        6.01
  -2.025753 |          1        0.32        6.33
         -2 |          1        0.32        6.65
  -1.955655 |          1        0.32        6.96
  -1.903069 |          1        0.32        7.28
  -1.831282 |          1        0.32        7.59
  -1.649209 |          2        0.63        8.23
  -1.634995 |          1        0.32        8.54
  -1.634909 |          1        0.32        8.86
  -1.596719 |          1        0.32        9.18
  -1.585271 |          1        0.32        9.49
  -1.566236 |          1        0.32        9.81
  -1.555384 |          1        0.32       10.13
  -1.537735 |          1        0.32       10.44
  -1.493188 |          1        0.32       10.76
  -1.439476 |          1        0.32       11.08
  -1.427657 |          1        0.32       11.39
  -1.417243 |          1        0.32       11.71
  -1.404352 |          1        0.32       12.03
  -1.383464 |          1        0.32       12.34
  -1.326808 |          1        0.32       12.66
  -1.264231 |          1        0.32       12.97
   -1.17204 |          1        0.32       13.29
  -1.167796 |          1        0.32       13.61
  -1.137557 |          1        0.32       13.92
  -1.134293 |          1        0.32       14.24
   -1.12872 |          1        0.32       14.56
   -1.10274 |          1        0.32       14.87
  -1.075049 |          1        0.32       15.19
  -1.054591 |          1        0.32       15.51
         -1 |          4        1.27       16.77
  -.9750423 |          1        0.32       17.09
   -.969398 |          1        0.32       17.41
  -.9556546 |          1        0.32       17.72
  -.9537494 |          1        0.32       18.04
  -.9522696 |          1        0.32       18.35
   -.951488 |          1        0.32       18.67
  -.9420571 |          1        0.32       18.99
  -.8934367 |          1        0.32       19.30
    -.88956 |          1        0.32       19.62
  -.8750191 |          1        0.32       19.94
  -.8257241 |          1        0.32       20.25
   -.812191 |          2        0.63       20.89
  -.8112097 |          1        0.32       21.20
  -.8090169 |          1        0.32       21.52
  -.7716649 |          1        0.32       21.84
  -.7536856 |          1        0.32       22.15
  -.7449098 |          2        0.63       22.78
  -.7234042 |          1        0.32       23.10
  -.7146988 |          1        0.32       23.42
  -.7119927 |          1        0.32       23.73
  -.6878091 |          1        0.32       24.05
  -.6798892 |          1        0.32       24.37
  -.6783295 |          1        0.32       24.68
  -.6321157 |          1        0.32       25.00
  -.6312127 |          1        0.32       25.32
  -.5950899 |          1        0.32       25.63
  -.5718479 |          1        0.32       25.95
  -.5543437 |          1        0.32       26.27
  -.5469182 |          2        0.63       26.90
  -.5375745 |          2        0.63       27.53
  -.5375247 |          1        0.32       27.85
  -.5105078 |          1        0.32       28.16
  -.5067863 |          1        0.32       28.48
        -.5 |          1        0.32       28.80
  -.4789619 |          1        0.32       29.11
  -.4624753 |          1        0.32       29.43
  -.4530818 |          1        0.32       29.75
  -.4522696 |          1        0.32       30.06
  -.4451427 |          1        0.32       30.38
   -.375886 |          1        0.32       30.70
  -.3581326 |          1        0.32       31.01
  -.3339872 |          1        0.32       31.33
  -.3297133 |          1        0.32       31.65
   -.314255 |          1        0.32       31.96
  -.3092194 |          1        0.32       32.28
  -.3068383 |          1        0.32       32.59
  -.2921069 |          1        0.32       32.91
  -.2880073 |          4        1.27       34.18
  -.2827001 |          1        0.32       34.49
  -.2772493 |          1        0.32       34.81
  -.2643118 |          1        0.32       35.13
  -.2475606 |          1        0.32       35.44
  -.2176123 |          1        0.32       35.76
  -.2026165 |          1        0.32       36.08
  -.1952477 |          1        0.32       36.39
  -.1917531 |          1        0.32       36.71
  -.1878091 |          3        0.95       37.66
   -.187809 |          1        0.32       37.97
  -.1795244 |          1        0.32       38.29
  -.1661649 |          1        0.32       38.61
  -.1621699 |          1        0.32       38.92
  -.1455429 |          1        0.32       39.24
  -.1331315 |          1        0.32       39.56
  -.1212449 |          1        0.32       39.87
  -.1142721 |          1        0.32       40.19
  -.0637829 |          1        0.32       40.51
  -.0528626 |          1        0.32       40.82
  -.0469182 |          1        0.32       41.14
  -.0375247 |          1        0.32       41.46
  -.0143466 |          1        0.32       41.77
   -.009017 |          1        0.32       42.09
          0 |         16        5.06       47.15
   .0030543 |          1        0.32       47.47
   .0054178 |          1        0.32       47.78
   .0107082 |          1        0.32       48.10
   .0266417 |          1        0.32       48.42
   .0476177 |          1        0.32       48.73
   .0477304 |          1        0.32       49.05
   .0481241 |          1        0.32       49.37
   .0526569 |          1        0.32       49.68
   .0540397 |          1        0.32       50.00
   .0548166 |          1        0.32       50.32
    .067454 |          1        0.32       50.63
    .075467 |          1        0.32       50.95
   .0976237 |          1        0.32       51.27
   .1106623 |          2        0.63       51.90
   .1178576 |          1        0.32       52.22
    .118309 |          1        0.32       52.53
   .1212438 |          1        0.32       52.85
   .1244781 |          1        0.32       53.16
   .1356764 |          1        0.32       53.48
   .1370368 |          1        0.32       53.80
   .1440226 |          1        0.32       54.11
   .1702244 |          1        0.32       54.43
   .2013917 |          1        0.32       54.75
   .2283351 |          1        0.32       55.06
    .245888 |          1        0.32       55.38
   .2477008 |          1        0.32       55.70
   .2550902 |          2        0.63       56.33
   .2565105 |          1        0.32       56.65
   .2743866 |          1        0.32       56.96
   .2787519 |          1        0.32       57.28
   .2808121 |          2        0.63       57.91
   .2826161 |          1        0.32       58.23
   .3033302 |          1        0.32       58.54
   .3488529 |          1        0.32       58.86
   .3678843 |          1        0.32       59.18
   .3750601 |          1        0.32       59.49
   .3872889 |          1        0.32       59.81
    .389805 |          1        0.32       60.13
    .392158 |          1        0.32       60.44
   .3982257 |          1        0.32       60.76
   .4138461 |          2        0.63       61.39
   .4208249 |          1        0.32       61.71
   .4263287 |          1        0.32       62.03
   .4370852 |          2        0.63       62.66
   .4457822 |          1        0.32       62.97
   .4530818 |          2        0.63       63.61
   .4624753 |          1        0.32       63.92
   .4758821 |          1        0.32       64.24
   .4786229 |          1        0.32       64.56
   .4867259 |          1        0.32       64.87
   .4999187 |          1        0.32       65.19
         .5 |          3        0.95       66.14
   .5375247 |          2        0.63       66.77
   .5469182 |          2        0.63       67.41
   .5595419 |          1        0.32       67.72
    .568801 |          1        0.32       68.04
   .5737091 |          1        0.32       68.35
   .6180339 |          1        0.32       68.67
   .6382093 |          1        0.32       68.99
   .6411021 |          1        0.32       69.30
   .6882687 |          1        0.32       69.62
   .6909831 |          1        0.32       69.94
   .7068813 |          1        0.32       70.25
   .7119927 |          1        0.32       70.57
    .719188 |          1        0.32       70.89
   .7227507 |          1        0.32       71.20
   .7768848 |          1        0.32       71.52
   .7912552 |          1        0.32       71.84
    .793129 |          1        0.32       72.15
   .8047528 |          1        0.32       72.47
   .8090169 |          1        0.32       72.78
   .8121909 |          2        0.63       73.42
   .8364789 |          1        0.32       73.73
   .8479342 |          1        0.32       74.05
    .885113 |          1        0.32       74.37
   .9061637 |          1        0.32       74.68
   .9076388 |          1        0.32       75.00
   .9146494 |          1        0.32       75.32
   .9278018 |          1        0.32       75.63
   .9293137 |          1        0.32       75.95
   .9767632 |          1        0.32       76.27
   .9814183 |          1        0.32       76.58
          1 |          5        1.58       78.16
   1.045246 |          1        0.32       78.48
    1.04773 |          1        0.32       78.80
   1.065284 |          1        0.32       79.11
   1.074859 |          1        0.32       79.43
   1.182224 |          1        0.32       79.75
   1.183962 |          1        0.32       80.06
   1.208874 |          1        0.32       80.38
   1.211993 |          1        0.32       80.70
   1.212936 |          1        0.32       81.01
   1.247991 |          1        0.32       81.33
   1.337549 |          1        0.32       81.65
   1.384189 |          1        0.32       81.96
   1.462475 |          1        0.32       82.28
   1.500019 |          1        0.32       82.59
   1.577806 |          1        0.32       82.91
   1.602243 |          1        0.32       83.23
   1.711993 |          1        0.32       83.54
   1.950415 |          1        0.32       83.86
   1.994859 |          1        0.32       84.18
          2 |          2        0.63       84.81
   2.035041 |          1        0.32       85.13
   2.045246 |          2        0.63       85.76
   2.048124 |          1        0.32       86.08
   2.059016 |          1        0.32       86.39
   2.071957 |          1        0.32       86.71
   2.076388 |          1        0.32       87.03
   2.154109 |          1        0.32       87.34
   2.313917 |          1        0.32       87.66
   2.482158 |          1        0.32       87.97
   2.522593 |          1        0.32       88.29
   2.860928 |          1        0.32       88.61
   2.876683 |          1        0.32       88.92
   2.911584 |          1        0.32       89.24
   2.970689 |          1        0.32       89.56
          3 |          1        0.32       89.87
   3.049055 |          1        0.32       90.19
   3.055729 |          1        0.32       90.51
   3.194383 |          1        0.32       90.82
   3.708118 |          1        0.32       91.14
   3.826478 |          1        0.32       91.46
   3.835852 |          1        0.32       91.77
   3.992284 |          1        0.32       92.09
   4.039225 |          1        0.32       92.41
   4.131917 |          1        0.32       92.72
   4.161896 |          1        0.32       93.04
   4.273609 |          1        0.32       93.35
   4.367885 |          1        0.32       93.67
   4.777671 |          1        0.32       93.99
   4.904396 |          1        0.32       94.30
    5.06983 |          1        0.32       94.62
   5.244784 |          1        0.32       94.94
    5.54798 |          1        0.32       95.25
   5.913167 |          1        0.32       95.57
   6.236609 |          1        0.32       95.89
   6.592893 |          1        0.32       96.20
   6.944773 |          1        0.32       96.52
          7 |          1        0.32       96.84
   7.616261 |          1        0.32       97.15
   8.533092 |          1        0.32       97.47
   9.963737 |          1        0.32       97.78
   10.79093 |          1        0.32       98.10
   10.79216 |          1        0.32       98.42
   11.96002 |          1        0.32       98.73
   13.04525 |          1        0.32       99.05
   15.67544 |          2        0.63       99.68
   35.04525 |          1        0.32      100.00
------------+-----------------------------------
      Total |        316      100.00

. sum base_pf_dens_gap2 end_pf_dens_gap2

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
base_pf_de~2 |       319    1.515789    6.027293        -10   50.04525
end_pf_den~2 |       316    .5430522    3.493964  -10.65975   35.04525

. gen base_gap_sign = base_pf_dens_gap / abs(base_pf_dens_gap)
(15 missing values generated)

. gen end_gap_sign = end_pf_dens_gap / abs(end_pf_dens_gap)
(23 missing values generated)

. gen base_gap_sign2 = base_pf_dens_gap2 / abs(base_pf_dens_gap2)
(18 missing values generated)

. gen end_gap_sign2 = end_pf_dens_gap2 / abs(end_pf_dens_gap2)
(25 missing values generated)

. replace base_gap_sign = 0 if abs(base_pf_dens_gap)<1
(160 real changes made)

. replace end_gap_sign = 0 if abs(end_pf_dens_gap)<1
(180 real changes made)

. replace base_gap_sign2 = 0 if abs(base_pf_dens_gap2)<1
(168 real changes made)

. replace end_gap_sign2 = 0 if abs(end_pf_dens_gap2)<1
(189 real changes made)

. tab base_gap_sign end_gap_sign

base_gap_s |           end_gap_sign
       ign |        -1          0          1 |     Total
-----------+---------------------------------+----------
        -1 |        35         27          7 |        69 
         0 |        16        130         14 |       160 
         1 |         7         23         55 |        85 
-----------+---------------------------------+----------
     Total |        58        180         76 |       314 


. tab base_gap_sign end_gap_sign if meana!=meanb | base_pf_RI1!=end_pf_RI1

base_gap_s |           end_gap_sign
       ign |        -1          0          1 |     Total
-----------+---------------------------------+----------
        -1 |        22         27          7 |        56 
         0 |        16         88         14 |       118 
         1 |         7         23         49 |        79 
-----------+---------------------------------+----------
     Total |        45        138         70 |       253 


. tab base_gap_sign2 end_gap_sign2

base_gap_s |          end_gap_sign2
      ign2 |        -1          0          1 |     Total
-----------+---------------------------------+----------
        -1 |        29         28          3 |        60 
         0 |        19        132         17 |       168 
         1 |         5         28         53 |        86 
-----------+---------------------------------+----------
     Total |        53        188         73 |       314 


. tab base_gap_sign2 end_gap_sign2 if mediana!=medianb | base_pf_RI1!=end_pf_RI1

base_gap_s |          end_gap_sign2
      ign2 |        -1          0          1 |     Total
-----------+---------------------------------+----------
        -1 |        16         28          3 |        47 
         0 |        19         90         17 |       126 
         1 |         5         28         47 |        80 
-----------+---------------------------------+----------
     Total |        40        146         67 |       253 


. restore

. 
. ************************
. ** PART FIVE: FIGURES **
. ************************
. 
.         * * * * *
.         * Fig 2 *
.         * * * * *
. 
. use Output\mastercleaned, clear

. 
. ** Create non-zero perception gaps for control group respondents
. replace gap = 1.39 - belief_food if randomscreen2==4 & food==1
(41 real changes made)

. replace gap = 1.96 - belief_spf if randomscreen2==4 & food==0 
(40 real changes made)

. 
. ** Collapse by deciles
. drop if abs(gap)>=20
(10 observations deleted)

. forval p = 10(10)90 {
  2.         egen p`p' = pctile(gap), p(`p') by(treatment question control)
  3.         }

.         
. ** Fit splines
. locpoly pf_revision1 gap if /* gap+3.25>=p10 & gap-3.25<=p90 & */ ri & control & spf, gen(xfit_CRcont yfit_CRcont) width(3.25) nograph adoonly

. locpoly pf_revision1 gap if /* gap+3.25>=p10 & gap-3.25<=p90 & */ ri & control & food, gen(xfit_FRcont yfit_FRcont) width(3.25) nograph adoonly

. locpoly pf_revision1 gap if /* gap+3.25>=p10 & gap-3.25<=p90 & */ ri & !control & spf, gen(xfit_CRinfo yfit_CRinfo) width(3.25) nograph adoonly

. locpoly pf_revision1 gap if /* gap+3.25>=p10 & gap-3.25<=p90 & */ ri & !control & food, gen(xfit_FRinfo yfit_FRinfo) width(3.25) nograph adoonly

. 
. ** Save spline fits
. preserve

. keep *fit* treatment question control

. gen id = _n

. reshape long xfit yfit, i(id) j(tmt) string
(note: j = _CRcont _CRinfo _FRcont _FRinfo)

Data                               wide   ->   long
-----------------------------------------------------------------------------
Number of obs.                      315   ->    1260
Number of variables                  12   ->       7
j variable (4 values)                     ->   tmt
xij variables:
xfit_CRcont xfit_CRinfo ... xfit_FRinfo   ->   xfit
yfit_CRcont yfit_CRinfo ... yfit_FRinfo   ->   yfit
-----------------------------------------------------------------------------

. replace control = .
(1260 real changes made, 1260 to missing)

. replace control = strpos(tmt,"cont")>0
(1260 real changes made)

. replace question = ""
(1260 real changes made)

. replace question = "RI" if strpos(tmt,"CR")>0 | strpos(tmt,"FR")>0
(1260 real changes made)

. * replace question = "PP" if strpos(tmt,"CP")>0 | strpos(tmt,"FP")>0
. replace treatment = ""
(1260 real changes made)

. replace treatment = "spf" if strpos(tmt,"CR")>0 /* | strpos(tmt,"CP")>0 */
(630 real changes made)

. replace treatment = "food" if strpos(tmt,"FR")>0 /* | strpos(tmt,"FP")>0 */
(630 real changes made)

. keep xfit yfit treatment question control

. save Output\Fig2splines, replace
file Output\Fig2splines.dta saved

. restore

. drop *fit*

. 
. ** Label deciles
. gen err_decile = .
(315 missing values generated)

. replace err_decile = 1 if gap < p10
(21 real changes made)

. replace err_decile = 2 if gap >= p10 & gap < p20
(17 real changes made)

. replace err_decile = 3 if gap >= p20 & gap < p30
(38 real changes made)

. replace err_decile = 4 if gap >= p30 & gap < p40
(14 real changes made)

. replace err_decile = 5 if gap >= p40 & gap < p50
(45 real changes made)

. replace err_decile = 6 if gap >= p50 & gap < p60
(4 real changes made)

. replace err_decile = 7 if gap >= p60 & gap < p70
(63 real changes made)

. replace err_decile = 8 if gap >= p70 & gap < p80
(23 real changes made)

. replace err_decile = 9 if gap >= p80 & gap < p90
(12 real changes made)

. replace err_decile = 10 if gap >= p90 & gap != .
(78 real changes made)

. keep gap pf_revision1 err_decile treatment question control

. collapse gap pf_revision1, by(err_decile treatment question control)

. 
. ** Create frequency weights for the scatterplots, for number of overlapping deciles
. egen decile_group = group(treatment question control)

. sort decile_group err_decile

. gen count_duples = err_decile[_n]-err_decile[_n-1] if decile_group[_n]==decile_group[_n-1]
(4 missing values generated)

. replace count_duples = err_decile[_n] if decile_group[_n]!=decile_group[_n-1]
(4 real changes made)

. 
. ** Create shading for quadrants 1 & 3, with different shading (i.e. different axes ranges) for each panel in the figure
. ** Note, the "expand 4" command doesn't change the appearance of the scatterplots, but is needed when we collapse by deciles
. expand 4
(81 observations created)

. sort treatment question

. by treatment question: gen counter=sum(1)

. gen y_shade = .
(108 missing values generated)

. gen x_shade = .
(108 missing values generated)

. egen y_shade_max = max(pf_revision1), by(treatment question)

. egen y_shade_min = min(pf_revision1), by(treatment question)

. egen x_shade_max = max(gap), by(treatment question)

. egen x_shade_min = min(gap), by(treatment question)

. replace y_shade = 5*floor(y_shade_min/5) if counter==1 | counter==2
(4 real changes made)

. replace y_shade = 5*ceil(y_shade_max/5) if counter==3 | counter==4
(4 real changes made)

. replace y_shade = 5 if y_shade < 5 & (counter==3 | counter==4)
(2 real changes made)

. replace x_shade = 5*floor(x_shade_min/5) if counter==1
(2 real changes made)

. replace x_shade = 0 if counter==2
(2 real changes made)

. replace x_shade = 0 if counter==3
(2 real changes made)

. replace x_shade = 5*ceil(x_shade_max/5) if counter==4
(2 real changes made)

. replace x_shade = 5 if x_shade < 5 & counter==4
(1 real change made)

. append using Output\Fig2splines

. egen xmax = max(x_shade), by(treatment question)

. egen xmin = min(x_shade), by(treatment question)

. egen ymax = max(y_shade), by(treatment question)

. egen ymin = min(y_shade), by(treatment question)

. replace y_shade = 6 if y_shade > 0 & (treatment!="spf" | question!="RI")
(676 real changes made)

. replace y_shade = -6 if y_shade < 0 & (treatment!="spf" | question!="RI")
(2 real changes made)

. 
. ** Make plots 
. twoway area y_shade x_shade if treatment=="food" & question=="RI", color(gs14) || line yfit xfit if treatment=="food" & question=="RI" & control==0 & yfit>=ymin & y
> fit<=ymax & xfit>=xmin & xfit <=xmax || scatter pf_revision1 gap [fw=count_duples] if treatment=="food" & question=="RI" & control==0, msize(small) msymbol(circle_h
> ollow) mcolor(gs2) ytitle("Revision") xtitle("Perception Gap") legend(off) yscale(range(-6(1)6)) ylabel(-5 0 5,angle(0))

. graph save Output\agg_updating_visFR_nocon, replace
(file Output\agg_updating_visFR_nocon.gph saved)

. twoway area y_shade x_shade if treatment=="spf" & question=="RI", color(gs14) || line yfit xfit if treatment=="spf" & question=="RI" & control==0 & yfit>=ymin & yfi
> t<=ymax & xfit>=xmin & xfit <=xmax  || scatter pf_revision1 gap [fw=count_duples] if treatment=="spf" & question=="RI" & control==0, msize(small) msymbol(circle_hol
> low) mcolor(gs2) ytitle("Revision") xtitle("Perception Gap") legend(off) ylabel(,angle(0)) /* yscale(range(-6(1)6)) ylabel(-5 0 5,angle(0)) */

. graph save Output\agg_updating_visCR_nocon, replace
(file Output\agg_updating_visCR_nocon.gph saved)

. twoway area y_shade x_shade if treatment=="food" & question=="RI", color(gs14) || line yfit xfit if treatment=="food" & question=="RI" & control==1 & yfit>=ymin & y
> fit<=ymax & xfit>=xmin & xfit <=xmax  || scatter pf_revision1 gap [fw=count_duples] if treatment=="food" & question=="RI" & control==1, msize(small) msymbol(circle_
> hollow) mcolor(gs2) ytitle("Revision") xtitle("Perception Gap") legend(off) yscale(range(-6(1)6)) ylabel(-5 0 5,angle(0))

. graph save Output\agg_updating_visFR_con, replace
(file Output\agg_updating_visFR_con.gph saved)

. twoway area y_shade x_shade if treatment=="spf" & question=="RI", color(gs14) || line yfit xfit if treatment=="spf" & question=="RI" & control==1 & yfit>=ymin & yfi
> t<=ymax & xfit>=xmin & xfit <=xmax  || scatter pf_revision1 gap [fw=count_duples] if treatment=="spf" & question=="RI" & control==1, msize(small) msymbol(circle_hol
> low) mcolor(gs2) ytitle("Revision") xtitle("Perception Gap") legend(off) ylabel(,angle(0)) /* yscale(range(-6(1)6)) ylabel(-5 0 5,angle(0)) */

. graph save Output\agg_updating_visCR_con, replace
(file Output\agg_updating_visCR_con.gph saved)

. 
. graph combine Output\agg_updating_visCR_nocon.gph Output\agg_updating_visCR_con.gph, rows(1) title("       SPF", size(medium))

. graph save Output\agg_updating_visCR_both, replace
(file Output\agg_updating_visCR_both.gph saved)

. graph export Output\agg_updating_visCR_both.emf, replace
(file H:\RESEARCH\PriceIsRight\Publication\Output\agg_updating_visCR_both.emf written in Enhanced Metafile format)

. 
. graph combine Output\agg_updating_visFR_nocon.gph Output\agg_updating_visFR_con.gph, rows(1) title("       Food", size(medium))

. graph save Output\agg_updating_visFR_both, replace
(file Output\agg_updating_visFR_both.gph saved)

. graph export Output\agg_updating_visFR_both.emf, replace
(file H:\RESEARCH\PriceIsRight\Publication\Output\agg_updating_visFR_both.emf written in Enhanced Metafile format)

. 
. graph combine Output\agg_updating_visCR_both.gph Output\agg_updating_visFR_both.gph, title({bf:   Treatment Group                Control Group}, size(medium)) rows(
> 2) ysize(20) xsize(15.5) imargin(0 0 0 0) t2title(" ", size(minuscule))

. graph export Output\Fig2.emf, replace
(file H:\RESEARCH\PriceIsRight\Publication\Output\Fig2.emf written in Enhanced Metafile format)

. 
. use Output\mastercleaned, clear

. 
.         * * * * *
.         * Fig 3 *
.         * * * * *
.         
. preserve

. 
. gen base_pf_RI1_flag0_5 = 0

. gen base_pf_RI1_flag5_10 = 0

. gen base_pf_RI1_flag10_15 = 0

. gen base_pf_RI1_flag15_50 = 0

. egen ri_N = total(ri==1 & control==0)

. 
. replace base_pf_RI1_flag0_5 = 100/ri_N if base_pf_RI1>=0 & base_pf_RI1<5
(204 real changes made)

. replace base_pf_RI1_flag5_10 = 100/ri_N if base_pf_RI1>=5 & base_pf_RI1<10
(60 real changes made)

. replace base_pf_RI1_flag10_15 = 100/ri_N if base_pf_RI1>=10 & base_pf_RI1<15
(29 real changes made)

. replace base_pf_RI1_flag15_50 = 100/ri_N if base_pf_RI1>=15 & base_pf_RI1<.
(32 real changes made)

. gen end_pf_RI1_cat = .
(325 missing values generated)

. replace end_pf_RI1_cat = 1 if end_pf_RI1>=0 & end_pf_RI1<5
(244 real changes made)

. replace end_pf_RI1_cat = 2 if end_pf_RI1>=5 & end_pf_RI1<10
(45 real changes made)

. replace end_pf_RI1_cat = 3 if end_pf_RI1>=10 & end_pf_RI1<15
(25 real changes made)

. replace end_pf_RI1_cat = 4 if end_pf_RI1>=15 & end_pf_RI1<.
(11 real changes made)

. label define pf_cats 1 "[0,5)" 2 "[5,10)" 3 "[10,15)" 4 "15+"

. label values end_pf_RI1_cat pf_cats

. 
. graph bar (sum) base_pf_RI1_flag* if ri==1 & control==0, over(end_pf_RI1_cat, relabel(1 "[0,5)" 2 "[5,10)" 3 "[10,15)" 4 "15+") sort(end_pf_RI1)) legend(col(1) pos(
> 3) subtitle("Baseline" "Expectations") order(4 3 2 1)) asyvars stack b1title("Final Expectations") l1title("Percent", angle(0)) bar(1, fcolor(gs14) lcolor(black) lw
> idth(vvthin) fintensity(90)) bar(2, fcolor(gs10) lcolor(black) lwidth(vvthin) fintensity(80)) bar(3, fcolor(gs7) lcolor(black) lwidth(vvthin) fintensity(80)) bar(4,
>  fcolor(gs4) lcolor(black) lwidth(vvthin) fintensity(90)) lintensity(0)

. graph export Output\updating_dist.emf, replace
(file H:\RESEARCH\PriceIsRight\Publication\Output\updating_dist.emf written in Enhanced Metafile format)

. 
. *Who moves where?
. tab end_pf_RI1_cat if base_pf_RI1_flag5_10 > 0

end_pf_RI1_ |
        cat |      Freq.     Percent        Cum.
------------+-----------------------------------
      [0,5) |         29       48.33       48.33
     [5,10) |         25       41.67       90.00
    [10,15) |          5        8.33       98.33
        15+ |          1        1.67      100.00
------------+-----------------------------------
      Total |         60      100.00

. **How much updating is there within each of the bars?
. gen non_rev = pf_revision1==0 if control==0 & pf_revision1!=.
(81 missing values generated)

. *Stayed in bar 1:
. tab non_rev if base_pf_RI1_flag0_5 > 0 & end_pf_RI1_cat == 1

    non_rev |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |         60       39.74       39.74
          1 |         91       60.26      100.00
------------+-----------------------------------
      Total |        151      100.00

. *Stayed in bar 2:
. tab non_rev if base_pf_RI1_flag5_10 > 0 & end_pf_RI1_cat == 2

    non_rev |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |          6       33.33       33.33
          1 |         12       66.67      100.00
------------+-----------------------------------
      Total |         18      100.00

. *Stayed in bar 3:
. tab non_rev if base_pf_RI1_flag10_15 > 0 & end_pf_RI1_cat == 3

    non_rev |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |          6      100.00      100.00
------------+-----------------------------------
      Total |          6      100.00

. *Stayed in bar 4:
. tab non_rev if base_pf_RI1_flag15_50 > 0 & end_pf_RI1_cat == 4

    non_rev |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |          4       44.44       44.44
          1 |          5       55.56      100.00
------------+-----------------------------------
      Total |          9      100.00

. 
. restore

. 
.         * * * * *
.         * Fig 4 *
.         * * * * *
. 
. preserve

. 
. ** Sample selection: treatment only, [-15,15] for perception gap, [-5,15] for baseline expectation
. keep if control==0
(81 observations deleted)

. keep if abs(gap) <=15
(21 observations deleted)

. keep if base_pf_RI1<=15 & base_pf_RI1>=-5 
(11 observations deleted)

. 
. ** "mywt" will become the density for the contour plots
. gen mywt = 1

. collapse (sum) mywt, by(gap base_pf_RI1 food spf)

. 
. ** Create "dummy" observations with zero density that extend the support of the joint density to the four corners of the plot
. count
   97

. local newobs = r(N)+8

. set obs `newobs'
obs was 97, now 105

. replace mywt = 0 if mywt==.
(8 real changes made)

. gen newobsindex = _n - (`newobs'-8) if mywt == 0 
(97 missing values generated)

. replace food = 0 if newobsindex>=1 & newobsindex<=4
(4 real changes made)

. replace spf = 1 if newobsindex>=1 & newobsindex<=4
(4 real changes made)

. replace food = 1 if newobsindex>=5 & newobsindex<=8
(4 real changes made)

. replace spf = 0 if newobsindex>=5 & newobsindex<=8
(4 real changes made)

. replace gap = -15 if newobsindex==1 | newobsindex==2 | newobsindex==5 | newobsindex==6
(4 real changes made)

. replace gap = 15 if newobsindex==3 | newobsindex==4 | newobsindex==7 | newobsindex==8
(4 real changes made)

. replace base_pf_RI1 = -5 if newobsindex==1 | newobsindex==3 | newobsindex==5 | newobsindex==7
(4 real changes made)

. replace base_pf_RI1 = 15 if newobsindex==2 | newobsindex==4 | newobsindex==6 | newobsindex==8
(4 real changes made)

. 
. ** Normalize the density to integrate to 1 in each treatment:
. egen wt_denom = total(mywt), by(spf)

. replace mywt = mywt / wt_denom
(97 real changes made)

. 
. ** Put the dummy densities just slightly below zero to give them a different color in the plot
. replace mywt = -.001 if mywt==-.1
(0 real changes made)

. drop newobsindex wt_denom

. label variable mywt "Density"

. 
. ** Finally, make the plots:
. ** There are occasional bugs with the Stata "contour" command so the following lines are commented out for stability:
. *twoway contour mywt gap base_pf_RI1 if spf, minmax yscale(range(-15(5)15)) ylabel(-15(5)15) xtitle("Baseline Expectation") ytitle("Perception Gap") ccuts(.0 .01 .0
> 2 .03 .04 .05 .06 .07 .08 .09 .10 .11 .12) ccolors(gs16 gs15 gs14 gs12 gs11 gs10 gs9 gs8 gs7 gs6 gs5 gs4 gs2 gs0) title("SPF Treatment")
. *graph save Output\contour_spf, replace
. *graph export Output\contour_spf.emf, replace
. *twoway contour mywt gap base_pf_RI1 if food, minmax yscale(range(-15(5)15)) ylabel(-15(5)15) xtitle("Baseline Expectation") ytitle("Perception Gap") ccuts(.0 .01 .
> 02 .03 .04 .05 .06 .07 .08 .09 .10 .11 .12) ccolors(gs16 gs15 gs14 gs12 gs11 gs10 gs9 gs8 gs7 gs6 gs5 gs4 gs2 gs0) title("Food Treatment")
. *graph save Output\contour_food, replace
. *graph export Output\contour_food.emf, replace
. 
. restore

. 
. 
end of do-file

. help contour

. help locpoly

. log close
      name:  <unnamed>
       log:  H:\RESEARCH\PriceIsRight\Publication\masterlog.log
  log type:  text
 closed on:  18 Aug 2014, 14:41:11
----------------------------------------------------------------------------------------------------------------------------------------------------------------------
